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Abstract 


An  experiment  to  measure  the  variation  of  wall  shear  stress  and  Reynolds  stress  over  a 
•  flat  plate  downstream  of  a  honeycomb  manipulator  was  conducted.  Two  velocity 
components  of  the  flow  were  measured  with  an  X-configured  hot-wire  probe,  and  wall 
shear  stress  was  obtained  from  a  surface  differential  pressure  gauge,  or  surface  "fencd*T 
Unmanipulated  boundary  layer  data  was  collected  at  eight  different  &ee  stream  speeds  to 
establish  boundary  layer  characteristics  and  to  calibrate  the  surface  differential  pressure 
gauge  by  making  use  of  Coles’  law  of  the  wall  for  die  mean  velocity  profile. 

Manipulated  boundary'  layer  data  was  obtained  at  three  free  stream  speeds  and  six 
manipulator  positions.  Almougn  manipulator  hole  size  and  length  remained  constant, 
height  was  varied  to  determine  the  influence  of  this  length  scale  upon  the  surface  friction 
coefficient. 


Manipulated  stream  and  vertical  mean  and  fluctuating  velocity  profiles  and  Reynolds 
stress  profiles  were  compared  to  their  unmanipulaied  counterparts  to  map  the  influence  of 
the  manipulator  in  the  streamwise  direction.  Immediately  dow  nstream  of  the 
manipulators  the  transport  of  turbulent  energy  to  the  wall  via  the  working  of  Reynolds 
stresses  was  sharply  reduced.  Complete  restoration  of  Reynolds  stress  profiles  to 
unmanipulated  form  was  not  observed  out  to  the  farthest  downstream  measurement 
position.  Local  drag  reduction  up  to  40%  was  verified  via  comparison  of  the 
unmanipulated  and  manipulated  friction  coefficients  obtained  from  the  wall  shear  stress 
measurements.  Mid-layer  peaking  of  the  Reynolds  stress  profiles  well  downstream  of  the 
manipulator  wai  accompanied  by  a  near-wall  depression  of  Reynolds  stress  levels  when 
compared  to  th^  unmanipulated  case.  A  decrease  in  the  magnitude  of  wall  shear  stress 
could  thereby  he  attributed  to  a  reduction  in  turbulent  energy  transport  toward  the  wall. 
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Nomenclature 


a,b 

^ ' f,m »  ^ ' fum 


q 


2 


U  V 


u 


u 

u 

or  u 


constant  coefficients  in  polynomial  or  other  equations 

manipulated  and  unmanipulated  wall  friction  coefficients  where 
Cf  =  2u2JUl 

sensor  diameter 

instantaneous  voltage  across  sensor 
fence  height  in  viscous  units  =  huj\ 
height  of  the  Boundary  Layer  Manipulator 
sensor  length 
exponent 

mean  pressure  at  a  given  position 
fluctuating  pressure  at  a  given  position 

rms  wall  pressure  fluctuation 
turbulent  kinetic  energy  =  «'2  +  v'2  +  w'2 
time  average  of  the  product  of  u’  and  v' 

instantaneous  velocity  in  stream  direction  at  height  y  above  wall 
(u  =«+«') 

mean  velocity  in  stream  direction  at  height  y  above  wall 
instantaneous  fluctuating  velocity  in  the  stream  direction 

rms  value  of  fluctuating  velocity  in  the  stream  direction 
non-dimensional  stream  velocity  -u!ux 


uz  friction  velocity  =  (tw/p): 

Um  or  UQ  velocity  in  stream  direction  outside  the  boundary  layer 
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H  > 


V 


^  or  v'm 

w 

w' 

y+ 

•*BLM 

8 

5' 

5‘und 

A P 

M- 

v 

P 


instantaneous  velocity  normal  to  the  wall  at  height  y  above  wall 
(v  =  v+v') 

mean  velocity  normal  to  wall  and  stream  velocity 
instantaneous  fluctuating  velocity  in  the  normal  direction 

rms  value  of  fluctuating  velocity  normal  to  the  wall  surface 

mean  velocity  parallel  to  the  plate  and  normal  to  the  stream  flow 

instantaneous  fluctuating  velocity  normal  to  the  plane  formed  by 
u'  and  v' 

yux 

non-dimensional  height  above  wall  =  — 

upstream  distance  of  the  Boundary  Layer  Manipulator  from  the 
measurement  position 

boundary  layer  thickness,  experimentally  defined  where 
u  =  .995 1/„, 

displacement  thickness,  for  incompressible  flow 

5'=IK> 

displacement  thickness  for  the  undisturbed  flow  for  a  given  free 
stream  speed  at  a  position  downstream  of  the  virtual  origin  where 
the  Boundary  Layer  Manipulator  is  placed 

differential  pressure  across  the  surface  fence 

mean  and  turbulent  viscous  dissipation 

absolute  viscosity  of  fluid 

kinematic  viscosity  of  fluid 

fluid  density 

total  shear  stress  =  p.(dM/9y)  -  pu  V 
x/p  =  v(du/dy)-u'v' 
wall  shear  stress 
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1  Introduction 


1.1  Background 

Both  the  reduction  of  skin  friction  drag  and  the  reduction  of  wall  pressure 
fluctuations  via  alteration  of  the  turbulent  flow  structures  within  the  boundary  layer 
surrounding  c  body  or  object  are  ongoing  efforts  in  the  fields  of  iluid  dynamics  and 
acoustics.  By  reducing  drag,  improved  speed  performance  and  fuel  economy  are 
achieved  for  aero-  and  hydro-vehicles.  By  minimizing  wall  pressure  fluctuations, 
improved  passive  SONAR  performance  is  achievable.  A  major  focus  of  the  research  has 
been  the  development  and  analysis  of  devices  to  break  up  the  large  eddy  structures  in  the 
turbulent  regime  which  are  believed  to  be  the  primary  sources  of  increased  drag  and  wall 
excitation.  These  devices  are  referred  to  as  LEBU’s  (Large  Eddy  Break-up  Devices)  or, 
more  recently,  as  Boundary  Layer  Manipulators  (herein  BLM’s).  The  latter  term  is 
preferred  since  it  does  not  presuppose  the  physical  mechanism  involved. 

To  date  inconsistent  results  have  been  obtained  concerning  the  effectiveness  of 
BLM’s.  Compounding  the  problem  are  the  various  geometries  of  the  BLM  devices  used. 
Beeler  [1]  found  a  12.5%  reduction  in  yjp'l,  using  two  rigid  ribbon  manipulators  (thin 
parallel  strips)  mounted  in  tandem  in  the  boundary  layer,  but  measurement  was  restricted 
to  only  one  position  downstream  of  the  BLM.  In  a  more  detailed  study  Nguyen  et  al  [2] 
used  single  and  triple  ribbons  to  measure  skin  friction  reductions  at  various  locations 
downstream  of  these  BLM’s  with  a  peak  reduction  of  35%.  A  departure  from  and  then 
re-establishment  of  Coles’  logarithmic  law  of  the  wall  [3]  as  measurements  were 
extended  downstream  of  the  manipulators  was  also  documented.  Hefner  et  al  [4]  used  a 
sawtooth  serrated  device  to  obtain  a  peak  skin  friction  reduction  of  24%,  and  also 
compared  the  boundary  layer  profiles  for  jj-  and  ~~  for  the  manipulated  and 
unmanipulated  cases  at  different  stations  downstream  of  the  BLM.  Honeycomb  BLM’s 
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were  utilized  by  Moller  and  Leehey  [5]  to  make  point  measurements  of  mean  wall  shear 
and  wall  pressure  spectra  at  various  b rations  downstream  of  the  manipulator.  Roth  and 
Leehey  [6]  varied  honeycomb  height  and  size  to  describe  their  influence  on  the  boundary 
layer  profiles  for  u  and  and  to  measure  in  the  streamwise  direction. 

All  of  the  preceding  investigations  were  conducted  in  wind  tunnel  facilities  with 
measurements  conducted  over  a  flat  plate  or  wall.  Flow  speeds  ranged  from 
approximately  2  -  35  meters  per  second.  As  the  data  accumulates,  researchers  will  focus 
their  attentions  on  determining  what  types  and  sizes  of  BLM’s  are  appropriate  to 
accomplish  a  particular  objective,  whether  for  overall  drag  reduction  on  a  vehicle  or  local 
wall  pressure  fluctuation  reduction  for  SONAR  applications.  One  critical  parameter  has 
received  little  attention,  the  Reynolds  stress  (-pit  V  for  incompressible  flow).  By 
measuring  Reynolds  stress  profiles  at  various  streamwise  locations  downstream  of  a 
BLM  and  comparing  these  to  the  unmanipulated  flow,  a  better  understanding  of  the 
mechanisms  involved  in  altering  the  turbulent  flow  structures  should  emerge,  particularly 
in  the  transport  of  turbulent  energy  within  the  boundary  layer.  This  appears  the  next 
logical  step  in  understanding  the  factors  determining  manipulated  turbulent  flow 
structures.  Westphal  has  reported  alterations  of  the  Reynolds  stress  profile  at  various 
stations  downstream  of  single  and  tandem  strip  or  plate  manipulators  [7].  Westphal 
found  large  reductions  in  Reyrolds  stresses  immediately  downstream  of  these 
manipulators,  followed  by  a  substantial  increase  in  these  stresses  in  the  mid-layer  farther 
downstream.  Westphal’s  findings  provide  a  basis  of  comparison  for  the  measurements 
undertaken  and  are  discussed  in  more  detail  later. 

Mathematical  representations  of  energy  transport  within  a  turbulent  boundary  layer 
have  been  forwarded  in  order  to  explain  measured  characteristics  of  the  boundary  layer. 
In  particular,  the  influence  of  shearing  stresses  in  the  sub-layers  has  been  examined. 
Townsend  has  proposed  energy  conservation  relationships  for  boundary  layers  in  channel 
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flow  and  over  flat  plates  which  include  the  energy  flow  resulting  from  the  working  of 

Reynolds  and  viscous  stresses  within  the  inner  and  outer  layers  [8].  In  Townsend’s 

representation,  the  gradient  normal  to  the  wall  surface,  of  the  product  of  the  local  mean 

flow  and  the  total  shear  stress  in  the  outer  layer  provides  an  energy  flow  into  the 

near-wall  or  viscous  sub-layer.  Equation  (1.1)  provides  symbolic  representation  of  this 

term  based  on  a  coordinate  system  where  x  is  in  the  stream  direction,  y  is  normal  to  the 

wall  surface  (y  =  0),  and  z  is  parallel  to  the  plate  surface  (cross  channel). 

rate  of  energy  transfer 
perunitmass 

In  wind  tunnel  experiments,  Klebanoff  found  that  approximately  85%  of  the  total  energy 
dissipation  occurs  in  the  near-wall  region  within  a  y+  not  exceeding  30  [9].  The 
dissipation  has  both  viscous  and  turbulent  components.  Klebanoff  concluded,  albeit  with 
an  incomplete  set  of  measurements,  that  turbulent  energy  production  and  dissipation 
peaked  in  the  near-wall  layer  and  that  pressure  forces  effected  an  inward  flow  of  energy 
toward  the  wall  [9].  In  Townsend’s  model  these  pressure  forces  arise  due  to  the  action  of 
the  shearing  stresses  (x)  represented  in  Equation  (1.1),  which  are  dominated  by  the 
Reynolds  stress  outside  the  viscous  sub-layer. 

The  presence  of  a  manipulator  in  the  turbulent  boundary  layer  will  influence  energy 
transport  mechanisms,  depending  on  the  size  and  geometry  of  the  BLM.  By  measuring 
the  Reynolds  stress  profile  in  the  manipulated  layer,  it  should  be  possible  to  determine 
what  type  of  influence  a  given  manipulator  has  in  altering  the  energy  transport  to  the 
wall.  If  the  manipulator  reduces  the  turbulent  energy  received  at  the  wall,  then  the 
root-mean-square  wall  pressure  or  wall  shear  stress  should  also  decrease.  Such  an  effect 
will  certainly  contribute  to  improvement  in  SONAR  self-noise  for  any  transducers 
mounted  on  the  wall  immediately  downstream  of  the  manipulator(s).  Whether  an  overall 
drag  reduction  is  realized  depends  primarily  on  the  magnitude  of  the  self-drag  of  the 
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manipulator.  The  combined  optimization  of  both  self-drag  and  self-noise  reduction  may 
not  be  achievable.  An  obvious  long-term  goal  in  the  study  of  manipulator  influence  on 
altering  energy  transport  processes  in  the  boundary  layer  is  to  determine  which  type  of 
manipulator  best  accomplishes  this  end.  Since  no  Reynolds  stress  profile  measurements 
downstream  of  a  honeycomb  manipulator  have  been  documented,  this  type  was  selected 
in  the  experiments  undertaken.  Measurement  of  the  x  and  y  velocity  components  at 
specific  and  numerous  locations  within  a  turbulent  boundary  layer  and  measurement  of 
wall  shear  stress  with  a  surface  differential  pressure  gauge  or  similar  device  will  provide 
the  necessary  data  to  establish  the  influence  of  the  manipulator  on  energy  transport  to  the 
wall. 

Hot-wire  anemometry  is  the  most  frequently  used  measurement  method  of  turbulent 
boundary  layer  velocity  profiles.  Small  single  (or  U-configuration)  and  X-configuration 
sensors  are  common  in  wind  tunnel  applications  with  unidirectional  mean  flow.  The 
X-probe  is  able  to  resolve  two  velocity  components  of  the  flow  at  the  position  of  the 
sensors.  Numerous  calibration  techniques  have  been  devised  for  both  U-  and  X-probes 
for  use  over  a  narrow  temperature  band  and  over  the  range  of  velocities  anticipated.  To 
ensure  the  highest  accuracy,  many  alternative  calibration  schemes  were  investigated. 
X-probe  calibration  schemes  can  be  grouped  into  two  major  categories,  calibration  based 
upon  heat  transfer  principles  and  flow  geometry  and  calibration  based  upon  geometric 
considerations  only.  Both  approaches  rely  upon  voltage  data  from  the  probe  sensors  to 
reflect  the  instantaneous  local  flow  conditions.  Boundary  layer  velocity  data  were 
obtained  using  an  X-probe  calibrated  by  a  geometric  procedure  refined  extensively  to 
ensure  the  highest  accuracy  of  measurement. 

Wall  shear  stress  was  measured  by  a  surface  differential  pressure  gauge  or  surface 
fence  and  followed  comprehensive  techniques  developed  by  Gur  and  Leehey  of  the 
Acoustics  and  Vibration  Laboratory  of  the  Massachusetts  Institute  of  Technology.  The 
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fence  was  calibrated  at  the  fence  position  by  obtaining  mean  velocity  profiles  in  the 
unmanipulated  flow  with  the  same  X-probe  used  to  obtain  the  manipulated  profile  data. 
By  numerically  fitting  the  profiles  to  Coles’  "law  of  the  wall"  data  [3],  the  value  of  uz 
and,  hence,  xw  for  specific  flow  conditions  was  determined.  The  A P  across  the  fence  then 
becomes  the  measurable  quantity  from  which  is  extracted. 

1.2  Objectives 

Objectives  of  the  present  research  fall  into  two  main  categories,  those  involving  the 
measurement  process  and  those  involving  the  measured  parameters.  Regarding  the 
measurement  process,  the  objectives  are  to: 

(1)  investigate  experimentally  various  calibration  procedures  in  X-probe  hot-wire 
anemometry 

(2)  refine  the  calibration  process  to  obtain  the  highest  measurement  accuracy  with  the 
greatest  computational  efficiency 

(3)  adapt  and/or  devise  hardware  and  software  components  to  support  the  data 
acquisition  required  by  the  refined  calibration  procedure. 

Regarding  the  parameters  to  measure,  wall  shear  stress  and  Reynolds  stress,  the 
objectives  are  to: 

(1)  measure  the  variation  of  wall  shear  stress  downstream  of  a  honeycomb  boundary 
layer  manipulator,  of  fixed  length  and  hole  size,  at  different  free  stream  velocities 

(2)  concurrently  with  wall  shear  stress  measurements,  measure  the  variation  in 
Reynolds  stress  profiles  downstream  of  the  manipulators 

(3)  compare  measured  parameters  with  results  of  previous  investigators  and  interpret 
the  findings 

(4)  evaluate  manipulator  influences  on  Reynolds  and  wall  shear  stresses  and  assess 
the  correspondence  between  the  influences  at  each  measurement  location. 
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2  Equipment  Description 


2.1  Facility  and  Equipment 

2.1.1  Wind  Tunnel 

All  data  were  collected  in  the  low-turbulence,  low-noise  wind  tunnel  in  the  MIT 
Acoustics  and  Vibration  Laboratory  illustrated  in  Figure  2.1.  Documentation  on  the 
design  and  construction  of  the  wind  tunnel  facility  may  be  found  in  Hanson  [10].  A 
General  Electric  Model  DC-300  blower  motor  draws  air  into  the  tunnel  test  section 
through  inlet  honeycomb  matrices  and  filter  screens.  A  General  Electric  adjustable  speed 
drive  in  combination  with  a  Red  Lion  Controls  Model  DT3D  controller  regulates  motor 
speed  to  within  ±  0.2  revolutions  per  minute  (RPM)  with  digital  display  accurate  to  0.1 
RPM. 

The  test  sections  are  inside  the  blockhouse  and  consist  of  a  closed  duct  and 
semi-open  jet  configuration  with  no  side  walls.  The  closed  section  extends 
approximately  nine  feet  from  the  inlet  end  of  the  blockhouse  while  the  semi-open  section 
is  4.25  feet  long.  The  zero  point  for  the  traversing  system  is  situated  3.3  feet  from  the 
inlet  end  of  the  blockhouse  and  ranges  over  90%  of  the  semi-open  section.  The  internal 
cross  section  of  the  duct  is  square,  fifteen  inches  on  a  side  (38  cm  x  38  cm).  The  duct 
portion  has  smooth,  masonite  side  walls  with  several  access  windows  mounted  along  its 
length.  Rectangular  plexiglass  test  plates  are  mounted  on  support  rails  to  form  the  lower 
boundary  in  the  semi-open  section.  The  upper  boundary  is  wooden  and  adapted  to 
support  the  probe  traversing  system  described  below.  The  plexiglass  plates  were  joined 
together  edge-to-edge  with  gaps  filled  by  putty.  The  putty  seams  were  then  sanded  with 
600  grit  paper  to  form  a  continuous  smooth  surface.  Boundary  layer  measurements  were 
made  above  this  surface. 
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General  Specifications 

Contraction  Ratio:  20  to 
Test  Section:  30  cm  x  3£ 


In  order  to  ensure  a  stable  boundary  layer  at  the  location  of  the  surface  fence  at  all 
free  stream  speeds,  a  seven  centimeter  (streamwise)  #36  grit  sandpaper  trip  was  attached 
to  the  lower  wall  near  the  entrance  to  the  tunnel  test  section.  The  trip  spanned  the  entire 
width  of  the  test  section.  Manipulators  were  placed  well  downstream  of  the  trip. 

2.1.2  Traverse 

Positioning  of  X-probes  in  the  boundary  layer  or  in  the  stream  is  accomplished  by  a 
traversing  system  specifically  designed  for  the  wind  tunnel  by  the  former  MIT 
Instrumentation  Laboratory  (Cambridge,  Massachusetts).  Two  separate  motor 
sub-systems  control  the  position  of  the  probe  vertically  above  the  test  plate  and 
horizontally  along  the  plate  in  the  stream  direction.  For  horizontal  control  a  motor  drives 
a  cable  system  in  forward  or  reverse  to  apply  force  to  the  traverse  cart.  The  cart’s  wheels 
ride  on  tracks  mounted  to  the  upper  boundary  of  the  test  sections.  Sensitivity  in  the 
stream  direction  is  0.05  inches. 

A  second  motor  mounted  on  the  cart  drives  a  rotating  disk  which  forces  a  threaded 
rod  vertically  upward  or  downward  with  0.0005  inch  sensitivity.  Attached  to  the  rod  is  a 
foil-shaped,  hollow  aluminum  arm  which  penetrates  into  the  tunnel  test  section  through  a 
padded  slot  in  the  upper  boundary  of  the  test  section.  A  dual  sensor  probe  support  rod 
(TSI  Model  1 155)  is  mounted  inside  the  aluminum  fairing  at  foil  maximum  diameter  and 
protrudes  at  the  lower  end  of  the  fairing  for  easy  probe  mounting.  Dual  coaxial  electrical 
leads  run  through  the  probe  support  from  two  BNC  connectors  to  a  dual  sensor  female 
connector  at  the  probe  end.  Figure  2.2  illustrates  the  traverse  cart  and  mounted  probe 
situated  in  the  tunnel  test  section. 
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Plexiglass  Plate  Surface  Fence 


Figure  2.2  Traverse  System  Schematic 

2.2  Velocity  Probes 

The  X-probe  used  in  the  experiments  was  aTSI  model  1243-T1.5  with 
platinum-coated  tungsten  hot  wire  sensors.  The  dual  sensor  cylindrical  probe  is  specially 
designed  for  boundary  layer  measurements  since  the  probe  body  and  sensor  support 
needles  form  a  "J"  shape.  The  design  provides  minimum  disturbance  to  the  flow  near 
solid  boundaries.  Probe  and  sensor  characteristics  are  summarized  in  Table  2.1  [11]. 

The  sensor  wires  have  enlarged  stems  at  both  ends  for  soldering  to  the  probe 
support  needles.  When  placed  in  the  flow,  one  sensor  slopes  at  45"  to  and  in  the  flow 
direction  while  the  second  sensor  is  aligned  perpendicular  to  the  first  such  that  it  too 
slopes  at  45°  to  and  in  the  flow  direction  in  the  opposite  sense.  The  two  sensors  lie  in 
parallel  planes  1.0  millimeter  (0.04  inches)  apart,  and  when  placed  in  the  flow,  their 
centers  are  the  same  height  above  a  horizontal  reference  plane.  Since  the  sensors  are 
spatially  separate,  use  of  an  X-probe  configured  in  the  described  manner  assumes  that  the 
flow  structure  at  a  given  height  above  the  lower  wall  surface  does  not  vary  in  the 
spanwise  direction. 
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Table  2.1  X-Probe  Characteristics 


Parameter 

Value 

Diameter  (D)  of  Sensing  Area,  pm 

(in) 

4 

(0.00015) 

Length  (L)  of  Sensing  Area,  mm 

(in) 

1.25 

(0.050) 

Aspect  Ratio  (L/D) 

312.5 

Distance  Between  Supports,  mm 

(in) 

1.5 

(0.06) 

Velocity  Range,  m/sec 
(ft/sec) 

0.15-200 

(0.5-600) 

Maximum  Ambient  Temperature  (°C) 

150 

Maximum  Operating  Temperature  (°C) 

300 

Recommended  Operating  Temperature  (°C) 

250 

Resistance  Temperature  Coefficient(12/°C) 

0.0042 

Upper  Frequency  Response  @100  m/sec,  kHz 
(Constant  Temperature  Mode) 

600 

Nominal  Ambient  Resistance,  (12) 

6 

Nominal  Operating  Resistance,  (12) 

10.8 

2.3  Pressure  Probes 
2.3.1  Pitot 

A  right-angle  Pitot  tube,  inserted  into  the  flow  near  the  center  of  the  tunnel  test 
section,  provided  a  pressure  signal  from  which  streamwise  velocity  was  calculated.  The 
X-probe  was  calibrated  using  the  velocities  obtained  from  the  Pitot.  The  Pitot  was 
connected  to  a  Validyne  DP15-TL  pressure  transducer  which  provided  a  proportional 
electrical  signal  to  a  Validyne  CD23  Pressure  Transducer.  Further  discussion  of  the 
Pitot’s  role  in  the  measurement  process  appears  in  Sections  2.5  and  3.1  below. 
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2.3.2  Surface  Differential  Pressure  Gauge 

Wall  shear  stress  measurements  were  made  with  a  surface  differential  pressure 
gauge  known  as  a  surface  fence.  The  fence  is  a  ridge  transverse  to  the  flow  protuding 
into  the  flow  approximately  0.002  inches.  This  height  resulted  in  an  h*  ranging  from 
approximately  1.8  at  the  lowest  speed  used  in  manipulated  measurements  to  a  value  of 
2.54  at  the  highest  speed.  These  values  of  h+  permitted  the  use  of  a  direct  calibration 
approach  as  outlined  below.  The  low  fence  height  resulted  in  negligible  disturbance  of 
the  mean  flow  since  the  fence  remained  within  the  buffer  region  of  the  boundary  layer. 
Also,  fence  aspect  ratio  (=  250),  length  to  height,  is  sufficient  to  ensure  negligible  fence 
tip  effects.  The  fence  illustrated  in  Figure  2.4  was  flush-mounted  into  the  plexiglass  plate 
in  a  fixed  position  one  foot  upstream  of  the  semi-open  tunnel  test  section.  Flow  over  the 
fence  results  in  a  pressure  drop  between  the  upstream  and  downstream  sides.  Two  taps 
on  either  side  of  the  fence  transmit  the  pressures  to  a  Validyne  CD23  Pressure 
Transducer.  The  proportional  signal  from  the  Validyne  enters  a  MASSCOMP  A/D 
converter  channel.  Software  transforms  the  input  signal  into  a  pressure  drop  in  Pascals. 

The  pressure  chamber  gaps  on  either  side  of  the  fence  are  0.002  inches  wide  by  0.12 
inches  in  length.  The  gap  length  in  viscous  units  ranged  from  108  at  the  lowest 
manipulated  measurement  speed  (15  m/sec)  to  150  at  the  highest  (23  m/sec).  Likewise, 
the  fence  length  of  0.5  inches,  in  viscous  units  ranged  from  450  to  630.  Other  relevant 
dimensions  are  provided  in  Figure  2.3. 

The  wall  shear  stress  at  the  surface  fence  can  be  determined  from  the  pressure  drop 
across  the  fence  [12].  Calibration  of  the  fence  was  accomplished  by  obtaining 
unmanipulated  mean  velocity  profiles  with  the  X-probe  over  a  range  of  stream  velocities 
and  numerically  fitting  each  to  Coles’  law  of  the  wall  relationship  to  obtain  the  friction 
velocity.  The  wall  shear  stress  for  each  profile  follows  directly  (Tw  =  pu<).  Leehey  and 
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Ports 


DIMENSIONS  (in) 

h  =  0.002 
w  =  0.002 

g  =  0.002 

d  =  1.26 
L  =  0.50 
L  =  0.12 


Figure  2.3  Surface  Differential  Pressure  Gauge 


Gur  used  this  approach  to  obtain  a  functional  relationship  between  differential  pressure 


across  the  fence  and  wall  shear  stress  given  by 


AP=a-C  (*+>  1) 


(2.1) 


where  the  exponent  n  is  approximately  1.5.  With  the  fence  calibrated,  wall  shear  stress 


under  manipulated  flow  conditions  was  measured.  Gur  and  Leehey  make  use  of  an 


alternative  non-dimensional  expression  for  the  response  of  the  surface  fence. 


(2.2) 


For  the  surface  fence  data  collected  during  the  boundary  layer  measurements,  Equation 


(2.2)  was  also  verified  to  be  consistent  with  (2.1)  and  with  the  findings  of  Gur  and 


Leehey. 
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2.4  Manipulators 

The  manipulators  chosen  were  the  honeycomb  variety  with  a  length  of  0.625  inches 
and  hole  size  of  0.125  inches.  In  viscous  unit  the  hole  size  {(T)  ranged  from  1 13  to  158. 
The  hole  size  and  length  were  based  on  the  findings  of  Roth  and  Leehey  who  found  these 
dimensions  to  reduce  wall  shear  stress  significantly  to  a  distance  of  at  least  300  5’^ 
downstream  of  the  manipulator  [6].  Since  the  differences  in  manipulator  heights 
considered  for  the  experiment  did  not  vary  widely,  no  alteration  of  the  length  scale  was 
justified.  The  length  selected  sufficiently  refined  the  flow  structure  such  that  the 
manipulator’s  influence  extended  well  downstream.  Two  views  of  a  manipulator  are 
provided  in  Figure  2.3.  The  manipulators  were  constructed  of  aluminum  with  a  web 
thickness  of  approximately  0.003  inches. 


Cross-Sectional  View 

Figure  2.4  Honeycomb  Manipulator  Schematic 


To  maintain  unidirectional  mean  flow,  the  manipulator  extended  the  full  width  of 
the  tunnel  section  and  was  positioned  such  that  the  side  walls  of  the  honeycomb  passages 
were  parallel  to  the  mean  flow.  The  manipulators  were  fixed  to  the  lower  tunnel  wall 
using  a  thin  layer  of  silicone  rubber  adhesive.  A  manipulator  of  a  specified  height  was 
positioned  upstream  of  the  surface  fence.  Surface  fence  and  X-probe  data  were  collected 
at  three  different  speeds  ((/„  =  15, 19,  and  23  m/sec).  The  manipulator  was  then 
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removed,  and  the  surface  cleaned  smooth.  An  identically  configured  manipulator  of 
different  height  (one  height  was  repeated)  was  then  attached  at  a  different  position 
upstream  of  the  fence  so  that  "h  =  8",  and  the  measurement  procedure  repeated.  The 
distance  between  the  manipulator  and  the  surface  fence  (or  the  stream  position  of  the 
profile  measurement)  is  represented  by  the  symbol  jcblm.  In  this  study  *BLM  had  values  of 
4, 12,  25,  65, 100,  and  150  centimeters. 

Manipulator  height  was  selected  to  approximate  the  undisturbed  boundary  layer 
thickness  based  on  its  position  relative  to  a  virtual  origin.  Indications  from  the  work  of 
Roth  and  Leehey  are  that  the  manipulator  effectiveness  is  optimized  under  this  condition 
[6],  For  flow  over  a  flat  plate,  the  relationship  governing  the  variation  of  turbulent 
displacement  thickness  8*  with  distance  x  downstream  from  a  virtual  origin  is  determined 
experimentally  as 

f  U _x  V5' 

8*(jc)  =  .046*  —  ,  (2.3) 

Ivy 

by  approximating  the  displacement  thickness  as  one-eighth  the  boundary  layer  thickness 
[13].  By  calculating  the  displacement  thickness  of  the  unmanipulated  boundary  layers  at 
a  given  free  stream  speed,  the  virtual  origin  of  the  flow  was  determined.  The  virtual 
origins  of  the  tripped,  unmanipulated  boundary  layers  at  15,  19,  and  23  meters  per  second 
were  respectively  3.52,  3.95,  and  4.67  meters  upstream  of  the  measurement  location  (the 
fence).  The  center  free  stream  speed  of  19  meters  per  second  was  the  reference  speed  for 
the  determination  of  8  where  each  BLM  was  located.  By  determining  the  position  of  a 
given  manipulator  relative  to  the  virtual  origin,  Equation  (2.3)  provided  8  (=  88‘). 
Manipulator  height  and  positions  are  summarized  in  Table  2.2. 
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Table  2.2  Manipulator  Positions  and  Heights 


XBLM 

(cm) 

^BLM 

(cm) 

Number  of  Holes 

4 

3.5 

11 

12 

5.1 

16 

25 

5.1 

16 

65 

4.1 

13 

100 

3.2 

10 

150 

3.8 

12 

The  trend  of  increasing  manipulator  height  as  the  measurement  position  from  the 
manipulator  increases,  as  dictated  by  Equation  (2.3),  was  intentionally  altered  at  the 
nearest  and  farthest  measurement  positions.  The  intent  was  to  augment  the  range  of  the 
ratio  of  manipulator  height  to  incident  boundary  layer  thickness  in  order  to  assess 
the  relationship  between  this  ratio  and  wall  shear  stress  (or  friction  coefficient).  A  more 
comprehensive  approach  would  set  manipulator  height  optimally  to  the  boundary  layer 
thickness  for  the  given  free  stream  speed  and  distance  from  the  virtual  origin.  These 
measurements  would  be  followed  by  two  more  sets  of  measurements,  one  with  a  constant 
height  manipulator  which  keeps  the  BLM  in  the  lower  50%  of  the  boundary  layer  at  all 
measurement  positions,  another  with  a  constant  height  BLM  always  exceeding  the 
boundary  layer  thickness  at  the  measurement  point.  In  this  way  a  broad  range  in 
hBLM/S'und  can  be  obtained. 

The  intent  in  using  a  honeycomb  manipulator  with  a  height  approximating  6  at  the 
manipulator  location  is  to  reduce  the  magnitude  of  the  fluctuating  velocity  components 
throughout  a  major  portion  of  the  boundary  layer  thickness.  Single  or  tandem  strip 
manipulators  should  be  less  effective  in  accomplishing  this  end.  The  honeycomb 
network  will  readily  remove  the  large  eddy  flow  structures  also.  The  wake  of  the 
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manipulator  will  initially  contain  small  scale  vortical  structures  at  different  levels  in  the 
boundary  layer  having  rotations  in  opposite  directions  in  the  x-y  plare  (one  direction  only 
shown  in  Figure  2.4)  [5].  Superposition  of  these  small  structures  immediately 
downstream  of  the  manipulator  should  effect  a  reduction  in  mean  and  fluctuating 
velocities  and  impede  the  transport  of  turbulent  energy  toward  the  wall.  Farther 
downstream  of  the  manipulator  the  small  scale  structures  will  merge  and  restore  a  larger 
scale  structure  to  the  boundary  layer  flow. 

2.5  Data  Acquisition  System 

A  MASSCOMP  Model  5400  32-bit  computer  is  the  core  of  the  data  acquisition 
system.  The  computer  system  contains  the  Model  AD12FA  Analog-to-Digital  (A/D) 
twelve  bit  converter  with  a  one  megahertz  throughput.  The  AD12FA  is  a  versatile  unit 
consisting  of  16  channels  with  an  input  voltage  range  from  -10  to  +10  volts.  Two  twelve 
bit  D/A  converters  provide  output  over  the  same  range.  Gain  is  programmable  across  the 
channels  separately. 

The  pressure  signal  from  the  Pitot  is  translated  into  a  proportional  electrical  signal 
by  a  Validyne  variable  reluctance  differential  pressure  transducer  (Model  DP15-TL)  in 
conjunction  with  a  Validyne  Model  CD23  Digital  Transducer  Indicator.  The  CD23 
produces  a  voltage  signal  equal  to  one-tenth  the  absolute  air  pressure  in  millimeters  of 
water.  A  voltage  divider  further  reduces  CD23  output  voltage  by  one-half  to 
accommodate  the  MASSCOMP  AD12FA  input  voltage  range.  The  Validyne 
components  are  regularly  calibrated  to  a  Betz  water  manometer. 

Since  the  X-probe  is  operated  as  a  constant  temperature  anemometer,  each  sensor 
forms  a  resistance  element  of  a  Wheatstone  Bridge.  Two  Dantec  Model  56C01/17 
Constant  Temperature  Anemometers  are  separately  balanced  to  each  sensor  so  that  each 
is  operated  at  the  desired  overheat  ratio  (OHR).  To  ensure  maximum  probe  sensitivity  an 
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OHR  of  1.9  was  selected  to  determine  bridge  settings.  This  overheat  value  corresponds 
to  an  operating  temperature  near  the  recommended.  Output  signals  from  the  bridges 
nominally  range  between  1. 0-4.0  volts  depending  upon  local  speed  and  the  type  and  size 
of  the  sensor. 

Sensor  frequency  response  was  evaluated  by  performing  a  Square  Wave  Test  with 
flow  incident  upon  the  probe  at  the  stream  speed  for  each  measurement.  A  square  wave 
generator  internal  to  the  Dantec  anemometers  provided  the  test  signals  with  response 
output  displayed  on  an  oscilloscope.  Measured  in  this  way,  each  sensor’s  response  was 
50  kHz. 


27 


3  Experiment  Description 

3.1  X-Probe  Calibration 

The  first  X-probe  calibration  procedures  applied  King’s  Law  of  Cooling  to 
determine  sensor  electrical  response  [14,15].  A  detailed  account  of  such  an  approach  was 
published  by  Schubauer  and  Klebanoff  in  work  performed  for  the  National  Advisory 
Committee  for  Aeronautics  (NACA)  [16].  To  obtain  high  accuracy  with  a  cooling  law 
calibration  approach,  precise  knowledge  of  each  sensor’s  alignment  to  the  flow  is 
required.  By  pitching  the  probe  relative  to  the  flow  for  different  calibration  speeds,  the 
alignment  of  each  sensor  to  the  flow  at  zero  pitch  angle  can  be  determined  [17,18].  Since 
other  calibration  approaches  exist  which  pitch  the  probe  during  the  calibration  procedure 
but  do  not  require  precise  knowledge  of  sensor  alignment,  a  cooling  law  approach  to 
calibration  offered  no  advantages.  The  overriding  consideration  in  rejecting  a  cooling 
law  calibration  approach  in  the  present  application  was  the  inherent  assumption  in  the 
procedure  that  v  is  negligible.  The  assumption  is  valid  for  unmanipulated  unidirectional 
mean  flow  but  not  for  manipulated  flows,  especially  near  the  manipulator. 

Recently,  X-probe  calibration  techniques  have  been  devised  which  eliminate 
assumptions  regarding  sensor  cooling  laws  and  specific  angular  alignment  between  the 
sensor  support  needles.  The  procedure  involves  pitching  the  probe  relative  to  the  flow  at 
each  calibration  speed  and  recording  the  speed,  angle,  and  sensor  voltages.  The  geometry 
for  the  procedure  is  illustrated  in  Figure  3.1.  The  probe  axis  bisects  the  projected 
intersection  of  the  two  sensors  such  that  the  free  stream  velocity  would  normally  parallel 
the  axis.  During  calibration  a  is  the  pitch  angle  of  the  probe  axis  relative  to  the  free 
stream.  During  boundary  layer  measurements  Q  represents  the  relative  velocity  with 
respect  to  the  probe  axis  at  any  instant  of  time  from  which  u  and  v  can  be  calculated. 
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Figure  3.1  X-Probe  Velocity  Geometry 

The  underlying  concept  behind  the  calibration  approach  is  that  a  unique  voltage  pair 
exists  across  the  sensors  for  a  specific  pair  of  velocity  components  u  and  v.  For  the 
concept  to  apply,  the  velocity  along  each  sensor  must  be  spatially  constant,  flow  must  lie 
in  the  plane  formed  by  the  two  sensors  (actually,  a  theoretical  plane  onto  which  both 
sensors  are  projected),  and  the  flow  vector  Q  lies  between  the  angle  formed  by  normals  to 
the  sensors,  nl  and  n2  (see  angle  (p.  Figure  3.1)  or  by  the  intersection  angle  of  the  sensors 
themselves,  whichever  is  smaller  [19].  Different  numerical  schemes  may  be  adopted  to 
structure  the  calibration  data  for  use  in  subsequent  boundary  layer  measurements. 

When  Willmarth  and  Bogar  used  the  velocity-angle-voltage  (VAV)  calibration 
approach,  they  computed  velocity  components  u  and  v  from  the  calibration  values  of  Q 
and  a  and  matched  these  to  corresponding  voltage  pairs  [19].  A  calibration  table  or  grid 
was  developed  separately  for  u  and  v  and  stored  for  later  use.  During  measurements  in 
the  boundary  layer  instantaneous  voltages  from  the  sensors  would  in  principle  fall  within 
a  "voltage  rectangle”  in  the  grid,  from  which  the  values  of  u  and  v  are  obtained  by 
interpolation.  Near  the  wall  (y+  <  65)  voltage  fluctuations  frequently  extended  beyond 
the  grid  limits  resulting  in  serious  errors  in  values  for  ^unw,  yjv'l,  and  u'v'  compared 
with  data  obtained  from  X-probes  calibrated  according  to  cooling  laws  [9].  Willmarth 
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and  Bogar  speculate  on  the  causes  and  significance  of  this  problem  and  conclude  that 
small  scale  turbulent  flow  structures  produce  different  flow  conditions  at  each  sensor 
since  the  sensors  are  spatially  separate  [19].  Thus,  although  the  calibration  voltage  grid 
encompasses  the  voltage  range  anticipated  for  the  intended  measurements,  velocities 
along  and  at  each  sensor  may  induce  a  voltage  response  uncharacteristic  of  the  calibration 
grid.  The  effect  may  be  compounded  by  velocity  components  in  the  third  dimension. 
Such  problems  are  inherent  in  the  assumptions  and  physical  characteristics  of  the  probes 
used  in  hot-wire  anemometry  and  would  invalidate  both  cooling  law  and  VAV  calibration 
approaches.  Since  other  factors  such  as  sensor  size  affect  probe  response  under  turbulent 
conditions,  all  must  be  considered  before  evaluating  the  appropriateness  of  a  given 
calibration  scheme.  For  example,  Willmarth  and  Bogar  used  sensors  much  smaller  than 
those  commercially  available  or  those  used  by  previous  researchers  such  as  Klebanoff,  a 
length  and  spacing  of  100  pm  and  diameter  of  0.5  pm  (see  Table  2.1  for  comparison). 
Lueptow  used  equivalently  small  probes  and  the  VAV  calibration  method  and  likewise 
documented  "out-of-grid"  voltages  near  the  wall  [20].  However,  in  comparing  their 
Reynolds  stress  measurements  near  the  wall  to  those  found  by  Klebanoff  and  others, 
Willmarth’s  and  Bogar’s  values  were  much  greater  while  those  of  Lueptow  were  much 
less. 

The  difficulties  encountered  by  Willmarth  and  Bogar  and  Lueptow  near  the  wall 
were  likewise  encountered  in  the  present  work.  The  poor  near-wall  data  obtained  with 
dual  sensors  can  be  attributed  to  sensor  spatial  separation,  flow  interference  due  to  sensor 
supports,  and  heat  transfer  to  the  cooler  wall.  Radiation  losses  appear  to  dominate  the 
heat  transfer  process  [14].  By  placing  an  X-probe  at  the  surface  in  a  no  flow  condition 
and  stepping  away  vertically  in  one -thousandth  inch  increments,  a  definite  voltage  or 
temperature  gradient  was  observed  when  using  both  a  standard  TSI  probe  and  a  TSI 
miniature  dual  sensor  probe  (2.5  pm  diameter,  0.5  mm  active  length).  The  temperature 
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gradient  became  flat  when  the  lower  supports  were  approximately  0.025  inches  from  the 
wall  or  when  the  probe  axis  was  approximately  .045  inches  from  the  wall.  Since 
radiative  heat  transfer  is  insensitive  to  flow  effects,  the  gradient  will  exist  under  flow 
conditions  and  introduce  voltage  errors  into  the  measurements.  The  errors  in  v  were 
determined  as  high  as  20%  when  the  voltage  error  was  only  3%.  This  is  due  both  to  the 
magnitude  of  v  (|  v  |  <  1.5m /sec)  and  the  slope  of  the  v-calibration  surface  for  near  zero 
velocities  (see  Figure  3.3  below).  These  errors  would  particularly  influence  the  values  of 

and  u  V.  No  suitable  correction  scheme  could  be  devised  to  circumvent  poor  probe 
performance  in  the  near  wall  region.  Thus,  the  near  wall  values  of  Reynolds  stress  could 
not  be  used  to  infer  and/or  confirm  the  value  of  friction  velocity  (or  wall  shear  stress)  for 
the  different  flow  conditions  encountered.  The  constant  stress  layer  approximation  given 
by  Equation  (3.1),  where  the  slope  of  the  mean  velocity  profile  becomes  negligibly  small 
in  comparison  to  the  Reynolds  stress  term,  failed  to  be  of  use  [21]. 

T^s-pwV  +  p^s-pMV  (3.1) 

dy 

Even  the  calculational  aspects  of  the  calibration  process  were  examined  in  an 
attempt  to  improve  accuracy  in  the  near  wall  region.  The  numerical  schemes  documented 
in  the  literature  vary  widely  in  complexity  and  computational  efficiency.  Willmarth  and 
Bogar  calculated  partial  derivatives  of  velocities  with  respect  to  voltages  at  nodes  of  the 
calibration  voltage  grid  in  order  to  interpolate  voltages  for  conversion  to  velocities  [19]. 
Johnson  and  Eckelmann  preferred  a  two-dimensional  Taylor  series  expansion  of  voltages 
to  obtain  first  and  second  order  partial  derivatives;  a  total  of  twelve  calibration  voltages 
and  partial  derivatives  were  required  for  each  point  of  the  calibration  grid  [22].  Lueptow 
et  al  constructed  velocity  look-up  tables  suitable  for  two  dimensional  linear  interpolation 
but  only  after  conducting  polynomial  and/or  spline  fits  relating  resultant  velocities  Q  and 
angles  a  to  voltage  pairs  in  a  finely  divided  voltage  grid  [23]. 
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Despite  the  evolution  of  the  aforementioned  numerical  techniques,  a  more  direct 


approach  was  adopted.  From  the  calibration  data  of  u,  v,  £0 ,  and  £0j  a  linear  least 
squares  fit  was  performed  to  determine  separately  the  coefficients  for  a  two  variable 
polynomial  in  E0  and  £q2,  one  polynomial  in  u,  the  other  in  v.  Oster  and  Wygnanski 
used  third-degree,  ten  coefficient  two-dimensional  polynomials  [24],  The  approach  used 
herein  expanded  upon  the  polynomial  strategy  by  raising  the  two-dimensional  polynomial 
to  degree  four  and  including  all  twenty-five  exponent  pair  combinations  of  £0  and  £0j. 
The  resulting  expressions  for  u  and  v  appear  in  equations  (3.2)  and  (3.3). 


u 


v 


=  i  i «*« 

j.Oi-O  1  2 


=  i  i 

jmOi-O  1  2 


(3.2) 

(3.3) 


The  twenty-five  coefficients  determined  for  each  velocity  component  were 
incorporated  into  a  conversion  routine  to  calculate  instantaneous  velocity  components  u 
and  v  from  voltages  obtained  from  the  sensors  during  boundary  layer  measurements.  By 
applying  this  procedure  separate  calibration  surfaces  for  u  and  v  were  constructed. 

Actual  calibration  surfaces  with  the  functional  form  of  Equations  *,3.2)  and  (3.3)  appear 
in  Figures  3.2  and  3.3,  respectively.  The  instantaneous  voltages  encountered  during 
boundary  layer  measurements  lie  within  the  calibration  voltage  range  of  each  sensor, 
thus,  the  two  surfaces  obtained  from  the  least  squares  fit  provide  all  possible  velocity 
components  measured  as  the  X-probe  traverses  the  boundary  layer.  The  u  and  v  velocity 
components  are  calculated  directly  from  the  appropriate  functions  vice  interpolated  from 
a  look-up  table.  Both  the  functional  and  Leuptow  look-up  table  methods  were  evaluated 
prior  to  boundary  layer  measurements  to  ascertain  the  accuracy  of  each.  The  functions 
specified  in  Equations  (3.2)  and  (3.3)  were  equivalently  accurate  to  the  look-up  table 
method  when  comparing  calculated  or  interpolated  velocities  to  original  calibration 
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velocities  for  a  given  voltage  pair.  In  terms  of  computation  time  during  data  acquisition, 
the  conversion  functions  proved  twice  as  fast  as  the  look-up  table  approach. 
Computational  speed  was  a  major  consideration  since  small  distance  steps  were  planned 
and  implemented  in  measuring  the  boundary  layer  profiles  to  obtain  a  detailed  picture  of 
small  scale  disturbances. 


Figure  3.2  Stream  Velocity  Calibration  Surface 
In  practice  the  calibration  procedure  requires  that  either  the  tunnel  speed  be  kept 
fixed  while  the  probe  is  pitched  at  various  angles  to  the  flow  or  that  the  probe  be  exposed 
to  the  full  range  of  calibration  speeds  at  a  given  angle,  which  is  repeatedly  re-set  to 
completion.  The  facilities  required  that  the  latter  procedure  be  followed.  Calibration 
angle  was  manually  re-set  by  positioning  the  traverse  cart  onto  a  step  support  positioned 
for  equal  positive  and  negative  calibration  angles.  Figure  3.4  depicts  the  physical 
arrangement  utilized  to  determine  the  calibrations  surfaces.  Seven  probe  pitch  angles  (O', 
±9",  ±18°,  ±36°)  were  set  by  using  the  step  device  illustrated.  With  a  given  angle  set,  the 
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Figure  3.3  Normal  Velocity  Calibration  Surface 
wind  tunnel  velocity  was  varied  from  zero  (to  establish  reference  Pitot  pressure)  to  a 
value  a  minimum  of  10%  above  the  maximum  free  stream  velocity  used  for  boundary 
layer  measurements.  Sufficient  time  was  allowed  for  the  flow  to  stabilize  before  the 
voltages  across  each  sensor  were  measured.  The  procedure  was  repeated  for  each  angle 
ensuring  that  the  temperature,  measured  at  the  wind  tunnel  settling  chamber  inlet  (see 
Figure  2.1),  remained  within  a  ±3°F  band.  By  adjusting  the  traverse  position,  the  probe 
remained  within  the  central  core  region  of  the  flow  away  from  the  tunnel  test  section 
walls.  The  final  calibration  measurement  was  taken  at  the  0°  position  with  the  traverse 
locked  in  position  for  boundary  layer  measurements.  Voltages  for  stream  speeds  less 
than  1.5  meters  per  second  were  discarded  before  determining  the  calibration  surfaces 
since  these  were  found  to  be  influenced  by  mutual  heat  transfer  between  the  sensors. 
Boundary  layer  measurements  near  the  wall  were  subsequently  limited  to  flow  conditions 
corresponding  to  a  local  mean  velocity  in  which  the  minimum  instantaneous  velocity  is 
1.5  meters  per  second.  Since  the  probe  axis  of  the  TSI  X-probe  used  in  this  study  cannot 
be  positioned  closer  than  0.375  millimeters  of  the  wall,  near-wall  measurements  are  valid 
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when  the  combination  of  free  stream  speed  and  streamwise  position  of  the  probe  relative 
to  the  virtual  origin  of  the  boundary  layer  result  in  no  instantaneous  local  speeds  below 
1 .5  meters  per  second. 


Traverse 


Plexiglass  Plate 

Figure  3.4  Wind  Tunnel  Calibration  Configuration 
After  discarding  the  calibration  data  corresponding  to  velocities  along  the  probe 
axis  less  than  1.5  meters  per  second,  two  sets  of  eighty-four  voltage-voltage-velocity 
triplets  remained  for  determination  of  the  calibration  surfaces.  Although  Johnson  and 
Eckelmann  recommend  fitting  such  surfaces  only  for  the  calibration  velocities 
corresponding  to  a  speed  range  to  which  the  probe  is  subjected  during  experimental 
conditions  [22],  the  coefficients  determined  were  fit  for  a  speed  range  of  2-35  meters  per 
second.  Justifications  for  this  approach  are  many.  First,  the  minimum  mean  local  speed 
experienced  by  the  probe  near  the  wall  was  approximately  six  meters  pier  second. 

Second,  in  attempting  linear  least  squares  fits  over  restricted  speed  ranges  with  third  and 
fourth  order  two  dimensional  polynomials  (sums  running  to  3  and  4  in  Equations  (3.7) 
and  (3.8)),  the  variance  between  measured  speeds  and  speeds  calculated  from  the  derived 
functions  remained  essentially  constant.  Third  and  most  importantly,  the  magnitude  of 
the  relative  error  between  measured  and  calculated  values  of  u  and  v  for  the  speed  range 
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encountered  in  the  boundary  layers  (5-23  meters  per  second)  was  typically  1%  or  less, 
with  some  higher  deviations  of  2%-4%,  primarily  in  r,  at  a  few  select  points,  particularly 
at  high  calibration  pitch  angles  where  the  influence  of  the  sensor  supports  begin  to  disturb 
the  flow  experienced  by  the  probe.  The  fourth  order  two-dimensional  polynomial 
functions  had  slighdy  lower  variances  and  relative  errors  than  the  third  order  functions 
and  were  used  in  the  calibration  procedure.  Restricting  the  fit  to  a  narrower  speed  range, 
however,  had  no  influence  on  the  relative  error  on  the  whole  in  the  speed  range  of 
interest.  An  individual  relative  error  between  a  measured  and  calculated  velocity  may 
increase  from  0.8%  to  1.1%  when  comparing  two  possible  calibration  ranges,  while 
another  may  decrease  from  -0.4%  to  -0.2%.  Thus,  with  no  advantage  in  selecting  a 
specific  calibration  range  based  on  the  speed  range  anticipated  during  boundary  layer 
measurements,  one  calibration  surface  each  for  u  and  v  applied  to  all  boundary  layer 
measurements.  Differences  in  u  between  measured  and  calculated  values,  using 
coefficients  from  the  least  squares  fit,  resulted  in  an  error  standard  deviation  equal  to 
0.044  meters  per  second.  The  error  standard  deviation  in  v  was  0.05 1  meters  per  second. 
Although  the  largest  magnitude  velocities  during  calibration  occur  in  the  stream 
direction,  due  to  the  large  pitch  angles  of  the  probe,  large  values  of  v  (up  to  23  meters  per 
second)  are  also  observed. 

3.2  Boundary  Layer  Measurements 

Once  the  X-probe  was  calibrated,  boundary  layer  measurements  proceeded.  The 
probe  sensors  were  positioned  just  forward  of  the  fence  ridge  with  the  probe  lowered  to 
the  surface.  The  first  step  was  to  obtain  velocity  profiles  for  the  unmanipulated  turbulent 
boundary  layer  at  the  speeds  of  interest.  Additional  unmanipulated  profiles  were 
obtained  at  speeds  above  23  meters  per  second  and  below  15  meters  per  second  to 
accomplish  a  proper  calibration  of  the  surface  fence.  The  data  acquisition  program  used 
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the  X-probe  calibration  coefficients  to  provide  an  output  of  the  following  quantities  at 
each  measurement  position  in  the  boundary  layer:  position  y,  mean  stream  velocity  u , 
mean  normal  velocity  v,  root  mean  square  fluctuating  stream  velocity  's/u'*,  root  mean 
square  fluctuating  normal  velocity  and  Reynolds  stress  mV.  A  total  of  150,000 
voltage  pairs  were  sampled  at  8000  Hertz  for  a  total  sampling  time  of  18.75  seconds  at 
each  of  some  150  positions  to  a  distance  of  approximately  eight  centimeters  from  the 
wall.  The  initial  position  of  the  probe  was  established  by  lowering  it  with  the  traverse 
such  that  the  lower  sensor  supports  were  just  in  contact  with  the  wall  surface.  The 
traverse  was  then  stepped  away  from  the  wall  0.001  inches  and  the  traverse  vertical 
position  indicator  zeroed.  The  estimated  distance  of  the  probe  axis  from  the  wall  was 
determined  by  extracting  a  "y-offset"  based  on  a  Coles  fitting  routine  to  be  discussed 
forthwith. 


Mean  stream  velocity  at  each  measurement  position  above  the  wall  surface  and  the 
corresponding  position  relative  to  the  arbitrary  zero  for  each  unmanipulated  profile 
formed  a  data  set  for  determining  the  friction  velocity  and  probe  offset  distance  from  the 
wall.  Experimental  data  was  fit  to  actual  data  obtained  by  Coles  for  an  equilibrium 
boundary  layer  [3].  An  iterative  least  squares  fitting  routine  solved  simultaneously  for 
friction  velocity  and  offset  distance.  Coles’  data  conforms  to  the  law  of  the  wall  as 
expressed  by  Equation  (3.4)  where  A  is  5.75  and  B  is  5.10  [3]. 

u*  =  AlogJy*)+B  (3.4) 


The  fit  was  accomplished  over  a  y+  ranging  from  20  to  200.  Figure  3.5  illustrates  the  law 
of  the  wall  fit  of  experimental  data  for  a  free  stream  speed  of  19  meters  per  second  for  an 
unmanipulated  boundary  layer;  other  speeds  are  found  in  Appendix  A.  With  uT  and  offset 
determined  for  each  profile,  a  calibration  curve  for  the  surface  fence  was  determined  (see 
Section  3.2.3)  and  boundary  layer  parameters  were  accurately  determined. 
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Unmanipulated  Mean  Velocity  Profile 


Figure  3.5  Law  of  the  Wall  Fit  (Unmanipulated,  U„=19  m/sec) 
Manipulated  profile  data  was  obtained  at  three  different  speeds  with  a  BLM  in  six 
different  positions.  All  data  acquisition  parameters  and  probe  positioning  routines  were 
identical  to  those  for  the  unmanipulated  profiles.  Surface  fence  differential  pressure  was 
converted  into  wall  shear  stress  and  friction  velocity  by  employing  the  fence  calibration 
curve.  In  order  to  determine  boundary  layer  parameters  accurately,  the  sensor  offset 
distance  from  the  wall  was  determined  by  using  a  modified  version  of  the  law  of  the  wall 
fitting  routine.  With  ur  known  for  each  stream  speed  a  value  for  the  y-offset  was  entered 
into  the  analysis  program  until  the  slope  of  mean  velocity  profile  in  the  range 
20  <  y  +  <  200  equaled  that  for  the  law  of  the  wall  profile.  The  value  of  y-offset  thus 
determined  was  added  to  each  raw  data  value  of  y,  in  appropriate  units,  to  establish 
relationships  for  calculating  boundary  layer  parameters  and  producing  accurate  profile 
plots.  Tables  3.1  through  3.4  summarize  the  major  boundary  layer  parameters  obtained 
for  the  three  free  stream  speeds  used  with  manipulators.  The  symbol  H  represents  the 
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shape  factor,  the  ratio  between  the  displacement  thickness  and  momentum  thickness  of 
the  boundary  layer,  and  is  presented  for  completeness  only.  Figures  3.6  (a)  through  (f) 
show  the  manipulated  boundary  layer  law  of  the  wall  fit  for  a  free  stream  speed  of  19 
meters  per  second  with  the  manipulator  at  each  of  its  six  distances  from  the  measurement 
position.  The  unusual  form  of  the  mean  velocity  profiles,  particularly  when  the 
measurements  are  near  the  BLM,  are  discussed  in  following  ^ub-section. 

Table  3.1  Unmanipulated  Boundary  Layer  Parameters 


u„ 

(m/sec) 

(m/sec) 

(Coles) 

(Pa) 

8* 

(10'3m) 

Ret 

Cf 

HO'3) 

H 

15 

0.55 

0.366 

8.08 

7700 

2.87 

1.43 

19 

0.67 

0.533 

8.69 

10408 

2.63 

1.44 

23 

0.77 

0.710 

9.35 

13478 

2.41 

1.43 

Table  3.2  Manipulated  Boundary  Layer  Parameters  (U„=15  m/sec) 


XBLM 

(cm) 

^BLM 

(cm) 

S’una 

(10‘3m) 

ut 

(m/sec) 

(fence) 

% 

(Pa) 

8' 

(10‘2m) 

Re&. 

Q 

(103) 

H 

4 

3.5 

7.99 

0.44 

0.233 

1.64 

15596 

1.82 

12 

5.1 

7.86 

0.41 

0.205 

1.75 

16484 

1.62 

25 

5.1 

7.60 

0.46 

0.259 

1.67 

15901 

2.03 

■Ml 

65 

4.1 

6.80 

0.47 

0.269 

1.59 

15707 

1.96 

1.62 

100 

3.2 

6.18 

0.51 

0.309 

1.28 

12148 

2.45 

1.52 

150 

3.8 

5.13 

0.49 

0.292 

1.59 

15654 

2.15 

1.56 
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Table  3.3  Manipulated  Boundary  Layer  Parameters  (U„=19  m/sec) 


XBLM 

(cm) 

^BLM 

(cm) 

S'und 

(10'3m) 

■ 

ut 

(m/sec) 

(fence) 

(Pa) 

5* 

(10'2m) 

Reb. 

cf 

(103) 

H 

4 

3.5 

8.36 

■sa 

0.342 

1.52 

18288 

1.65 

12 

5.1 

8.21 

mm 

0.324 

1.65 

20073 

1.54 

25 

5.1 

7.97 

mPM 

0.377 

1.57 

19263 

1.78 

1.49 

65 

4.1 

7.23 

0.60 

0.429 

1.39 

17611 

1.92 

1.57 

100 

3.2 

6.67 

0.62 

0.456 

1.19 

14327 

2.19 

1.51 

150 

3.8 

5.72 

0.60 

_ 

0.423 

1.55 

19130 

1.96 

1.55 

Table  3.4  Manipulated  Boundary  Layer  Parameters  (U„=23  m/sec) 


XBLM 

(cm) 

^BL-M 

(cm) 

5‘uno 

(10  3m) 

(m/sec) 

(fence) 

(Pa) 

5* 

(10'2m) 

Ret 

cf 

(10-3) 

H 

4 

3.5 

9.17 

0.64 

0.487 

■El 

22367 

1.53 

1.58 

12 

5.1 

9.03 

0.63 

0.466 

wBm 

25222 

1.45 

1.52 

25 

5.1 

8.84 

0.67 

0.537 

mm 

23225 

1.74 

1.50 

65 

4.1 

8.20 

0.69 

0.559 

Bitl 

20535 

1.79 

1.57 

100 

3.2 

7.64 

0.73 

0.628 

1.15 

16808 

2.08 

1.51 

150 

3.8 

6.79 

0.70 

0.583 

1.42 

20860 

1.90 

1.52 
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Manipulated  Mean  Velocity  Profile 


Figure  3.6(a)  Law  of  the  Wall  Fit  (hBLM=3.5  cm,  U«=19  m/sec) 
Manipulated  Mean  Velocity  Profile 


Figure  3.6(b)  Law  of  the  Wall  Fit  (hBLM=5.1  cm,  U„=19  m/sec) 


Manipulated  Mean  Velocity  Profile 


Figure  3.6(c)  Law  of  the  Wall  Fit  (hBLM=5.1  cm,  U„=19  m/sec) 


Manipulated  Mean  Velocity  Profile 


Figure  3.6(d)  Law  of  the  Wall  Fit  (hBLM=4.1  cm,  U_=19  m/sec) 


As  seen  in  Figures  3.6  (a)-(f),  Coles'  law  of  the  wall  does  not  apply  to  manipulated 
boundary  layers.  When  the  slope  matching  is  accomplished  to  determine  the  value  of  the 
y-offset,  the  mean  velocity  profile  for  a  manipulated  boundary  layer  parallels  the  Coles 
plot  for  an  unmanipulated  boundary  layer  for  values  of  y+  ranging  from  20  to  200  and 
above;  however,  the  manipulated  profile  has  a  values  of  u+  either  less  than  or 
approximately  equal  to  those  for  the  unmanipulated  profile.  This  contradicts  the  findings 
of  Roth  and  Leehey  where  the  manipulated  values  of  if  everywhere  exceeded  those  of 
the  unmanipulated  boundary  layer  [6].  To  determine  the  validity  of  the  present  findings, 
a  consistency  check  of  mean  stream  velocity  ratios  between  the  manipulated  and 
unmanipulated  cases  was  performed  for  a  given  y*  based  on  the  relative  magnitudes  of 
friction  velocity  in  each  case.  No  discrepancy  was  detected  nor  could  any  physical 
argument  be  found  to  justify  one  finding  in  preference  to  another. 

3.2.1  Unmanipulated  Velocity  Profiles 

Since  boundary  layer  measurements  were  taken  over  a  relatively  narrow  speed 
range,  graphical  representation  of  data  is  limited  to  measurements  in  the  boundary  layer 
with  a  free  stream  speed  of  19  meters  per  second.  Results  of  measurements  with  free 
stream  speeds  at  15  and  23  meters  per  second  closely  match  those  of  19  meters  per 
second.  For  completeness  Appendix  A  contains  relevant  plots  of  data  at  these  other  free 
stream  speeds. 

The  unmanipulated  mean  velocity  profiles  obtained  for  the  current  measurements 
adhere  in  form  and  proportion  to  the  extensively  documented  equilibrium  turbulent 
boundary  layer  profiles  measured  over  flat  plates  under  similar  conditions.  Figure  3.7  is 
one  such  profile  with  a  boundary  layer  thickness  of  approximately  six  centimeters  at  the 
point  of  measurement.  The  free  stream  speed  of  19  meters  per  second  is  a  rounded  figure 
based  on  the  averaging  of  all  free  stream  speeds  for  both  the  manipulated  and 
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unmanipulated  cases  taken  with  the  tunnel  fan  speed  set  to  a  specific  RPM 

(500  RPM  <=s>  18.6  ml  sec).  The  same  approximation  applies  to  the  free  stream  speeds 

of  15  and  23  meters  per  second. 


Figure  3.7(a)  Unmanipulated  Mean  Velocity  Profile  (U„=19  m/sec) 

To  specify  turbulent  flow  structure,  knowledge  of  the  root-mean-square  (rms) 
fluctuating  velocity  components  is  required.  For  equilibrium  turbulent  boundary  layers 
above  flat  plates  the  profile  data  obtained  by  Klebanoff  has  been  recognized  as  the 
standard  for  comparison  [16].  Klebanoff  obtained  his  data  with  single  and  double  sensor 
hot-wire  probes  appropriately  oriented  to  the  unidirectional  mean  flow  to  obtain  the  rms 
fluctuating  velocities.  In  comparing  the  rms  streamwise  fluctuating  velocity  profile 
obtained  for  unmanipulated  flow  to  that  of  Klebanoff  (Figure  3.7(b)),  a  close 
correspondence  exists  throughout  the  boundary  layer  except  near  the  wall.  This 
discrepancy  is  due  to  the  different  measuring  devices  employed.  Near  wall  hot-wire 
measurement  of  the  streamwise  fluctuating  velocity  component  can  only  be  effectively 
accomplished  with  a  single  sensor  probe  as  used  by  Klebanoff.  An  X-probe  of  the  type 
and  size  used  for  the  measurements  cannot  be  brought  close  enough  to  the  wall  to  resolve 
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the  near-wall  peak  in  Additionally,  the  physical  separation  of  X-probe  sensors  and 
heat  transfer  effects  between  the  wall  and  sensors  introduce  voltage  errors  which  corrupt 
velocity  measurements  in  the  near-wall  regime.  For  the  free  stream  speed  of  19  meters 
per  second  all  X-probe  data  is  considered  valid  when  y/8  attains  a  value  of  approximately 
0.03.  As  seen  in  Figure  3.7(b),  the  correspondence  between  yju^ obtained  with  the 
X-probe  and  that  obtained  by  Klebanoff  s  single  sensor  probe  begins  essentially  at  this 
point.  Thus,  the  unmanipulated  turbulent  boundary  layer  established  in  the  wind  tunnel 
conforms  parametrically  to  the  turbulent  boundary  layers  used  in  many  experiments 
designed  to  determine  turbulent  boundary  layer  characteristics. 


V 

Figure  3.7(b)  — vs  j  (Unmanipulated,  U„=19  m/sec) 


Similarly,  the  profile  for  compares  favorably  with  historical  data  (Figure 
3.7(c)).  Of  interest,  however,  is  the  large  near-wall  discrepancy  in  the  values  of  Vv'*. 
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Although  not  relevant  to  the  specific  objectives  set  forth,  the  near-wall  deviation  between 
the  current  measurement  of  > fv *  and  the  historical  data  is  discussed  to  reinforce  issues 
raised  previously  regarding  heat  transfer  effects  with  hot-wire  sensors  near  a  wall.  The 
TSI  X-probe  indicates  a  peaking  of  Vv*  similar  to  that  for  yfu^  obtained  by  Klebanoff 
and  others.  Use  of  a  TSI  X-probe  with  smaller  sensors  resulted  in  identical  profiles  for 
as  those  obtained  with  the  TSI  probe  having  larger  sensors,  as  in  Figure  3.7(c).  The 
near  wall  measurements  have  already  been  deemed  invalid  below  a  y/5  of  0.03,  but, 
unlike  the  data  comparison  for  a  dual  sensor  hot-wire  configuration  must  be  used  to 
resolve  Vv^.  Klebanoff  s  data  was  also  obtained  with  X-probes  and  was  subject  to 
similar  sources  of  error.  Reasons  for  the  near  wall  discrepancy  in  from  that  of 
Klebanoff  include  differences  in  wall  material  (plexiglass  versus  aluminum  plate)  and  in 
type  of  sensor  supports  (TSI’s  curved  bronze  prongs  versus  Klebanoff’ s  fine  jeweler’s 
broaches).  These  would  produce  different  local  heat  transfer  effects  which  ultimately  are 
reflected  as  velocity  errors. 

Voltage  measurements  of  the  X-probe  sensors,  as  the  probe  was  stepped  away  from 
the  plexiglass  wall  in  a  no-flow  condition,  indicated  an  exponential-like  drop  in  voltage 
from  the  wall  surface  to  a  distance  where  the  lower  sensor  supports  were  0.030  inches 
above  the  wall  surface;  the  voltages  of  each  sensor  stabilized  beyond  this  distance  and 
differed  from  the  wall  value  by  approximately  0.07  volts.  Such  a  procedure  was  not 
conducted  above  an  aluminum  surface  but  should  be  accomplished  for  comparison  and  to 
determine  the  validity  of  earlier  measurements.  The  small  voltage  changes  observed  in 
the  no-flow  condition  as  the  wall  is  approached  have  the  greatest  influence  on  the  value 
of  the  instantaneous  vertical  velocity  component,  if  the  effect  persists  under  flow 
conditions.  If  radiative  heat  transfer  is  the  dominant  mechanism  responsible  for  the 
voltage  changes  observed,  then  voltage  errors  are  introduced  into  the  flow  measurements 
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out  to  a  distance  from  the  wall  of  0.030  inches.  Radiation  from  a  hot-wire  probe  to  a 
nearby  surface  has  been  found  a  significant  heat  transfer  mechanism  [30].  This 
phenomenon  could  explain  the  sharp  rise  in  Vv*  near  the  wall. 


V 

Figure  3.7(c)  —  vs  j  (Unmanipulated,  U„=19  m/sec) 

In  summary,  the  unmanipulated  velocity  profiles  obtained  in  the  wind  tunnel  at  free 
stream  speeds  of  15, 19,  and  23  meters  per  second  indicate  the  existence  of  an 
equilibrium  turbulent  boundary  layer  parametrically  equivalent  to  the  turbulent  boundary 
layer  found  above  a  flat  plate  over  a  similar  range  of  Reynolds  numbers  based  on  length 
along  the  plate  in  the  stream  direction. 

3.2.2  Manipulated  Velocity  Profiles 

Interesting  phenomena  are  observed  in  comparing  the  manipulated  profiles  to  the 
unmanipulated  profile  as  both  height  and  upstream  position  of  the  BLM  are  varied. 
Figure  3.8(a)  shows  a  well-defined  mean  stream  flow  structure  exiting  the  downstream 
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end  of  the  BLM.  The  10-1 1  holes  from  the  bottom  to  top  of  the  manipulator  are  clearly 
discemable  as  the  profile  undergoes  rapid  oscillations  due  to  the  honeycomb  hole 
structure.  The  profile  is  noticeably  suppressed  in  magnitude  up  to  a  distance  from  the 
wall  approximately  that  of  the  manipulator  height.  The  same  trends  persist  farther 
downstream  from  the  BLM  as  shown  in  Figure  3.8(b)  except  that  the  flow  structure  has 
recovered  to  a  more  normal  form,  albeit  still  suppressed.  Also,  the  5.1  centimeter 
manipulator  height  is  revealed  by  the  form  of  the  profile.  Progressing  still  farther 
downstream  (i.e.  manipulator  is  positioned  farther  upstream  from  the  measurement 
position),  the  profile  acquires  a  different  form,  first  indicated  in  Figure  3.8(d)  where  xBLM 
equals  65  centimeters.  Stream  speed  appears  to  increase  linearly  with  height  above  the 
wall  from  a  position  very  close  to  the  wall  to  the  height  corresponding  to  the  approximate 
boundary  layer  thickness  where  u  ~  U„.  At  even  greater  distances  downstream  from  the 
BLM,  as  in  Figures  3.8(e)  and  3.8(f),  the  approximately  linear  profile  form  endures. 

A  noticeable  deviation  in  Um  between  the  manipulated  and  unmanipulated 

measurements  occurs  when  xBLM  is  65  centimeters.  Since  the  manipulator  both  increases 
the  drag  of  the  test  section  and  decreases  local  wall  friction  by  its  presence,  the  net  drag 
change  of  the  test  section  was  estimated  to  be  negligibly  small.  The  estimate  was  based 
on  the  work  of  Hanson  who  determined  that  the  test  section  only  contributed  16%  to  the 
overall  drag  of  the  wind  tunnel  [10].  Thus,  wind  tunnel  blower  motor  RPM  was  selected 
as  an  appropriate  means  of  establishing  the  mean  free  stream  velocity  in  the  test  section 
since  the  blower  would  not  experience  any  appreciable  loading  due  to  the  presence  of  the 
manipulator.  The  discrepancy  in  mean  free  stream  velocity,  then,  is  not  attributed  to  any 
influence  of  the  BLM  but  rather  to  acquiring  data  with  the  X-probe  when  the  largest 
acceptable  temperature  difference  (±3'F)  between  calibration  and  measurement  occurred. 
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Stream  Speed  (m/sec) 


Figure  3.8(a)  Velocity  Profile  Comparison  (U„=19  m/sec) 
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Figure  3.8(b)  Velocity  Profile  Comparison  (U„=19  m/sec) 
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Figure  3.8(c)  Velocity  Profile  Comparison  (U„=19  m/sec) 


xiLH«65  cm 


Figure  3.8(d)  Velocity  Profile  Comparison  (U„=19  m/sec) 
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Stream  Speed  (m/sec) 


Figure  3.8(e)  Velocity  Profile  Comparison  (U„=19  m/sec) 
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Figure  3.8(f)  Velocity  Profile  Comparison  (U«=19  m/sec) 
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A  significant  yet  unexpected  finding  during  the  boundary  layer  measurements  was 
the  influence  of  the  manipulator  on  the  mean  vertical  velocity.  Fortunately,  the 
calibration  technique  selected  for  the  X-probe  permitted  the  determination  of  this 
parameter.  In  an  unmanipulated  boundary  layer  v  «  u,  and  measurement  of  v  with  an 
X-probe  calibrated  with  the  velocity-angle-voltage  procedure  readily  confirms  this  fact. 

In  fact,  misalignment  of  the  probe  in  an  unmanipulated  flow  can  be  detected  by  a  value  of 
v  exceeding  approximately  1%  of  the  free  stream  speed.  Fluid  loading  on  the  probe 
sensors  and  supports  will  introduce  a  small  bias  into  the  measurement  resulting  in  slightly 
positive  values  of  v . 

A  typical  vertical  velocity  profile  for  an  unmanipulated  layer  appears  in  Figure 
3.9(e),  which  is  in  fact  a  measurement  for  the  manipulated  profile  150  centimeters 
downstream  of  the  manipulator  where  the  flow  has  essentially  recovered  equilibrium. 
Positive  velocity  indicates  velocity  away  from  and  normal  to  the  wall.  The  profile  is 
essentially  restored  to  its  unmanipulated  form,  a  straight  line  where  v  »  0.  (The  near-wall 
measurements  are  to  be  ignored  for  reasons  outlined  previously).  Moving  progressively 
towards  the  manipulator,  the  vertical  velocity  field  becomes  increasingly  disturbed  as 
indicated  in  Figures  3.9  (d),  (c),  (b),  and  (a).  Immediately  downstream  of  the 
manipulator  the  mean  vertical  velocity  field  oscillates.  Shown  well  by  Figure  3.9(a)  also 
is  a  small  positive  mean  vertical  flow  near  the  outer  edge  of  the  boundary  layer 
seemingly  originating  near  the  top  of  the  3.5  centimeter  manipulator.  A  similar  but  less 
pronounced  profile  pattern  occurs  above  the  5.1  centimeter  manipulator  further 
downstream  (Figure  3.9(b)).  A  large  streamwise  mean  velocity  component  must  be 
added  to  this  vertical  component  to  depict  the  resultant  mean  velocity  flow  structure. 
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Figure  3.9(d)  Vertical  Velocity  Profile  (U„=19  m/sec) 
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Vertical  Speed  (m/sec) 


Figure  3.9(e)  Vertical  Velocity  Profile  (U„=19  m/sec) 

Manipulated  boundary  layer  profiles  for  Vm'*,  normalized  with  respect  to  the  free 


stream  speed,  indicate  a  reduction  of  turbulent  intensity  in  the  near-wall  region  of  the 
boundary  layer  extending  outward  a  distance  approximately  the  height  of  the  individual 
BLM.  Very  near  the  BLM  the  strong  influence  of  the  manipulator  structure  is  evident  as 
yu*  varies  rapidly  with  position  (Figure  3.10(a)).  As  distance  downstream  from  the 
manipulator  increases,  the  rapid  variation  of  “V [u  *  dies  out,  and  its  value  is  some 
25%-40%  of  the  unmanipulated  profile  up  to  a  distance  from  the  wall  approximating  the 
manipulator  height  (Figures  3.10  (b)  and  (c)).  Farther  downstream  the  profile  is  restored 
in  the  outer  portion  of  the  boundary  layer,  but  near  the  wall  the  Vw'*  profile  remains 
essentially  flat  rather  than  increasing  as  the  wall  is  approached. 
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Figure  3.10(a)  -77-vs^  (U„=19  m/sec,xBLM=4  cm) 


*au<‘12cm 


Figure  3.10(b)  -^-vs^  (U.=19  m/sec, xBLM=12  cm) 


Figure  3.10(c)  vs  |  (U.=19  m/sec,xBLM=25  cm) 


*^-65  cm 


Figure  3.10(d) 


(IL=19  m/sec, xblm=65  cm) 
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Manipulated 

Unmanipulated 


Figure  3.10(e) 
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Vm  VS  8 


(U_=19  m/sec,xBLM=100  cm) 


Vu"150  cm 


Figure  3.10(f)  ^vsf  (U„=19  m/sec,xBUv(=150  cm) 
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The  effect  of  the  manipulators  on  the  rms  vertical  fluctuating  velocity  component 
follows  a  similar  pattern  to  that  of  although  the  profile  recovers  more  completely 
farther  downstream.  Near  the  manipulator  the  honeycomb  structure  dictates  the  profile  as 
evidenced  in  Figure  3.1 1(a).  The  oscillating  profile  pattern,  however,  significantly 
exceeds  the  unmanipulated  profile  values  of  yjv^.  Figure  3.1 1(a)  also  indicates  a 
suppression  of  the  near-wall  peak  in  Vv'*,  approaching  values  near  those  of  Klebanoff. 
Farther  downstream  of  the  manipulator  the  near-wall  peaking  in  Vv^  returns  and  persists 
(Figures  3.1 1  (b)-(f)).  At  intermediate  distances  downstream  of  the  manipulator  is 

reduced  in  magnitude  by  as  much  as  60%  and  remains  essentially  constant  from  a 
position  where  y/8  is  approximately  0.2  to  a  distance  above  the  wall  where  the  strong 
shearing  effects  at  the  top  of  the  manipulator  begin  to  influence  the  profile  (Figures  3.11 
(b)  and  (c)).  When  a  position  65  centimeters  downstream  of  a  BLM  is  reached,  the 
profile  in  Vv'*  is  partially  restored  to  its  unmanipulated  form  (Figure  3.1 1(d)).  The 
recovery  of  Jv*  (Figures  3.1 1  (e)  and  (f))  parallels  that  of  at  the  same  downstream 
measurement  positions  (Figures  10  (e)  and  (f))  in  that  values  of  >/v'*  are  lower  in  the 
manipulated  case  out  to  a  distance  from  the  wall  of  0.3y/8. 


V 

Figure  3.11(a)  —  vsj  (LL=19  m/sec,xBLM=4  cm) 


Figure  3.11(b)  ^  vs  j  (U«=19  m/sec, xBLM=12  cm) 
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3.3  Stress  Measurements 
3.3.1  Wall  Shear  Stress 

Wall  shear  stress  for  the  manipulated  boundary  layers  was  measured  directly  by  the 
calibrated  surface  differential  pressure  gauge.  Calibration  was  based  on  eight  data  points 
determined  from  boundary  layer  profile  data  for  free  stream  speeds  ranging  from  7-35 
meters  per  second.  By  applying  the  Coles  law  of  the  wall  fitting  routine  for  each  of  eight 
mean  stream  velocity  profiles,  friction  velocity  and  wall  shear  stress  were  determined. 
Wall  shear  stress  was  then  related  to  the  corresponding  pressure  drop  across  the  gauge 
measured  at  the  fence  for  each  measurement  speed.  The  calibration  was  accomplished 
over  a  range  of  h+  from  1.02  to  3.74.  Figures  3.12  (a)  and  (b)  provide  two  separate 
representations  of  the  calibration  relationship.  A  linear  least  squares  fit  of  the  log10(xw) 
and  log10(AP)  calibration  data  resulted  in  a  calibration  curve  represented  by  Equation 
(3.10)  and  plotted  in  Figure  3.12(a). 

xw  =  .28A Poa  (3.10) 

Equation  (3.10)  for  xw  agrees  well  with  the  findings  of  Gur  and  Leehey  using  a  similar 
gauge.  An  alternate  formulation  of  the  relationship  between  wall  shear  stress  and 
pressure  drop  across  the  surface  gauge  in  terms  of  h+  is  provided  by  Figure  3.12(b).  For  a 
given  fence  height  and  fluid,  the  formulations  are  equivalent. 
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3.3.2  Reynolds  Stress 

The  most  significant  influence  of  the  honeycomb  manipulators  appears  in  the 
Reynolds  stress  profiles  normalized  to  the  square  of  the  friction  velocity  (Figures  3.13  (a) 
-  (f)).  The  friction  velocity  was  obtained  from  surface  fence  measurements  at  the  same 
position  downstream  of  the  BLM  where  each  profile  was  measured.  The  profile  four 
centimeters  downstream  of  the  manipulator  (Figure  3.13(a))  dramatically  shows  the 
Reynolds  stress  term  obtaining  a  negative  value  (i.e.  u  V  >  0),  a  phenomenon  never 
occurring  within  an  unmanipulated  equilibrium  turbulent  boundary  layer.  Two 
phenomena  characterize  the  peaking  pattern  exiting  the  manipulator.  First  is  the  distinct 
influence  of  each  row  of  manipulator  holes  on  the  local  Reynolds  stress,  specifically  the 
oscillating  patten  observed  in  the  velocity  profiles.  Second  is  the  combined  influence  of 
all  rows  of  holes,  such  that  superimposed  upon  the  oscillating  Reynolds  stress  pattern  is  a 
curvature  resulting  in  the  largest  positive  values  of  «  V  at  approximately  the  mid-height 
of  the  manipulator  (shown  by  the  large  negative  peaks  in  near  the  manipulator 

“t 

mid-height  in  Figure  3.13(a)).  The  negative  peaks  in  —■  tend  to  decrease  in  magnitude  at 
measurement  positions  away  from  the  mid-height  toward  the  top  and  base  of  the  BLM. 

Farther  downstream  the  Reynolds  stress  becomes  positive  again  (i.e.  «  V  <  0)  but  is 

greatly  reduced  in  magnitude  over  a  depth  of  the  boundary  layer  corresponding  to  a  major 
fraction  of  the  manipulator  height  (Figures  3.13  (b)  and  (c)).  (Westphal  apparently  found 
a  similar  reduction  in  Reynolds  stresses  near  the  flat  plate  manipulators  used  in  his 
experiments  but  provided  no  graphical  representations  of  such  [7].)  The  individual 
manifestations  of  stresses  induced  by  flow  through  the  honeycomb  structure  have  merged 
due  to  diffusion  processes.  Likewise,  the  diffusion  of  the  strong  shearing  influence  of  the 
flow  over  the  top  boundary  of  the  manipulator  is  clearly  evident  as  the  fractional  width  of 
this  region  with  respect  to  the  total  boundary  layer  height  has  increased. 
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Still  faither  downstream  the  Reynolds  stress  profile  shows  signs  of  restoration  to  an 


unmanipulated  boundary  layer  form  but  takes  on  two  distinct  characteristics.  The  first  is 
a  minor  suppression  near  the  wall;  the  second  is  a  mid-layer  peaking  the  same  distances 
from  the  wall  where  farther  upstream  the  profile  was  greatly  suppressed  (Figures  3.13  (d) 
-  (f)).  As  mentioned  previously,  Westphal  found  similar  mid-layer  peaking  far 
downstream  from  single  and  tandem  flat  plate  manipulators  [7].  Westphal’s  profiles  lack 
the  fine  point-to-point  measurement  detail  of  the  present  work,  and  suppression  of  the 
Reynolds  stress  as  the  wall  is  approached  is  not  evident  at  his  downstream  measurement 
locations. 


Figure  3.13(a)  -t-vsj  (LL=19  m/sec) 
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Figure  3.13(e)  ^  vs  j  (U 
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,=19  m/sec) 
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4  Discussion  of  Results 


4.1  Manipulator  Influence  on  Velocity  Profiles 

In  Section  3.2.2  the  manipulated  mean  and  fluctuating  velocity  profiles  were 
individually  described.  In  this  section  a  more  comprehensive  view  of  the  data  is 
forwarded  which  integrates  these  individual  results.  The  mean  velocity  profiles  in  the 
stream  direction  generally  adhere  to  expectations  based  on  the  use  of  a 
honeycomb-structured  manipulator  (see  Section  2.4).  All  near-manipulator  profiles 
indicate  that  energy  has  been  extracted  from  the  mean  flow  and  that  the  scale  of  the  flow 
has  been  reduced.  Farther  downstream,  dispersion  results  in  a  the  loss  of  local  flow 
identity,  and  the  mean  flow  energy  is  gradually  restored  due  to  diffusion  from  the  outer 
regions  of  the  boundary  layer  to  the  inner.  Full  restoration  of  the  mean  velocity  profile  in 
the  stream  direction  was  not  observed  at  the  farthest  measurement  position  (xBLM=150 
cm).  The  mean  vertical  velocity  profiles  are  best  explained  by  considering  the  influence 
of  the  manipulator  channels  on  the  flow  at  different  channel  levels  of  the  manipulator.  A 
horizontal  flow  in  principle  should  exit  these  channels  but  due  to  the  different  incident 
velocities  at  the  different  channel  levels  and  due  to  minor  manufacturing  differences  of 
the  channels  an  irregular  vertical  pressure  distribution  will  exist  at  the  manipulator 
channel  exits.  The  pressure  differences  will  introduce  a  small  vertical  velocity 
component  onto  the  mean  flow  which  quickly  dissipates  farther  downstream. 

The  same  general  trends  discussed  for  the  mean  velocity  profiles  applies  to  the 
fluctuating  components.  Strong  variations  in  the  magnitude  of  the  rms  fluctuating  stream 
and  vertical  velocities  near  the  manipulator,  indicative  of  the  honeycomb  cell  structure, 
diffuse  and  superpose  to  reduce  sharply  the  magnitude  of  the  these  velocities  throughout 
a  major  fraction  of  the  boundary  layer  thickness.  The  mechanisms  responsible  for  these 
effects  and  the  partial  recovery  of  the  profile  to  an  unmanipulated  form  are  explored  in 
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the  next  sub-section.  It  is  more  instructive  to  examine  the  results  of  the  present 
measurements  in  comparison  to  similar  measurements  with  honeycomb  manipulators 
placed  in  a  free  stream.  Particular  attention  is  given  to  the  fluctuating  stream  component. 
The  works  of  Lumley  and  McMahon  [26],  Loerhke  and  Nagib  [27],  and  Batchelor  [28] 
provide  the  bases  of  comparison. 

Due  to  the  great  differences  in  the  flow  structure  incident  upon  a  honeycomb 
manipulator  occupying  the  cross-section  of  a  flow  channel  from  that  of  a  manipulator 
placed  within  a  boundary  layer,  direct  comparison  of  mean  and  fluctuating  velocity 
profiles  is  not  justified.  Streamwise  variation  of  these  velocities  is  the  primary  focus 
when  manipulating  the  free  stream  while  both  transverse  and  streamwise  variations  are  of 
interest  in  manipulated  boundary  layers.  What  is  of  interest  are  the  turbulent  energy 
decay  relationships  for  screens  and  grids,  which  have  also  been  applied  to  free  stream 
honeycomb  manipulators.  Batchelor  specifies  an  initial  decay  relationship  for  the 
turbulence  exiting  a  free  stream  manipulator  or  grid  such  that  un<*(x  -Xg)~\  where  Xg  is  a 
virtual  origin  approximately  equal  to  ten  times  the  honeycomb  cell  diameter  [28].  In 
general,  the  relationship  applies  for  approximately  100  cell  diameters  downstream  of  the 
manipulator,  which  is  approximately  32  centimeters  for  the  present  measurements.  With 
measurements  in  the  manipulated  boundary  layer  at  4, 12,  and  25  centimeters 
downstream  of  the  BLM,  no  such  decay  behavior  could  be  found.  Manipulated  boundary 
layer  energ>  decay  in  u'2  occurs  more  gradually  and  levels  off  in  the  region  12 
centimeters  downstream  of  the  manipulator.  If  the  peak  value  of  Reynolds  stress  at  four 
centimeters  from  the  manipulator  is  compared  to  the  mid-layer  value  at  12  centimeters, 
the  decay  is  comparable  to  that  found  for  grids,  but  between  12  and  25  centimeters,  the 
relationship  is  invalid.  Also,  the  turbulent  energy  decay  still  farther  downstream  of  grids 
varies  as  x5n  [26].  This  is  not  observed  for  boundary  layer  honeycomb  manipulators,  and 
the  turbulence  levels  rise  gradually  to  levels  observed  in  unmanipulated  boundary  layers 
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at  the  far  downstream  locations.  In  summary,  turbulent  energy  decay  relationships  and 
velocity  variations  in  the  stream  direction  for  grids  and  free  stream  honeycomb 
manipulators  have  little  bearing  on  the  present  investigation  of  turbulent  boundary  layer 
honeycomb  manipulators,  particularly  in  understanding  energy  transport  phenomenon. 
4.2  Manipulator  Influence  on  Measured  Stresses 

The  friction  coefficient  (Cf)  provides  a  measure  of  the  wall  shear  stress.  By  relating 
the  friction  coefficient  in  the  manipulated  boundary  layer,  at  a  specific  free  stream  speed 
and  location  relative  to  the  boundary  layer’s  virtual  origin,  to  the  friction  coefficient  of 
the  unmanipulated  boundary  layer  at  the  same  free  stream  speed  and  location,  the 
influence  of  the  manipulator  on  wall  shear  stress  is  quantified.  As  the  manipulator’s 
position  upstream  of  the  measurement  location  is  varied,  the  streamwise  influence  of  the 
manipulator  upon  wall  shear  stress  is  determined.  Figure  4.1  summarizes  manipulator 
influence  on  wall  shear  stress  as  a  result  of  varying  not  only  free  stream  speed  and 
manipulator  location  in  the  boundary  layer,  but  also  manipulator  height.  The  ratio 
between  the  manipulated  and  unmanipulated  friction  coefficients  is  plotted  as  BLM 
height  and  upstream  distance  from  the  measurement  location  are  varied,  both  of  which 
are  normalized  to  the  displacement  thickness  of  the  unmanipulated  boundary  layer  at  the 
location  of  the  BLM.  In  Figure  4.1  d*  represents  5*und,  and  xBLM  has  been  multiplied  by 
0.01  to  improve  overall  plot  scaling. 

As  shown  by  Figure  4.1,  a  35%-40%  reduction  in  the  friction  coefficient  is  evident 
just  downstream  of  the  manipulator;  this  reduction  is  relatively  insensitive  to  changes  in 
manipulator  height.  As  distance  downstream  from  the  manipulator  increases,  the  friction 
coefficient  ratio  becomes  relatively  constant  beyond  150  displacement  thicknesses  for  the 

^blm 

approximate  range  4  <  -7-  <  5.  Roth  and  Leehey  found  the  same  trend  in  the  behavior  of 

sw 

the  friction  coefficient  ratio  for  a  one-eighth-inch-cell  honeycomb  manipulator 
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Figure  4.1  Friction  Coefficient  Ratio 

approximately  1.3  inches  in  height  [6],  As  the  height-to-displacement  thickness  of  the 
manipulator  increases,  the  friction  coefficient  ratio  slopes  upward  toward  a  value  of  unity 
at  far  downstream  locations.  A  depression  in  the  friction  coefficient  ratio  is  observed  in 
the  range  6.0  <  <  6.5  out  to  a  distance  downstream  of  100  b*und.  By  referring  to 

6»u/ 

Tables  3.1  through  3.4,  the  depression  is  determined  to  occur  when  the  free  stream  speed 
is  19  meters  per  second  and  when  25cm  <  xBLM  <  65 cm  and  4cm  <  hBLM  <  5cm .  The  small 
number  of  data  points  used  in  developing  Figure  4.1  makes  discussion  of  the  causes  of 
this  surface  feature  premature;  it  is  also  premature  to  draw  conclusions  regarding 
optimization  of  parameters.  The  matter  warrants  further  detailed  investigation. 

By  examining  the  Reynolds  stress  profiles  at  all  measurement  speeds  and  the 
friction  coefficient  ratio  plot  at  measurement  positions  close  to  the  manipulator,  a 
reduction  in  the  magnitude  of  the  friction  coefficient  between  the  manipulated  and 
unmanipulated  boundary  layers  occurs  in  conjunction  with  a  reduction  of  the  Reynolds 
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stresses  in  the  boundary  layer,  excepting  the  large  shearing  peak  from  the  top  of  the 
manipulator.  Simply  stated,  closest  to  the  manipulator  the  Reynolds  stresses  are  least  and 
so  is  the  wall  shear  stress.  It  may  be  surmised  then  that  a  reduction  in  the  Reynolds  stress 
levels  in  the  boundary  layer  will  decrease  the  transport  of  turbulent  energy  toward  the 
wall,  which  in  turn  reduces  the  wall  shear  stress  or  friction  coefficient.  This  view  is 
overly  simplistic  and  cannot  be  accepted  at  measurement  positions  farther  downstream 
where  the  Reynolds  stress  distribution  obtains  an  enhanced  structure  in  the  mid-layer  (see 
e.g.  Figures  3.13  (d)  -  (f)).  Westphal  was  concerned  with  this  phenomenon  and  suggested 
that  a  net  reduction  in  skin  friction,  using  the  flat  plate  type  manipulator,  may  be  minimal 
due  to  the  rapid  recovery  of  the  Reynolds  stresses  through  the  diffusion  process  [7].  The 

Cfm 

fact  that  the  coefficient  of  friction  ratio  is  still  significantly  less  than  one  (7—  ~  0.8)  even 
at  downstream  positions  where  the  Reynolds  stress  profile  is  enhanced  indicates  other 
mechanisms  at  are  work. 

A  physical  description  of  how  the  honeycomb  manipulator  influences  the  flow 
structure  of  the  boundary  layer  is  best  obtained  by  examining  the  equations  for  the  mean 
and  turbulent  kinetic  energies.  Such  equations  are  specified  or  may  be  extracted  from 
analyses  found  in  the  works  of  Hinze  [29]  and  Townsend  [8].  The  framework  for 
understanding  manipulated  flow  is  established  by  applying  these  equations  to  the 
unmanipulated  turbulent  boundary  layer. 

The  mean  kinetic  energy  equation  in  standard  tensor  notation,  neglecting  thermal 
and  gravitational  effects,  is  given  by  Equation  (4.1). 


l3(U;)2  1_  d(uf  d{(p)(Uj)}  _  du'ju'j  _  dnl; 
2  dt  2  J  dxj  dxj  dxj  dxj 


(4.1) 


3  rtu . 

For  steady-state  conditions  (^  =  0)  and  recognizing  the  continuity  relationship  =  0), 


the  expression  is  simplified.  Also,  the  Reynolds  stress  and  viscous  stress  terms  can  be 
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combined  to  give  the  total  stress, 


(4.2) 


For  an  unmanipulated  boundary  layer  the  Reynolds  stress  term  is  the  dominant  stress 

d* 

term  except  near  the  wall  where  ^  is  large.  The  equation  for  mean  kinetic  energy  reduces 
to 


—  3  ui  —  dp  —  dz:j 

u~ — -  =  +  . 

1  dxj  2  1  dxj  dxj 


(4.3) 


If  homogeneity  in  the  cross-flow  direction  is  assumed,  u3  =  w  =  0  and  all  x3  =  z 
derivatives  vanish.  Then  (4.3)  becomes 
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By  consideration  of  velocity  and  length  scales,  albeit  with  different  approaches  in 
neglecting  product  and  derivative  terms  containing  v,  Hinze  and  Townsend  reduced  the 
mean  kinetic  energy  equation  to 


—  3  -  3  Y (w)2  (v)2  —dp  »du  dzu 
dx  dy  2  2  )  dx  dy  dy 


A  A  f  t 

where  x  =  z^  [8,29].  The  left-hand  side  of  (4.5)  represents  the  convection  of  the  mean 

kinetic  energy  within  the  boundary  layer.  The  pressure  term  on  the  right-hand  side 

describes  the  work  done  by  the  imposed  pressure  gradient  across  the  volume.  The  stress 
_ ^ 

terms  together  equal  -u  Thus,  energy  supplied  to  the  mean  flow  via  wind  tunnel 
blower  differential  pressure  is  transported  primarily  by  the  mean  stream  flow  and 
transformed  into  mean  kinetic  energy  and  work  done  by  the  stress  gradient  normal  to  die 
wall.  This  latter  effect  is  characterized  by  an  inward  transfer  of  stresr  energy  (^ (-zu)) 
toward  the  wall  and  conversion  of  mean  flow  energy  to  turbulent  energy  (-x^),  primarily 
in  the  near  wall  sub-layer  [9]. 
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The  effects  of  the  blower  in  supplying  energy  to  the  mean  flow  are  best  understood 
by  integrating  the  pressure  term  in  the  mean  kinetic  energy  equation  over  a  rectangular 
Cartesian  volume  element  within  the  boundary  layer,  differentially  expressed  as  dxdydz. 
The  original  formulation  of  the  pressure  gradient  as  found  in  (4.1)  can  be  used  since 
negligible  terms  or  terms  eliminated  by  the  continuity  relationship  will  not  survive  the 
analysis.  Thus,  the  energy  supplied  to  an  arbitrary  rectangular  volume  within  the 
boundary  layer  can  be  expressed  by 

(m)  (p)dxdy  ,  (4.6) 

V  s2  s, 


J^{(«)(p)}dWydz  =  j(u)(p)dxdy  -  J 


where  S,  and  S2  represent  the  upstream  and  downstream  faces  of  the  rectangular  solid, 
respectively,  and 

j^-{(v)(p)}dxdydz  =  J ( v)(p)dxdz  -  j ( v)(p)dxdz  ,  (4.7) 

v 0y  A?  x, 


where  A,  and  A2  are  the  near  and  far  surfaces,  respectively,  of  the  rectangular  volume 
parallel  to  the  wall.  The  integral  in  (4.7)  vanishes  for  unmanipulated  flow.  In  a  like 
manner  the  other  terms  in  (4.5)  may  be  integrated. 

The  steady-state  turbulent  kinetic  energy  equation,  neglecting  thermal  and 
gravitational  variations,  is  given  by 
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(4.8) 


By  applying  the  same  approximations  and  scales,  as  done  for  the  mean  equation,  to 
eliminate  negligible  terms,  (4.8)  becomes 
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Since  _  _^  +  v^i,  the  mean  and  turbulent  kinetic  energy  equations,  (4.5)  and  (4.9)  can 

3v 

be  added  to  obtain 

M+7£B{<I)!+(7),+7> 

where  J  =  mean  flow  viscous  dissipation. 

Using  an  integral  approach,  the  energy  provided  to  an  arbitrary  volume  by  imposed 
pressure  gradients  in  an  unmanipulated  boundary  layer  is  converted  to  mean  and 
turbulent  viscous  dissipation  (heat)  within  the  volume,  to  a  net  gain  or  loss  in  mean  and 
turbulent  kinetic  energy  within  the  volume,  depending  upon  the  transport  of  total  kinetic 
energy  across  the  boundaries  of  the  volume,  to  a  loss  of  kinetic  and  pressure  turbulent 
energy  (second  term  in  (4.10))  from  the  volume  (net  outward  diffusion),  and  to  a  loss  of 
stress  energy  from  the  volume  (net  inward  transport  toward  the  wall).  The  total  energy 
equation  (4.10),  then,  contains  only  one  significant  mode  of  energy  transfer  to  the  wall, 

3  * _ 

specifically  the  term  ~—(xu).  Between  the  viscous  sub-layer  and  the  outer  region  of  the 
boundary  layer,  mean  stress  energy  is  transported  to  the  wall,  primarily  via  the  Reynolds 
stress.  Analysis  of  the  same  mechanism  in  the  manipulated  boundary  layer  is  therefore 
important. 

In  the  manipulated  cases  two  streamwise  domains  are  of  interest,  that  immediately 
downstream  of  tire  BLM  where  the  flow  structure  retains  identity  with  an  individual 
honeycomb  cell  and  that  where  the  flow  sub-structures  have  fused.  Moller  and  Leehey 
found  that  honeycomb  manipulator  influence  on  wall  shear  and  pressure  is  likewise 
divided  into  streamwise  zones,  one  near  zone  where  the  manipulator  length  scale  is 
dominated  by  honeycomb  cell  diameter  and  a  downstream  zone  where  the  length  scale  is 
dominated  by  manipulator  height  [5].  For  the  experiments  undertaken,  the  boundary  of 
these  two  domains  lies  approximately  ten  centimeters  downstream  of  the  downstream 
face  of  the  BLM.  In  the  near-manipulator  domain  the  use  of  velocity  and,  particularly, 
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spatial  scale  relationships  cannot  be  continued.  The  honeycomb  manipulators  destroy  the 
homogeneity  of  the  mean  flow  on  a  scale  the  size  of  an  individual  cell,  but  not 
necessarily  on  the  whole.  Since  but  one  profile  measurement  was  taken  in  the 
near-manipulator  domain,  insufficient  data  exists  to  evaluate  energy  transport 
mechanisms  by  conducting  a  term-by-term  breakdown  of  Equations  (4.3)  and  (4.8). 
Based  on  the  manipulated  profiles  four  centimeters  downstream  of  the  BLM,  the  stress 

dt,y 

gradients  in  both  x  and  y  directions  (and  perhaps  z)  may  be  significant  in  this  domain. 
Since  Moller  and  Leehey  found  an  increase  in  the  spanwise  ("cross  channel")  coherence 
of  the  wall  pressure  cross-spectral  density  in  a  honeycomb-manipulated  boundary  layer 
[5],  the  manipulator  likely  transforms  the  one -dimensional  energy  transport  mechanism 
of  an  unmanipulated  boundary  layer  into  three-dimensional  ones.  Thus,  the  homogeneity 
of  the  cross-flow  is  likewise  lost.  This  effect  weakens  progressing  downstream  from  the 
manipulator. 

Beyond  the  near-manipulator  domain  the  flow  structure  "remembers"  the 
manipulator’s  effects  but  undergoes  restoration  to  a  form  approximating  the 
unmanipulated  profile.  A  most  dramatic  influence  of  the  manipulator  is  seen  in  the 
strong  suppression  of  the  Reynolds  stress  25  and  65  centimeters  downstream  of  the 
manipulator.  Also,  since  the  scaling  relationships  have  essentially  been  restored  (i.e. 
Equation  (4.5)  applies  again),  the  term  xu  becomes  the  dominant  mode  of  energy 
transport  toward  the  wall.  Comparisons  of  the  profiles  of  the  Reynolds  stress  and  t u 
between  the  manipulated  and  unmanipulated  cases  (Figures  3.13  (a)-(f)  and  4.2  (a)-(f), 
respectively)  provide  insight  into  the  decrease  in  or  Cf  in  the  manipulated  boundary 
layer. 

Progressing  from  the  closest  measurement  position  downstream  of  the  BLM  to  the 
farthest  measurement  position,  the  manipulator’s  influence  on  energy  transport  can  be 
deduced.  Near  the  manipulator  competing  regions  of  intense  Reynolds  and  viscous 
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stresses  are  introduced  into  the  downstream  flow  as  the  honeycomb  cells  destroy  the 
large  scale  turbulence  in  the  flow.  The  cyclic  alternation  of  these  stresses  normal  to  the 
wall  result  in  a  reduced  net  transport  of  energy  into  the  near-wall  layer.  Thus,  wall  shear 
stress  decreases.  As  mentioned  previously,  more  refined  measurements  at  various 
downstream  positions  in  the  near-manipulator  domain  should  indicate  which  terms  of  the 
total  stress  tensor  dominate  the  transport  of  energy.  Diffusive  and  particularly  convective 
transport  of  these  stresses  cause  the  regions  to  merge  such  that  positive  and  negative 
stresses  superpose  to  form  much  lower  Reynolds  stress  levels  compared  to 
unmanipulated  levels,  especially  in  the  mid-layer  region.  A  distinctive  Reynolds  stress 
profile,  with  characteristic  mid-layer  peaking,  emerges  as  turbulence  generated  by  the 
manipulator  increases  in  scale  due  to  convective  and  diffusive  transport.  Near  the  wall 
<  0.2)  far  downstream,  the  Reynolds  stress  and  xu  do  not  recover  to  their 
unmanipulated  levels.  Although  considered  a  "memory"  effect  of  the  manipulator,  the 
physical  explanation  of  the  low  near-wall  stress  levels  is  not  apparent  from  consideration 
of  terms  in  the  kinetic  and  turbulent  energy  equations  nor  from  examination  of  the 
velocity  or  stress  profiles. 

The  physical  explanation  for  the  reduced  wall  shear  stress  at  positions  downstream 
of  the  BLM  where  the  Reynolds  stress  profile  possesses  a  mid-layer  peak  may  be  found 
in  Townsend’s  interpretation  of  the  transport  mechanism  of  turbulent  energy  toward  the 
wall.  The  product  of  T  and  u  is  the  quantity  of  interest.  Figures  4.2  (a)-(f)  provide 
profiles  of  the  specific  stress  energy  transport  quantity  xu  normalized  to  the  product  u*Um 
for  each  measurement  position  at  the  free  stream  speed  of  19  meters  per  second.  The 
near  wall  magnitude  of  xu  in  the  manipulated  cases  is  significantly  less  than  the 
manipulated  case.  Thus,  near  the  wall  the  total  stress  energy  available  for  transport  to  the 
wall  is  reduced  by  the  influence  of  the  manipulator.  Hence,  wall  shear  stress  or  friction 
coefficient  at  these  locations  for  the  manipulated  boundary  layer  is  less  than  the 
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unmanipulated  layer.  The  result  then  of  reduced  wall  shear  stresses  at  the  far 
downstream  locations,  despite  the  mid-layer  peaking  of  the  Reynolds  stress,  is  supported 
by  this  analysis  of  the  data  in  terms  of  a  stress  energy  transport  model. 


Figure  4.2(a)  -j1-  vs  |  (U„=19  m/sec,  xBLM=4  cm) 
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Figure  4.2(d)  vs  j  (U„=19  m/sec,  xBLM=65  cm) 
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Figure  4.2(e)  -p-vs^  (U«=19  m/sec,  xBLM=100  cm) 
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5  Conclusions  and  Recommendations 

5.1  Conclusions 

Numerous  conclusions  follow  from  the  data  obtained  from  X-probe  and  surface 
pressure  measurements  in  honeycomb-manipulated  boundary  layers: 

(1)  Relationships  describing  the  turbulent  energy  decay  and  stieamwise  velocity 
variations  downstream  of  free  stream  honeycomb  manipulators  and  grids  do  not  apply  to 
wall-bounded  turbulence  manipulated  by  honeycomb  structures. 

(2)  Velocity  and  stress  profiles  downstream  of  a  honeycomb  manipulator  are 
characterized  by  two  domains,  a  near-manipulator  domain  where  the  profiles  possess  a 
structure  determined  by  individual  cells  of  the  manipulator  and  a  downstream  domain 
where  the  profiles  lose  the  cell-induced  sub-structures.  For  the  manipulators  used  in  the 
measurements,  the  near-manipulator  domain  ends  approximately  10  centimeters 
downstream  of  the  downstream  face  of  the  BLM. 

(3)  Honeycomb-manipulated  boundary  layers  exhibit  local  drag  reductions  up  to 
40%  compared  to  an  unmanipulated  boundary  layer  at  the  same  free  stream  speed.  The 
greatest  drag  reduction  or  reduction  in  wall  shear  stress  occurs  immediately  downstream 
of  the  manipulator,  and  drag  reduction  persists  150  6*^  downstream  at  a  level  15%  below 
the  unmanipulated  case. 

(4)  Honeycomb  manipulators  introduce  a  structure  into  the  downstream  flow  which 
strongly  suppress  the  transport  of  stress  energy  toward  the  wall,  particularly  via  the 
working  of  Reynolds  stresses.  The  wall-normal  gradient  of  the  transport  term  tm  , 
composed  primarily  of  the  Reynolds  stress  except  near  the  wall,  determines  the  flow  of 
energy  toward  the  wall  beyond  the  near-manipulator  domain.  The  reduction  in  stress 
energy  transported  to  the  wall  appears  to  be  the  principal  cause  of  the  decreased 
magnitude  of  the  wall  shear  stress/friction  coefficient. 
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(5)  Mid- layer  peaking  in  the  Reynolds  stress  levels  of  the  honeycomb-manipulated 
boundary  layers  well  downstream  of  the  manipulator  was  found  to  conform  to  previous 
work  with  flat  plate-type  manipulators.  Although  there  are  increased  levels  of  Reynolds 
stress  in  the  mid-layer,  there  is  less  stress  energy  near  the  wall  compared  to  the 
unmanipulated  boundary  layer,  resulting  in  decreased  wall  shear  stress  downstream  of  the 
manipulator. 

(6)  The  honeycomb  manipulator  produces  sizable  oscillations  in  the  mean  vertical 
velocity  profile  immediately  downstream  of  the  manipulator,  which  decay  significantly 
within  the  near-manipulator  domain.  The  same  behavior  is  exhibited  by  the  mean  stream 
velocity  profile  and  stress  profiles. 

(7)  Manipulated  boundary  layers  do  not  adhere  to  Coles’  law  of  the  wall  although 
some  profiles  were  found  to  closely  approximate  this  law. 

5.2  Recommendations 

Although  all  specified  objectives  were  accomplished,  the  results  obtained  from 
manipulated  boundary  layer  measurements  warrant  similar  investigations  of  expanded 
scale.  To  this  end,  the  following  recommendations  are  offered: 

(1)  Velocity  and  stress  profiles  should  be  measured  at  fine  streamwise  steps  in  the 
near-manipulator  domain  to  explore  in  detail  the  energy  transport  mechanisms 
responsible  for  the  strong  mid-layer  suppression  of  the  Reynolds  stress  in  the  downstream 
domain.  Likewise,  an  increased  number  of  velocity  and  stress  profiles  at  smaller 
streamwise  intervals  should  be  measured  to  determine  similar  mechanisms  accounting  for 
the  gradual  development  of  the  mid-layer  peak  in  Reynolds  stress  levels  farther 
downstream. 
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(2)  The  tops  of  the  manipulators  should  be  tapered  and  smoothed  in  an  attempt  to 
reduce  the  large  shearing  peak  introduced  into  the  velocity  and  stress  profiles  and  to 
determine  any  influence  of  this  phenomenon  upon  energy  transport  via  Reynolds  stress. 

(3)  The  speed  range  of  all  measurements  should  be  expanded  to  obtain  a  more . 
comprehensive  understanding  of  dependencies  of  the  velocity  and  stress  profiles  and  of 
energy  transport  on  Reynolds  numbers. 

(4)  A  consistent  "optimal"  manipulator  height  (hBLM°z x08)  should  be  maintained 

over  the  entire  speed  range  as  xBLM  varies,  and  a  constant  manipulator  height,  above  and 
below  optimal,  should  be  maintained  throughout  the  speed  range  as  xBLM  varies.  In  this 
way  the  variation  of  C(  JCf  um  with  Xb^/S*^  and  hg^/5*^  is  best  determined  (see  Fig. 
4.1).  Then,  based  on  these  findings,  changes  should  be  considered  separately  in 
manipulator  honeycomb  cell  size  and  length. 

(5)  Maintain  X-probe  calibration  using  the  velocity-angle-voltage  technique  to 
ensure  that  mean  vertical  velocity  variations  are  accurately  measured.  Measurement  of 
the  w-component  of  mean  and  fluctuating  velocities  ("cross  flow")  via  triple  hot-wire 
techniques  is  recommended  to  obtain  a  detailed  knowledge  of  all  velocity  and  stress 
influences.  Then,  the  effects  of  all  stress  components  (ztJ  and  u;Uj,  ij  -1,23)  in  energy 
transport  processes  can  be  deduced. 
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Appendix  A  Supplementary  Graphs 

The  body  of  this  document  contained  graphical  representations  of  data  at  the  speed 
of  19  meters  per  second  for  both  manipulated  and  unmanipulated  boundary  layer 
measurements.  For  completeness  and  comparison  the  data  collected  at  15  and  23  meters 
per  second  are  included  in  this  Appendix.  The  order  of  presentation  follows  that  for  19 
meters  per  second  and  includes  both  unmanipulated  and  manipulated  data  as  found  in 
Sections  3.2.1,  3.2.2,  and  3.3.2  only. 

A.l  Fifteen  Meters  Per  Second  Data 


Figure  A. 1.1  Unmanipulated  Mean  Velocity  Profile  (U„=15  m/sec) 
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Figure  A.1.2  —  vSg 


(Unmanipulated,  U„=15  m/sec) 


Figure  A. 13 


(Unmanipulated,  U„=15  m/sec) 
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Figure  A.1.4  Stream  Velocity  Profile  Comparison  (xBLM=4  cm) 


Stream  Velocity  (m/sec) 


Figure  A. 1.5  Stream  Velocity  Profile  Comparison  (xBLN1=12  cm) 
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Stream  Velocity  (m/sec) 
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Figure  A.1.6  Stream  Velocity  Profile  Comparison  (xBLM=25  cm) 


Figure  A.1.7  Stream  Velocity  Profile  Comparison  (xBLM=65  cm) 
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Figure  A. 1.9  Stream  Velocity  Profile  Comparison  (xBLM=150  cm) 
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Figure  A.1.14  (xBLM=12cm) 
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Figure  A.1.24  —  vs  j  (xBLM=150  cm) 


Figure  A. 1.25  —  vs  ^  (xBLM=4  cm) 
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Figure  A.l  JO  -^vs^  (xBLM=150  cm) 


A.2  Twenty-Three  Meters  Per  Second  Data 


0  5  10  15  20  25 

Stream  Velocity  (m/aec) 


Figure  A.2.1  Unmanipulated  Mean  Velocity  Profile  (U_=23  m/sec) 
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Figure  A.2.2  —  vs 


(Unmanipulated,  U„=23  m/sec) 
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Figure  A.23  —  vs  ^  (Unmanipulated,  U„=23  m/sec) 
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Figure  A.2.4  Stream  Velocity  Profile  Comparison  (xBLM=4  cm) 
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Figure  A.2.5  Stream  Velocity  Profile  Comparison  (xBLM=12  cm) 
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Figure  A.2.6  Stream  Velocity  Profile  Comparison  (xBLM=25  cm) 


Figure  A.2.7  Stream  Velocity  Profile  Comparison  (xBLM=65  cm) 


109 


.10 


Figure  A.2.8  Stream  Velocity  Profile  Comparison  (xBLM=100  cm) 


0  5  10  15  20  25 
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Figure  A.2.9  Stream  Velocity  Profile  Comparison  (xBLM=150  cm) 
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Distance  Above  Wall  (m) 


Vertical  Speed  (m/sec) 


Figure  A.2.10  Vertical  Velocity  Profile  (xBLM=4  cm) 


Figure  A.2.11  Vertical  Velocity  Profile  (xBLM=12  cm) 
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Figure  A.2.14  —  vSg  (xBLM=12cm) 
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Figure  A.2.15  —  vs  6  (xblm=25  cm) 
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Figure  A.2.24  -^-vs^  (xBLM=150cm) 


Figure  A.2.25  vs  ^  (xBLM=4  cm) 
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Figure  A.2.28  —7-  vs  ^  (x 
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Figure  A.2.29  (xBLM=100cm) 
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gureA.2.30  (xBLM=150  cm) 
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Appendix  B  Tabular  Data 

Raw  data  collected  for  the  calibration  of  the  X-probe  and  the  unmanipulated  and 
manipulated  turbulent  boundary  layers  is  not  included  in  this  document.  The  data  is 
stored  on  floppy  disks  in  the  Acoustics  and  Vibrations  Laboratory  at  the  Massachusetts 
Institute  of  Technology.  Copies  of  the  data  file  disks  are  also  held  by  the  author  and  his 
thesis  supervisor. 
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Appendix  C  Computer  Programs 


Computer  programs  for  the  collection  and  manipulation  of  X-probe  data,  including 
graphical  representation,  are  found  in  this  Appendix.  All  programs  are  in  the  FORTRAN 
computer  language.  The  programs  were  run  on  the  MASSCOMP  UNIX  system.  The 
data  acquisition  program  for  surface  gauge  differential  pressures  is  not  included.  Mr. 
Yuksel  Gur  of  the  Acoustics  and  Vibration  Laboratory  of  the  Massachusetts  Institute  of 
Technology  courteously  provided  the  program  and  gauge  for  use  in  the  measurements. 
Subroutines  formulated  by  Mr.  Kay  Herbert,  also  of  the  Acoustics  and  Vibrations 
Laboratory,  for  similar  applications  are  included  and  referenced  within  the  programs. 

The  programs  are  presented  in  the  logical  order  of  use.  First  the  X-probe  is 
calibrated  and  the  two-dimensional  polynomial  linear  least  squares  fits  for  u  and  v 
accomplished.  The  coefficients  from  these  fits  are  used  to  transform  X-probe 
instantaneous  voltages  into  instantaneous  velocities  in  the  data  acquisition  program  for 
the  boundary  layer.  Statistics  are  performed  in  the  data  acquisition  program  to  calculate 
mean  velocities  and  rms  values  of  both  fluctuating  components  and  the  covariance  of  the 
fluctuating  components,  the  Reynolds  stress.  The  remaining  programs  use  the  output 
data  of  boundary  layer  measurements  and  the  differential  pressures  from  the  surface 
differential  pressure  gauge  to  calculate  boundary  layer  parameters  and  present  the  data  in 
a  format  consistent  with  historical  representations. 

There  are  two  separate  programming  paths  to  follow  depending  upon  whether  the 
boundary  layer  was  manipulated  or  unmanipulated.  If  unmanipulated,  the  Coles  fitting 
algorithm  uses  mean  velocity  profile  data  to  determine  both  friction  velocity,  y-offset, 
and  other  parameters  such  as  boundary  layer  displacement  and  momentum  thicknesses. 
Then,  differential  pressures  are  used  in  conjunction  with  corresponding  friction  velocities 
to  determine  a  calibration  curve  of  wall  shear  stress  as  a  function  of  differential  pressure. 
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If  manipulated  data  is  obtained,  the  differential  pressure  gauge  readings  taken  with  each 

manipulated  profile  are  averaged  and  used  to  determine  friction  velocity.  The  friction 

velocity  is  used  in  an  interactive  Coles  fitting  algorithm  (visual  slope  matching  of 

historical  and  measurement  data)  to  determine  y-offset  of  the  manipulated  data  and 

subsequently  other  boundary  layer  parameters.  In  the  plotting  programs  graphical  profile 

representation  corresponds  to  the  format  used  by  Klebanoff  [9]. 

PROGRAM  XVANGCAL:  Data  Acquisition  of  Calibration  Voltages  (X-probe) 

and  Pitot  Pressure  (Converted  to  Velocity)  at  Each 
Calibration  Angle 

Input  Output 

Pitot  gain  velocity-angle-voltage  data  files 

temperature  (q-a-volt(l)-volt(2)) 

angle  velocity-velocity-voltage  data  files 

(u-v-volt(  1  )-volt(2)) 


c* **************************** ***************************** 
c  This  program  collects  pressure  and  voltage  data 
c  for  an  X-probe  hot-wire  pair  placed  in  a  stream, 
c  One  pressure  (chan2)  and  two  voltage  channels 
c  are  used,  and  the  pressure  data  is  convened  to 
c  velocities.  Data  formatted  into  four  columns 
c  is  placed  into  four  output  files,  two  in  q-a  format 
c  and  two  in  u-v  format.  Two  output  files  are  for 
c  quick  analysis  to  determine  calibration  coefficients, 
c  These  files  eliminate  low  velocity  data  points  where 
c  mutual  heat  transfer  between  sensors  distorts  voltage 
c  data.  To  include  these  points,  the  two  other  files 
c  can  be  modified  by  eliminating  non-entries, 
c  All  data  acquisition  terms  are  identified  in  the 
c  MASSCOMP  Data  Acquisition  Manual  with  example 

c  programs  provided. 

c  Note:  When  preparing  for  data  collection,  ensure 
c  that  the  X-probe  is  set  to  the  correct  angle  and 
c  small  increments  of  increasing  velocity  are  made 
c  for  the  first  few  data  points 
c 

c  Frank  Camelio  September  1989 

c** ************************* ********** ********************** 
program  xvangcal 
include  ’/usr/include/mr.f 
integer  RDWR,  NEARFRQ,  LOW,  NOTUSED 
integer  NFRAMES,  NCHAN 
parameter  (NFRAMES  =  80000) 
c  Number  selected  for  NFRAMES  can  be  increased  if  desired  but 
c  experience  has  shown  that  80000  sample  points  are  sufficient 
c  for  mean  velocity  data. 
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parameter  (NCHAN  =  3) 
parameter  (RDWR  =  0) 
parameter  (NEARFRQ  =  0) 
parameter  (FREQ  =  8000.0) 
parameter  (BRSFRQ  =  1000000.0) 
parameter  (LOW  =  0) 
parameter  (NOTUSED  =  0) 
character  frameclkdev*(*) 
character  burstclkdev*(*) 
character  addev*(*) 

parameter  (frameclkdev  =  ’/dev/dacpO/clkO’) 
parameter  (burstclkdev  =  ’/dev/dacpO/cIkl  ’) 
parameter  (addev  =  ’/dev/dacpO/adfO’) 
integer  adpn,  clkpnl,  clkpn2,  fchan,  incr,  gain 
integer  byteslocked 
integer  status(2) 

integer*2  rawdata(NCHAN*NFRAMES) 
real  rfreq.rbrst 

real*4  ang,q(25),vel(25),ve2(25),pe(25),u(25),v(25) 
c**** **************** ******************************************* 
c  ang=angle  manually  set  for  calibration  velocity  range 
c  q=tunnel  stream  velocity 
c  vel=voltage  measured  from  X-probe  channel  1 
c  ve2=voltage  measured  from  X-probe  channel  2 
c  pe=voltage  measure  from  Pitot  (Validyne) 
c  (converted  to  velocity  for  the  output  files) 
c  u=q*cos(ang),  velocity  component  parallel  to  probe  axis 
c  v=q*sin(ang),  velocity  component  normal  to  probe  axis 
c  an  parallel  to  probe  support  stem 

c  Configured  for  a  maximum  of  25  velocity  data  points 
c  for  each  angle  set  during  the  calibratioa 
c  For  each  angle  step  through  25  or  less  (recommend 
c  12  minimum)  RPM's  (velocities)  repeating  the  same 
c  RPM’s  at  the  other  angles. 

£***  ft  ♦  *  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ft*  ♦  ♦  ^  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ^  ^  It  ft  ♦  ft  ♦  ♦  ^  ft  ♦  ♦  ♦  ♦  ♦  ft  ♦  ^  ^  ♦  ft  ♦  ft**  ♦  ft  ft  ft  ♦ 

character*  10  date 

common  rawdata 

data  adpn/-l  / 

data  clkpnl  /- 1/ 

data  clkpn2  /- 1/ 

call  mrlock(0,0,byteslocked) 

write(6,*)  byteslocked,  ’  bytes  locked  in  memory’ 

write(6,*) 

c  frnt  =  format  files  containing  all  determined  parameters  &  text 
open(8,FILE=’xqacal.dat\STATUS=’fresh’) 
open(85  JTLE= ’xuvcal  .dat  ’  ,STATUS=’ fresh ’) 
open(9,FILE= ’  xqacal.fmt  ’  ,ST  ATU  S= ’  fresh’ ) 
open(95,FTLE=’xuvcal.fint\STATUS=’fresh’) 
write(9,l) 

1  format(3x,’q-vel’,4x, ’angle’, 4x,’el’,5x,’e2’) 
write(9,2) 

2  format(3x, ’m/sec’, 5x, ’rad’, 4x,’volts’,3x, ’volts’) 
write(9,*) 

write(95,3) 

3  format(3x,’u-ver,4x,’v-vel’,4x,’el’,5x,’e2’) 
write(95,4) 

4  format(3x, ’m/sec’, 4x, ’m/sec’, 3x,’volts’,3x,’volts’) 
write(95,*) 
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write(6.*)  ’chan2=pitot  chanO=hotwire_l  ’ 
write(6,*)  ’chanl=hotwire_2’ 
fchan  =0 
write(6,*) 

write(6 ,*)  ’Enter  date  as  a  10  character  string  (no  spaces).’ 
read(5,*)  date 
write(6  ,*) 

write(6,*)  ’Input  gain  setting  for  pitot  ’ 
read(5,*)  fgain 
write(6  ,*) 

write(6,*)  ’Enter  temperature  in  deg  Fahrenheit.’ 

read(5,*)  tempf 

write(6,*) 

temp=<tempf-32.)/l  .8 
iho=l  .293/(1  .+.00367*temp) 

1*0 

vavl=0 
vav2=0 
5  k=l 
1=1+1 
write(6,*) 

write^,*)  ’Secure  fan,  reset  RPM  to  0,  set  traverse  angle.’ 

write(6,*)  ’Enter  angle  setting  of  probe  relative  to  the’ 

write(6,*)  ’stream  flow  ( +  or  -  in  degrees).’ 

read(5,*)  ang 

ang=ang*3 . 14159/180 

write(6,*) 

write(6,*)  ’start  at  velocity=0  ’ 

10  write(6 ,*)  ’take  data,  re-set  angle,  or  quit?’ 

write(6,*)  ’(  >  0  =  data,  0=re-set,  <  0  =  quit)  ’ 
c* **************************** ****** ****** ************ ********** 
c  Each  time  re-set  is  chosen,  values  for  a  given  calibration 
c  run  are  written  to  file  and  the  angle  must  be  changed, 
c  After  calibration  on  final  angle  setting,  enter  -1  to 
c  write  to  file  and  terminate  program. 

c*************************************************************** 
read(5,*)  connum 
if  (connum.le.0)  goto  25 
write(6,*) 

call  mropen(adpn,addev,RDWR) 

call  mropen(clkpnl,frameclkdev,RDWR) 

call  mropen(clkpn2,burstclkdev,RDWR) 

call  mrclk2(clkpnl,  clkpn2,  NOTUSED,  NEARFRQ,  FREQ, 

&  rfreq,  NEARFRQ,  BRSFRQ,  rbrst,  NCHAN,  LOW) 

write(6,*)  ’freq  for  clk6  =  ’rfreq 

write(6,*)  ’ffeq  for  clk7  =  ’rbrst 

incr=l 

gain=0 

call  mradxin(adpn,clkpnl  ,clkpn2, fchan, NCHANrncr, gain, 

&  NFRAMES.rawdata) 
write(6,*)  ’transfer  starting ....’ 
call  mrevwt(adpn,  status, 80000) 
write(6,*)  ’transfer  complete’ 
call  mrclosall 

c  The  following  terms  are  pressure,  velocity  rawdata  voltage  sums. 
pvol=0. 
woll=0. 
wol2=0. 

do  20  i=l,NFRAMES*NCHAN-2,3 
vvol  1  =vvol  1  +float(rawdata(i )) 
wol2=vvol2+float(rawdata(i+l )) 

20  pvol=pvol+float(rawdata(i+2)) 
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pe(kl=pvoI*  .00244 1 4 1 /floaUNFR  AMES)/fgain 
ve  1  (k)=vvol  1  * . 00244 141 /float(NFRAMES ) 
ve2(k)=wol2*.00244141/float(NFRAMES) 
c  The  following  line  grouping  is  entered  to  detect  problems  in  the 
c  course  of  calibration.  Calibration  is  conducted  from  low  to 
c  high  values  of  pressure  and  velocity.  If  a  voltage  on  any 
c  channel  is  less  than  a  preceding  voltage  measurement,  the 
c  appropriate  warning  is  sent  to  screen.  For  example,  hot-wire 
c  resistance  changes  can  indicate  a  potential  failure,  or  more 
c  commonly  at  low  velocities  near  zero,  the  flow  may  not  have 
c  stabilized  and  pitot  voltage  is  erratic. 
write(6,*) 

write(6,*)  pe(k),vel(k),ve2(k) 
write(6,*) 
if(k.gt.l)  then 
if(pe(k).lt.pe(k-l»  then 

write(6,*)  ’CAUTION:  RE-SETTING  -  P  VALUE  <  PRECEDING’ 
1=1-1 
goto  5 
endif 

if(ve  1  (k).lt.ve  l(k- 1 ))  then 

write(6 ,*)  ’CAUTION:  RE-SETTING  -  VI  VALUE  <  PRECEDING’ 
endif 

if(ve2(k).lt.ve2(k-l))  then 

write(6,*)  ’CAUTION:  RE-SETTING  -  V2  VALUE  <  PRECEDING’ 
1=1-1 
goto  5 
endif 
endif 
k=k+l 
goto  10 

25  pe0=pe(l) 
k=k-l 
do  30  j=ljc 

ccc  Set  zero  offset  for  pressure  at  each  calibration  angle. 
q(j)=sqrt(  19.6 1 4*pe(j)/rho) 

u0)=tiCj)*cos(ans) 

v(j)=q(j)*sin(ang) 
if(j.eq.l.or.q(j).lt.l.5)  goto  28 
write(8,*)  q(j),ang,vel(j),ve2(j) 
write(85,*)  u(j),v(j),vel(j),ve2(j) 
write(6,*)  q(j),ang,vel(j),ve2(j) 

28  write(9,*)  q(j),ang,vel(j),ve2(j) 

30  write(95,*)  u(j),v(j),vel(j),ve2(j) 

vavl=vavl+vel(l) 
vav2=va  v2+ve2(  1 ) 
if  (connum.eq.0)  goto  5 
if  (connum.lt.0)  goto  50 
50  vavl=vavl/l 
vav2=vav2/l 

write(6,*)  0.0,0.0,vavl,vav2 
write(8  *)  0.0,0.0,vavl,vav2 
write(85,*)  0.0,0.0,vavl,vav2 
write(9,*) 
write(9,*) 

c  The  average  zero  velocity  voltage  is  calculated  for  info  only. 
write(9,"r)  ’eavl(0)=’,vavl,’  eav2(0)=’,vav2 
write(9,*)  ’Ch0=pitot  Chl=hotwirel  Ch2=hotwire2’ 
write(9,*)  ’Pitot  gain  =  ’.fgain 
write(9,*)  ’Temp  (F)  =  ’.tempf 
write(9,*)  ’Date:  ’.date 
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write(9,*)  'Filename:  xqacal.fmt' 

write(95,*) 

write(95,*) 

write(95,*)  ’eavl(0)=',vavl,’  eav2(0)=’,vav2 

write(95,*)  ’ChO=pitot  Chl=hotwirel  Ch2=hotwire2’ 

write(95,*)  'Pitot  gain  =  \fgain 

write(95,*)  'Temp  (F)  =  \tempf 

write(95,*)  'Date:  ’.date 

write(95,*)  'Filename:  xuvcal.fmt’ 

close(8) 

close(85) 

close(9) 

close(95) 

stop 

end 

c**************************************************************** 

PROGRAM  XCOEFF:  Linear  Least  Squares  Fit  of  Velocity  -  Angle 

Calibration  Data  to  Determine  Coefficients  of  a 
Two-Dimensional  Polynomial  Representation  u  and  v 
or  Q  and . 


Input 

data  file  name 
data  file  type 
number  of  coefficients 


Output 

coefficients  files 
coefficients  and  error  files 


q  *  *  *  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  *  %  *  *  **  *  He  4c  *  *  *  *  4c  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

c  This  program  reads  velocity-angle-voltage  data  or  stream 
c  and  vertical  velocity  data  and  fits  separately  a  double 
c  polynomial  surface,  one  in  velocity  and  one  in  angle  OR 
c  one  in  u-velocity  and  one  in  v-velocity,  to  voltage  data 
c  on  both  channels  of  an  X-probe.  The  data  is  obtained  from 
c  a  calibration  which  pivots  the  probe  at  different  positive 
c  and  negative  angles  with  respect  to  the  free  stream, 
c  Polynomial  order  is  either  3,  giving  16  coefficients,  or 
c  4,  giving  25  coefficients, 
c  The  2-d  polynomial  format  is: 
c  N  N 

c  f(el,e2)=SUM  SUM  a(i  j)’*tel**<J-l)*e2**(i-l)  [N=4  or  5] 
c  i=l  j=l 

c 

c  where  el,e2  are  channel  voltages  and  a(ij)  are  the  fitted 
c  coefficients.  For  simplicity  in  analysis  coefficients  are 
c  transferred  to  vector  format  vice  an  array. 

£*  ft  £  ♦  *  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  *  ♦  ★  ♦  ♦  ♦  ♦  ♦  ♦  •  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦ 

program  xcoeff 

implicit  real*8  (a-b,o-z) 

real*8  qu(200,l),av(200,l) 

real  *8  quest(200),qudel(200),avest(200),avdel(200) 

real*8  etqu(25,l),etav(25,l),qufest(200),avfest(200) 

real*8ee(200,25),eet(25,200),ete(25,25) 

real*8  qumaxerr,quaverT,qusumerT,quvarsum,quvar 

real*8  avmaxerr,avaverr,avsumerr,awarsum,awar 

real*8  qurel(200),avrel(200) 

c*************** *********************** ***************************** 

c  qu=q  or  u  velocity  obtained  from  calibration  file 
c  av=angie  or  v  velocity  obtained  from  calibration  file 
c  ...err=error  terms  in  comparing  fitted  to  measured  values 
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c  ...sum=term  used  in  determining  avg  error 
c  ...rel=relative  error  term 
c  ...del=difference  term 

c  ...fest=term  estimated  after  functional  calculation 
c  ...est=estimated  term  from  least  squares  fit 
c  ee=two  channel  voltage  product  in  polynomial 
c  eet,ete,etqu,etav=arrays  of  voltages  and  velocities  used 
c  in  matrix  multiplication  (t=transpose) 

c  qul6,qu25=resulting  functions  for  calculating  q  or  u 
c  avl6,av25=resulting  functions  for  calculating  angle  or  velocity 
c  ao,bo=coefficients  for  qu  and  av  functions  respectively 
c******************************************************************** 
character*15  name, fileout  l,fileout2 
common  /any/  xx(25,27)ix(25) 
common  /xco/  el(200),e2(200) 
common  /funcalc/  ao(25),bo(25) 
real  qul6,avl6,qu25,av25 
externa]  qu  1 6,av  1 6,qu25,av25 
50  write(6,*)  ’Enter  name  of  calibration  data  file.’ 
read(5,*)  name 
write(6,*) 

c*** ************************************************* ********* 
c  READ  VALUES  FROM  Q-A  OR  U-V  CALIBRATION  .DAT  FILE 
c*********** *************************************** *********** 
open(5 1  ,file=name,status=  ’old  ’ ) 
m=l 

100  read(51,*,end=150)  qu(m,l),av(m,l),el(m),e2(m) 

m=m+l 
goto  100 
150  close(51) 
m=m-l 

write(6,*)  ’Enter  1  if  data  is  q-a,  2  if  u-v.’ 

read(5,*)  num2 

write(6,*) 

200  write(6,*)  ’Enter  #  of  coefficients  to  fit  (INTEGER).’ 
read(5  *)numl 
write(6,*) 

if  (numl.eq.l6.and.num2.eq.l)  then 
fileoutl=’xl6coqa.dat’ 
fileout2=’x  1 6coqa.  fmt  ’ 
endif 

if  (numl.eq.l6.and.num2.eq.2)  then 
fileoutl=’xl6couv.dat’ 
fileout2=’x  1 6couv.fmt  ’ 
endif 

if  (numl.eq.25.and.num2.eq.l)  then 
fileout  1 = ’x25coqa.dat  ’ 
fileout2=’x25coqa.  frnt  ’ 
endif 

if  (numl.eq.25.and.num2.eq.2)  then 
fileout  1 = ’x25couv  ,dat  ’ 
fileout2=’x25couv.frnt* 
endif 

open(7,file=fileoutl  ,status=’ fresh  ’) 
open(8,file=fileout2,status=’fresh’) 
if(num2.eq.l)  write(8,*)  ’  q  coefF,’  acoeCT 

if(num2.eq.2)  write(8,*)  ’  u  coeff’,’  v  coeff 

300  if  (numl.eq.16)  call  coeffl6(m,ee) 
if  (numl.eq.25)  call  coeff25(m,ee) 
do  400  i=l,m 
do  400  j=l.numl 
400  eet(j4)=ee(i,j) 
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call  matmlt(eet,ee,ete.num  1  ,m,num  1 ) 
call  matmlt(eet,qu,etqu,num  1  ,m,  1 ) 
call  matmlt(eet,av,etav,numl,m,l) 
do  500  i=l,numl 
do  500  j=l,numl 
500  xx(ij)=ete(i,j) 
do  600  i=I,numl 
xx(i,num  1 +1  )=etqu(i,  1 ) 

600  xx(i,num  1  +2)=etav(i,  1 ) 
call  algrtnuml,2,l,e) 
do  700  i=l,numl 
ao(i)=xx(i  jium  1 +1 ) 
bo(i)=xx(i,num  1  +2) 
write(6,*)  ao(i),bo(i) 
write(7  ,*)  ao(i),bo(i) 

700  write(8,*)  ao(i),bo(i) 
write(8,*) 
write(8,*) 
close(7) 

£****************************♦******  **********^  ********** 
c  ERROR  ANALYSIS 

c********************************************************* 
write(8,*) 
if(num2.eq.l)  then 

write(8,*)  ’  q-act  q-fit  \’q-funct  ’ 
delta-qY  q  rel  err’ 
endif 

if(num2.eq.2)  then 

write(8,*)  ’  u-act  u-fit  Yu-funct  ’ 

&,’  delta-u’,’  u  rel  err’ 
endif 

do  800  i=l,m 

if(num  1  .eq.  1 6)  qufest(i)=qu  1 6(e  1  (i),e2(i)) 
if(num  1  .eq.25)  qufest(i)=qu25(e  1  (i),e2(i)) 
quest(i)=0 
do  750j=l,numl 

750  quest(i)=quest(i)+ao(j)*ee(i,j) 
qudel(i)=qu(i,l  )-quest(i) 
if(qu(i,l).eq.0.)  then 
qurel(i)=0. 
else 

qurel(i)=100*qudel(i)/qu(i,l ) 
endif 

if(qurel(i).eq.O.)  then 

write(8,760)  qu(i,l),quest(i),qufest(i),qudel(i) 

760  format(4f9.5,3x) 

else 

write(8,770)  qu(i,l),quest(i),qufest(i),qudel(i),qurel(i) 

770  format(5f9.5,3x) 

endif 

if(i.eq.l)  qumaxerT=abs(qudel(i)) 
qumaxerr=max(abs(qudel(i)),qumaxerr) 

800  continue 
write(8,*) 
qusumerr=0 
do  900  i=l,m 

900  qusumerr=qusumerr+abs(qudel(i)) 
quaverr=qusumerr/dble(m) 
quvarsum=0 
do  1000  i=l,m 

1000  quvarsum=quvarsum+(abs(qudel(i))-quaverr)**2 
quvar=quvarsum/dble(m) 
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write(8,*) 
if(num2.eq.l)  then 

write(8,*)  ’Max  q  vel  euor  magnitude  (m/sec)  =  \qumaxerr 
write(8,*)  ’Avg  q  vel  error  magnitude  (m/sec)  =  \quaverr 
write(8,*)  ’Var  q  vel  error  magnitude  (m/sec)  =  \quvar 
endif 

if(num2.eq.2)  then 

write(8,*)  ’Max  u  vel  error  magnitude  (m/sec)  =  \qumaxerr 
write(8,*)  ’Avg  u  vel  error  magnitude  (m/sec)  =  \quaverr 
write(8,*)  ’Var  u  vel  error  magnitude  (m/sec)  =  ’.quvar 
endif 
write(8,*) 
write(8,*) 
if(num2.eq.l)then 

write(8,*)  ’  a-act  ’,’  a-fit  ’,’  a-funct’,’  delta-a’ 
a  rel  err’ 
endif 

if(num2.eq.2)  then 

write(8,*)  ’  v-act  ’,’  v-fit  ’,’  v-fiinct’,’  delta-v’ 
v  rel  err’ 
endif 

do  1100  i=l,m 

if(num  1  .eq.  16)  avfest(i)=avl 6(el(i),e2(i» 
if(numl.eq.25)  avfest(i)=av25(el(i),e2(i)) 
avest(i)=0 
do  1050j=l,numl 

1050  avest(i)=avest(i)+bo(j)*ee(i,j) 
avdel(i)=av(i,  1  )-avest(i) 
if(av(i,l).eq.O.)  then 
avrel(i)=0. 
else 

avrel(i)=100*avdel(i)/av(i,l) 

endif 

if(avrel(i).eq.0.)  then 

write(8,760)  av(i,l),avest(i),avfest(i),avdel(i) 
else 

write(8,770)  av(i,l  ),avest(i),avfest(i),avdel(i),avrel(i) 
endif 

if(i.eq.l)  avmaxerr=abs(avdel(i)) 
avmaxerp=max(abs(avdel(i)),avmaxerr) 

1100  continue 
write(8,*) 
avsumerr==0 
do  1200  i=l,m 

1200  avsumerr=avsumerr+abs(avdel(i)) 
avaverr=avsumerr/dble(m ) 
awarsum=0 
do  1300i=l,m 

1 300  avvarsum=awarsum+(abs(avdel(j))-avaverT)**2 
awar=awarsum/dble(m) 
write(8,*) 
if(num2.eq.l)  then 

write(8,*)  ’Max  ang  error  magnitude  (rad)  =  ’.avmaxerr 
write(8,%)  ’Avg  ang  error  magnitude  (rad)  =  \avaverr 
write(8,*)  ’Var  ang  error  magnitude  (rad)  =  \awar 
endif 

if(num2.eq.2)  then 

write(8,*)  ’Max  v  vel  error  magnitude  (m/sec)  =  ’.avmaxerr 
write(8,*)  ’Avg  v  vel  error  magnitude  (m/sec)  =  \avaverr 
write(8,*)  ’Var  v  vel  error  magnitude  (m/sec)  =  \awar 
endif 
write(8,*) 
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write(6.*)  'Enter  3  to  change  #  coefficients.’ 
write(6.*)  'Enter  4  to  change  q-a  to  u-v  or  vice-versa.  ’ 
write(6,*)  ’Else  enter  any  other  integer.’ 
read(5,*)  num3 
write(6  ,*) 

if(num3.eq.3)  goto  200 

if(num3.eq.4)  goto  50 

close(8) 

stop 

end 

g****************************  ************************************* 

c  SUBROUTINES  -  PRODUCTS  OF  VOLTAGES  (16  AND  25  COEFF’s) 
£***** *************************************************** ********* 
subroutine  coeffl6(m,ee) 
implicit  real*8  (a-h,o-z) 
dimension  ee(200,25) 
common  /xco/  el(200),e2(200) 
do  10  i=ljn 
ee(i,l)=l 
ee(i^2)=e  1  (i) 
ee(i,3)=el(i)*el(i) 
ee(i,4)=e  1  (i)*ee(i,3) 
ee(i,5)=e2(i) 
ee(i,6)=e2(i)*el(i) 
ee(i,7)=e2(i)*ee(i,3) 
ee(i,8)=e2(i)*ee(i,4) 
ee(i,9)=e2(i)*e2(i) 
ee(i,10)=e2(i)*ee(i,6) 
ee(i,l  l)=e2(i)*ee(i,7) 
ee(i,  12)=e2(i)*ee(i,8) 
ee(i,13)=e2(i)*ee(i,9) 
ee(i,14)=e2(i)*ee(i,10) 
ee(i,  1 5  )=e2(i)*ee(i,  11) 
ee(i,  1 6)=e2(i)*ee(i,  12) 

10  continue 
return 
end 

£*****************************  ************************************** 
subroutine  coeff25(m,ee) 
implicit  real*  8  (a-h,o-z) 
dimension  ee(200,25) 
common  /xco/  el(200),e2(200) 
do  10i=l,m 
ee(i,l)=l 
ee(i,2)=el(i) 
ee(i,3)=el(i)*el(i) 
ee(i,4)=el(i)*ee(i,3) 
ee(i,5)=el(i)*ee(i,4) 
ee(i,6)=e2(i) 
ee(i,7)=e2(i)*el(i) 
ee(i,8)=e2(i)*ee(i,3) 
ee(i,9)=e2(i)*ee(i,4) 
ee(i,10)=e2(i)*ee(i,5) 
ee(i,ll)=e2(i)*ee(i,6) 
ee(i,12)=e2(i)*ee(i,7) 
ee(i,13)=e2(i)*ee(i,8) 
ee(i,l  4)=e2(i)*ee(i,9) 
ee(i,15)=e2(i)*ee(i,  10) 
ee(i,16)=e2(i)*ee(i,  1 1 ) 
ee(i,  1 7)=e2(i)*ee(i,  12) 
ee(i,18)=c2(i)*ee(i,13) 
ee(i,19)=e2(i)*ee(i,14) 


ee(i.20)=e2(i)*ee(i,15) 
ee(i  J21  )=e2(i)*ee(i,16) 
ee(i.22)=e2(i)*ee(i,17) 
ee(i,23  )=e2(i)*ee(i,  18) 
ee(i,24)=e2(i)*ee(i,  19) 
ee(i,25)=e2(i)*ee(i.20) 

10  continue 
return 
end 

g***************************** ************************************** 

c  SUBROUTINES  -  MATRIX  MULTIPLICATION,  SIMUL  EQN  SOLUTION 
c  (see  Carnahan,  Applied  Numerical  Methods,  Chapter  4) 

c*********** ****************************** ************************** 

subroutine  matmlt(aa,u,t,mm,nn,pp) 
real*8  aa(25,200),u(200 ,25)4(25.25) 
integer  pp 
do  1  i=l,mm 
do  1  j=l,pp 

1  t(i,j)=0 

do  2  i=l,mm 
do  2  j=l,pp 
do  2  k=l,nn 

2  t(i,j)=aa(ijk)*u(k,j)+t(i,j) 
return 

end 

c************* *************************************** *************** 

c  algr  solve  a  real  algebraic  system 
c 

c  id=l :  reduce,  solve  for  nr  right  hand  sides 
c  id=0,  nr>0:  solve  for  nr  right  hand  sides 
c  id=0,  m=0:  wielandt  correction  step 
c  id=- 1 :  reduce,  first  wielandt  step 
c 

c  Kay  Herbert  1987 

c***************************** ***************** ********************* 
subroutine  algr(n,nr,id,e) 
implicit  real*8  (a-h.o-z) 
common  /arry/  a(25,27),ix(25) 

m=n+nr 

ml=n-l 

nl=n+l 

if(id.eq.O)  goto  10 
do  1  k=l,ml 
xx=0. 
kx=k 

do  2  kk~k,n 

x=dabs(a(kkjc)) 

if  (x.le.xx)  go  to  2 

xx=x 

kx=kk 

2  continue 
ix(k)=kx 

if  (kx.eq.k)  go  to  3 
do  4  jj=k,m 
z=a(k,jj) 
a(k,jj)=a(kx,jj) 

4  a(kx,jj>=z 

3  continue 
kl=k+l 

do  1  kk=kl,n 

if  (dabs(a(kjc)).le.0.)  go  to  1 
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a(kk  M  )=a(kk,k)/a(  k.k ) 
do  5jj=kl,m 

5  a(kk,jj)=a(kk,jj)-a(kktk)*a(k,ij) 

I  continue 

if  (id.gt.O)  go  to  12 

a(njil)=l. 
do  6  kk=l,ml 
k=n-kk 

a(kjil)=-a(k,n) 
if  (k.eq.ml)  go  to  6 
kl=k+l 
do  7  l=kl  jnl 

7  a(k,nl)=a(k4il)-a(k,l)*a(l,nl) 

6  a(k,nl)=a(k,nl)/a(k40 
return 

10  if  (nr.le.O)  m=nl 
do  11  l=nl.m 
do  11  k=l,ml 
kx=ix(k) 
z=a(k,l) 
a(k,l)=a(kx,l) 
a(kx,l)=z 
kl=k+l 
do  11  kk=kl,n 

II  a(kk,l)=a(kk,l)-a(kkk)*a(k,l) 
if  (nr.gt.O)  go  to  12 

e=a(n,n)/a(n,m) 
do  13  k=l,n 

13  a(k,m)=e*a(k,m) 

12  if  (m.le.n)  return 
do  14  kk=l,n 
k=nl-kk 

if  (k.eq.n)  go  to  15 
kl=k+l 
do  20  l=nl,m 
do  20  U=kl,n 

20  a(kj)=a(k,l)-a(k41)*a(114) 

15  dol41=nl,m 

14  a(k4)=a(k,I)/a(k4c) 

return 

end 

£*****«************************************•********************** 
c  FUNCTIONS  TO  CONVERT  VOLTAGES  TO  VELOCITY  AND  ANGLE 
£*  £  £  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ^  ♦  ♦♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  *  ^  ♦*  ♦  ♦  ♦  ♦  ♦♦  ♦  ♦  ♦  ♦  ♦  ^  ♦  *  ♦  ♦  ♦  ♦  ^  ^  £  £  *  *  £  £ 

function  qul6(vl,v2) 

implicit  real*8  (e,f,v) 

common  /funcalc/  e(25),f(25) 

qu  1 6=e(  1  )+v  1  *(e(2)+v  1  *(e(3)+v  1  *e(4)))+ 

&v2*(e(5)+vl*(e(6)+vl*(e(7)+vl*e(8)»+ 

&v2*(e(9)+v  1  *(e(  1 0)+v  1  *(e(  1 1  )+v  1  *e(  1 2)))+ 

&v2*(e(  1 3  )+v  1  *(e(  1 4)+v  1  *(e(  1 5)+v  1  *e(  1 6) ))))) 

return 

end 

c********* ****************************** ************************** 

function  avl6(vl,v2) 
implicit  real*8  (e,f,v) 
common  /funcalc/  e(25),f(25) 
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av  1 6=f(  1  )+v  1  *(f(2  )+v  1  *(f(3  )+v  1  *f(4)))+ 

&v2*(f(5)+vl*(f(6)+vl*(f(7)+vl*f(8)))+ 

&v2*(f(9)+vl*(f(10)+vl*(f(l  l)+vl*f(12)))+ 

&v2*(f(  13)+v  1  *(f(  14)+v  1  *(f(  1 5  )+v  1  *f(  1 6)))))) 

return 

end 

£****************************************************************** 
function  qu25(vl,v2) 
implicit  real  *8  (g,h,v) 
common  /fiincalc/  g(25),h(25) 
qu25=g(  1  )+v  1  *(g(2)+v  1  *(g(3)+v  1  *(g(4)+v  1  *g(5))»+ 
&v2*(g(6)+vl*(g(7)+vl*(g(8)+vl*(g(9)+vl*g(I0)»)+ 

&v2*(g(  1 1  )+vl  *(g(  12)+v  1  *(g(  1 3)+v  1  *(g(  14)+v  1  *g(  15))))+ 

&v2*(g(  1 6)+vl  *(g(  1 7)+v  1  *(g(  1 8)+ v  1  *(g(  1 9)+ v  1  *g(20))))+ 

&v2*(g(21)+vl*(g(22)+vl*(g(23)+vl*(g(24>+vl*g(25)))))))) 

return 

end 

£**♦***♦**♦******♦********%****♦♦*♦******+****♦*******%*********** 

function  av25(vl,v2) 
implicit  real*8  (g,h.v) 
common  /funcalc/  g(25),h(25) 
av25=h(  1  >+v  1  *(h(2)+v  1  *(h(3)+v  1  *(h(4)+v  1  *h(5))))+ 
&v2*(b(6)+vl*(h(7)+vl*(h(8)+vl*(h(9)+vl*h(10))))+ 

&v2*(h(  1 1  )+v  1  *(h(  1 2)+v  1  *(h(  1 3)+v  1  *(h(  1 4)+v  1  *b(  15))))+ 

&v2*(h(  1 6)+v  1  *(h(  1 7)+vl*(h(  1 8)+vl  *(b(19)+vl*b(20))))+ 

&v2*(h(21)+vl*(h(22)+vl*(h(23)+vl*(h(24)+vl*h(25)))))))) 

return 

end 

c  ***  **  ***  *  *****  *  **  *  Hr  *  *  *  *  *  *  *  *  *  *  *  *  *  ***  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  **  *  ****  *  *  *  *  * 

PROGRAM  XTBLDATACO:  Boundary  Layer  Data  Acquisition  Program  Utilizing 

an  X-probe  and  Two-Dimensional  Polynomial 
Calibration  Functions  —  Outputs  Mean,  Fluctuating, 
and  Reynolds  Stress  Velocity  Values  at  Each 
Measurement  Position  in  the  Boundary  Layer 


Output 


Input 

coefficients  file  name 
temperature 
stream  location 
traverse  control  options 
maximum  distance  from  wall 
stream  velocity  estimate 


output  data  file 
distance  from  wall 
mean  stream  velocity 
mean  vertical  velocity 
rms  stream  fluctuating  velocity 
rms  vertical  fluctuating  velocity 
Reynolds  stress  (velocity  product) 


Q  *  *  *  *  *  *  ft  £  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ft  ♦  ♦  ♦  ♦♦  ♦  *  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ★  ♦  ♦  ♦  ♦  ♦  ♦  *  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  *  ♦  ^  ♦  ♦  ♦  ♦ 

c  This  program  collects  tbl  data  using  an  x-probe  calibrated 
c  with  a  25  coefficient  polynomial  surface  (voltage).  No  angle 
c  correction  is  used  in  the  calculations.  A  multi-loop  procedure 
c  is  employed  where  rawdata  is  transferred  to  a  different 
c  variable  name  so  that  data  acquisition  and  calculations  can  be 
c  simultaneously  performed.  Calcuational  speed  is  improved  by 
c  using  statistical  techniques  to  determine  fluctuating  values, 
c  Note  the  references  to  gain  using  the  A/D  converter  on  the 
c  MASSCOMP.  Program  readily  adjusted  to  use  16  coefficients. 
£**********************************************•**************  ****** 
program  xtbldatacq 


135 


include  '/usr/include/mr.f 

integer  RDWR.NEARFRQ.LOW.NOTUSED 

integer  NFRAMES.NCHAN  .RANGE 

real  FREQ3RSRFQ 

parameter  (RDWR=0) 

parameter  (NEARFRQ=0) 

parameter  (FREQ=8000.0) 

parameter  (BRSFRQ=1000000.0) 

parameter  (LOW=0) 

parameter  (NOTUSED=0) 

parameter  (NFRAMES=  10000) 

parameter  (NCHAN=2) 

£*************************************************************** 
c  A/D  CONVERTER  RANGE  IS  SET  FROM  0  TO  +1 0  VOLTS  BY  THE 
c  THE  FOLLOWING  PARAMETER  STATEMENT.  TO  CHANGE,  SEE 
c  DATA  ACQUISITION  MANUAL  UNDER  "MRADRAN"  AND  RE-COMPILE. 

c*************************************************************** 

parameter  (RANGE=1 ) 
character  £rameclkdev*(*) 
character  burstclkdev*(*) 
character  addev*(*) 
cnaracter*15  fname 

parameter  (£rameclkdev=7dev/dacp0/clk0’) 
parameter  (burstclkdev=7dev/dacp0/clkl  ’) 
parameter  (addev  =  7dev/dacp0/adf0’) 
integer  adpn.clkpn  1  ,cllcpn2,fchan,incr, gain 
integer  byteslocked 
integer  status(2) 

integer* 2  rawdata(NCHAN*NFRAMES) 
integer*2  rawdatl(NCHAN*NFRAMES) 
integer*2  cgarray  (2.NCHAN) 
real  rfreq.rbrst 
common  rawdata 
common  /ufunct25/  za(25) 
common  /vfunct25/  zb(25) 
external  u25,v25 

*****************************  ***  ******  ******  ******  *********** 
c  TRAVERSE  AND  PLOTTING  INITIALIZATION  AND  ALIGNMENT 
£******************************************************* *********** 
c  traverse  control 

integer  count_clk(2)jnt_clk,clkpath0,clkdapn 
integer  intcojiumcal 
integer*4  tcount.tcount  1 
integer*4  istep 
integer*2  countervail) 
c  plotting 

character*30  lpl  ,lp2,lp3 
character*10  Inl,ln2,ln3 
real*4  xl(300),yl(300),x2(300) 
integer*4  gh(1000) 
common  /plot2/  gls.ymaxx.umaxx 
c  traverse 

common  /wncom/  counterval 

external  bwintsr 

data  clkpathO  /- 1/ 

data  clkdapn  /-l/ 

data  count_clk(  1 )  /- 1  / 

data  count_cIk(2)  /-l/ 

data  int_clk  /- 1/ 

data  tcount  /Of 

data  istep  /0/ 

data  adpn  /- 1/ 
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data  clkpnl  1-1/ 
data  clkpn2  /- 1/ 

£*************************************************************** 
c  DC  GAIN  IS  SET  TO  4  BY  THE  FOLLOWING  DATA  STATEMENT.  SEE 
c  DATA  ACQUISITION  MANUAL  UNDER  "MRADRAN". 

£***** **4^**********************%******************************* 

data  cgarray  112,2,21 
call  systemf/etc/loaddacp  O’) 
call  mrlock(0,0,byteslocked) 
write(6,*)  byteslocked,’  bytes  locked  in  memory’ 
call  mropen(clkpatb0,7dev/dacp0/clk4',0) 
call  mrclksetterfclkpathO,  15,5, 3,0, 0,2,8) 
call  mrclkarm(l,cllq)ath0) 
c  plotting 

lpl=’X-PROBE  TBL  VELOCITIES’ 
lp2=’y  (in)’ 

lp3-’uy  10*sqrt(uprms)  (m/s)’ 
lpsl=’ VELOCITY  AUTO-SPECTRUM’ 
open(8,file=’xtbldatacq.dat’,status=’fresh’) 
c* ********************************** ************************* 
c  READ  CALIBRATION  DATA  FILE 

£*  £  ♦  4C  £  ^  ^  ♦  ♦  £  ♦  ♦  frit  *  ♦  ♦  ♦  £  ft  £  ♦  £  ♦  ♦  £  ♦  *♦ ♦  ♦  ♦  ^  ft*  £  £  ♦  *  ft  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  *  ♦  * 

write(6,*)  ’Enter  name  of  (25)  coefficients  file.’ 

read(5,*)  fhame 

write(6,*) 

open(51,file=fname,status=’old’) 
do  50 1=1,25 
read(51,*)  za(i),zb(i) 

50  continue 
close(51) 

^ftftftftftftftftftftftftftftftftftftftftftftftftftftft ********************************** 

c  INITIALIZE  AND  COLLECT  RAW  DATA 
c**  • ********************************************************** 
write(6,*) 

write(6,*)  ’Enter  temperature  (F)  =  ?  ’ 

read(5,*)  tempf 

write(6,*) 

write(6,*)  ’Ensure  probe  calibration  temperature  and  test’ 
write(6,*)  ’temperature  are  within  +/-  3F.  Else  RECALIBRATE’ 
write(6,*) 

write(6,*)  ’Put  traveree  in  desired  x-location’ 

write(6,*)  ’Enter  x-location  (in).’ 

read(5,*)  xposi 

write(6,*) 

write(6,*) 

write(6,*)  ’Ensure  the  following  electrical  set-up:’ 
write(6,*)  ’ChanO=HotWirel-dc  Chanl=HotWire2-dc’ 
write(6,*) 
c  traverse 

call  mropen(count_clk(l ),’/dev/dacp0/clk6’,0) 
call  mropen(count_clk(2),7dev/dacp0/clk5  ’,0) 
call  mropen(clkdapn,  7dev/dacpO/clkda3  ’  ,0) 
call  mropen(int_clk, ’/dev/da cp0/clk7’,^) 
thresh=2.0 

call  mrclkdaset(clkdapn.thresh) 

idd  =  0 

iplot  =  1 

tcount=0 

tcountl=0 

ibuff=l 
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fchan=16 

incr=l 

c  DC  gain  is  set  by  data  statement  for  cgairay.  If  changed,  the 
c  following  statement  must  be  changed. 
gain=2 

if(gain.eq.O)  dcgain=l. 
if(gain.eq.  1)  dcgain=2. 
if(gain.eq.2)  dcgain=4. 

140  write(6,*)  ’maximum  y  measurement  ?  y=  ?’ 
read(5,*)  ymaxx 
if  (ymaxx.gt.4.8)  goto  140 
write(6,*)  ’y  -  max  =  ymaxx 
write(6,*)l 

c  If  one  desires  to  use  the  traverse  in  manual  mode  enter  a  value 
c  for  ytravmax  a  distance  above  the  wall  desired  for  initializing 
c  this  mode.  To  maintain  auto  traverse  throughout,  enter  a  value 
c  for  ytravman  >  ymaxx. 

write(6,*)  ’maximum  y  measurement  for  auto  traverse?  y=  ?’ 

read(5,*)  ytravmax 

write(6,*) 

write(6,*)  ’maximum  u  velocity  based  on  tunnel  ipm?  u=  ?’ 

read(5,*)  umaxx 

umaxx=l  .25*umaxx 

write(6  .*) 

write(6,*) 

c  Recommend  ml  be  minimum  of  10. 

write(6,*)  ’Enter  INTEGER  #  of  data  loops,  10000  FRAMES  each.’ 

read(5,*)  ml 

write(6,*) 

write(6,*)  ’Install  traverse  control  in  output  of  elk  4’ 
write(6,*)  ’Install  traverse  count  in  s  of  elk  6’ 
write(6,*)  ’Move  hot-wire  to  wall  and  set  Controller  ’ 
write(6,*>  T.  switch  to  UP  movement !!!’ 
write(6,*)  ’2.  on  computer  control’ 
wiite(6,*)  ’3.  set  y-counter  to  0’ 
write(6,*) 
write(6,*) 

write(6,*)  ’When  ready  to  take  data,  press  return.’ 

read(5,*) 

write(6,*) 

call  mropen(adpn,addev,RDWR) 
call  mropen(clkpnl,frameclkdev,RDWR) 
call  mropen(clkpn2,burstclkdev,RDWR) 
call  mradmod(adpnJRANGE,0) 
call  mradran(adpn,2,cgarray) 

c* ******************************************************** 
c  LOOP  starts  here  -  for  each  vertical  position  in  b.l. 

150  continue 

£***** ******************** ******************************** 
yposi=.001  *float(tcount) 
puav=0. 
pvav=0. 
pup2=0. 
pvp2=0. 
pupvp=0. 

c* **************************************** **************** 
cc  data  acquisition 

call  mrclk2(clkpnl,clkpn2,NOTUSEDJ^EARFRQ,FREQjfreq, 

&  NEARFRQ.BRSFRQ,ibrst,NCHAN,LOW) 
call  mradxin(adpn,clkpnl  ,clkpn2,fchan,NCHAN,incr,gain, 

&  NFR  AMES  jawdata ) 

call  mrevwt(adpn,status, 100000) 
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c* ******************************************************** 

do  405  kl=l.ml 

do  155  i=l, NCHAN*NFR AMES 
155  rawdatl(i)=rawdata(i) 

c* ****************************** ********** ****** ********** 
cc  data  acquisition 
if(kl.eq.ml)  goto  175 

call  mrclk2(clkpnl,clkpn2,NOTUSEDJ^EARFRQJ7REQjfreq, 

&  NEARFRQ3RSFRQ,rbrst,NCHAN,LOW) 

call  mradxin(  adpn,clkpn  1  ,clkpn2,fchan,NCHAN,incr,gain, 

&  NFRAMES  jawdata) 

c************************************************************** 
c  The  statistical  technique  of  calculating  fluctuating 
c  velocities  are  standard  ones  and  outlined  in  Carnahan, 
c  Applied  Numerical  Methods,  Sections  8.5  &  8.6.  Note  that 
c  rawdata  is  converted  and  stored  as  rawdatl  so  that 
c  simultaneous  data  acquisition  and  calculations  occur.  The 
c  first  data  set  is  obtained  before  entering  the  "ml"  loop, 
c  and  no  data  is  collected  in  the  loop  when  i=ml.  The  letter 
c  p  before  a  term  refers  to  cumulative  values  for  i=l  to  ml. 
c  Other  terms  are  identified  below  unless  obvious: 
c  ui  instantaneous  u  velocity  component 

c  vi=instantaneous  v  velocity  component 

c  .av=avg  value  .av2=sums  of  squares  later  averaged 
c  uvav=sum  of  instantaneous  velocity  product 
c  up2=u  prime  squared  vp2=v  prime  squared 
c  upvp=Reynolds  stress  fluctuating  velocity  product 
c************************************************************* 
175  uav=0. 
vav=0. 
uav2=0. 
vav2=0. 
uvav=0. 

do  200  i=l,NFRAMES 
j=l+NCHAN*(i-l) 

vol0=.00244 1 4 1  *real(ra  wdat  1  (j))/dcgain 

vol  1=.00244 1 4 1  *real(rawdat  1  (j+1  ))/dcgain 

ui=u25(vol0,voll) 

vi=v25(vol0,voll) 

uav=uav+ui 

vav=vav+vi 

uav2=uav2+ui*ui 

vav2=vav2+vi  * vi 

uvav=uvav+ui*vi 

200  continue 

201  uav=uav/real  (NFRAMES) 
vav=vav/real(NFRAMES) 

up2=(uav2-real(NFRAMES)*uav*  *2)/real(NFRAMES- 1 ) 
vp2=(vav2-real(NFR  AMES)*vav**2)/real(NFRAMES- 1 ) 
uvav=uvav/real(NFRAMES) 
upvp=uvav-uav*vav 
puav=puav+uav 
pvav=pvav+vav 
pupvp=^)upvp+upvp 
pup2=pup2+up2 
pvp2=pvp2+vp2 

£+  +  +  +  +  +  +  +  *  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  *+  +  +  +  +  +  +  +  +  +  1'  +  +  +  + 

if(kl.eq.ml)  goto  405 

call  mrevwt(adpn,status,  100000) 

+  +  +  *  +  +  +  +  +  +  +  +  +  +  +  +  *  +  +  +  +  +  +  +  *  +  +  +  +  +  +  *  +  +  +  +  +  +  +  *  +  +  + 

405  continue 

puav=puav/real(ml) 
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pvav=pvav/real(ml ) 
pupvp=pupvp/real(ml) 
puprms=sqrt(pup2/real(ml)) 
pvprms=sqrt(pvp2/real(ml)) 

£********> ********************************************************** 
415  x2(iplot)=10.*puprms 
xl(iplot)=puav 
yl(iplot)=yposi 

c  call  plot(iplot,xl,yl  A24pl4p2,lp3add) 
iplot=iplot+l 
idd=l 

write(8,445)  yposi,puav,pvav, 

&puprms,pupvp,pvprms 

write(6,*) 

write(6,420) 

420  format(5x,'y\7x,’uy’,8x,’vy',5x,’uprms’,5x, 

&’upvp’,5x,’vptms’) 
write(6,445)  yposi, puav.pvav, 

&puprms,pupvp,pvprms 
445  format(6f9.5,3x) 
write(6,*) 

500  write^,*)  ’tcount  =’,tcount 
write(6,*) 

^******** ********************************  ************* 
c  traverse 

if(yposi.ge.ytravmax)  goto  650 
600  call  mrclkintgat  (int_clk,l) 

call  mrcounters  (2,count_clk,3,l,0) 
call  mrbufall  (count_clk, countervail, 2, 4) 
call  m  re  Ik  arm  (2,count_clk) 
call  mrxinq  (count_clk(l),  1,1,0) 
c 

if  (tcount ge.300)  istep=3600 
if  (tcount.It.300)  istep=3000 
if  (tcount.lt.  150)  istep=2400 
if  (tcount.It.80)  istep=1200 
if  (tcount.lt.30)  istep=400 
if  (tcount.It.15)  istep=150 
if  (yposi  .gt.ymaxx)  goto  700 
call  mrclksetter(clkpath0, 15,5  ,istep, 0,0,2, 8) 
call  mrclkarm(  1  .clkpathO) 
call  astpause$(0,3000) 
c 

call  mrcIkdis(I,clkpath0) 

call  mroneshot  (int_clk,0,50.,rdelay,0,100.,rwidth,0,0) 

call  m  re  Ik  arm  (l,int_clk) 

tcount=tcount+counterval(ibuff) 

if  (ibuff.eq.l)  then 

ibuff=3 

else 

ibuff=l 

endif 

if(yposi.lt.ytravmax)  goto  690 

£*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ♦  *  *  *  *  *  *  *  *  *  *  *  ife  *  *  *  *  *  *  *  4c  *  *  *  *  *  4c  *  *  4c  *  *  **  *  *  *  * 

650  if  (yposi. gt.ymaxx)  goto  700 

write(6,*)  ’Ensure  traverse  in  manual  for  remaining  data.  ’ 
write(6,*) 

write(6,*)  ’Move  traverse  to  desired  y-position.  Enter  value’ 

write(6,*)  ’accurate  to  3  decimal  places.’ 

read  (5,*)  yman 

tcount=int( 1 000.*yman) 

write(6,*) 
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£*********************************************************** 
c  LOOP  ends  here 
690  goto  150 

g**** *********************************** ******************** 

c  finish 
c 

700  write(8,*)  tempf,’  temperature  (F)’ 

write(8,*)  xposi,’  inches  downstream  of  reference  point’ 
close(8) 

call  mroneshot  (int_clk,0,50.  ,rdelay  ,0, 1 00.  ,rwidth,0,0) 

call  mrclkarm  (l,int_clk) 

write(6,*)  ’switch  traverse  control  to  manual’ 

write(6,*)  ’press  return’ 

read(5,*) 


idd=2 

iplot=iplot-l 

call  plot(iplot,xl  ,y  1  ,x2,lp  1  ,lp24p3  ,idd) 

call  mrclkintgat(int_clk,0) 

call  mrclosall 

stop 

end 

****************************************************************** 

ccc  FUNCTIONS  TO  CONVERT  VOLTAGES  TO  VELOCITY  AND  ANGLE 
****************************************************************** 

function  u25(el,e2) 
common  Aifunct25/  a(25) 

u25=a(  1  )+e  1  *  (a(2)+e  1  *(a(3  )+e  1  *(a(4)+e  1  *a(5))))+ 

&e2*(a(6)+e  l*(a(7)+el  *(a(8)+e  1  *(a(9)+el  *a(  10))))+ 

&e2*(a(  1 1  )+el  ’*‘(a(  1 2)+e  1 *(a( 1 3)+e  1  *(a( 14)+e  1  *a(  15))))+ 

&e2*(a(  1 6)+el *(a(  1 7)+e  1  *(a(  1 8)+el  *(a(  1 9)+e  1  *a(20))))+ 

&e2*(a(2 1  )+el  *(a(22)+e  1  *(a(23)+el  *(a(24)+e  1  *a(25)))))))) 

return 

end 

£*  %  *  *  *  *  *  ★  *  *  *  *  *  +  ★  *  *  *  *  *  *  *  ★  *  *  *  *  *  *  *  *  *  *  *  *  *  4c  *  *  ★  *  *  *  *  *  *  *  ★  *  *  *  *  *  *  *  *  *  *  *  **  *  * 

function  v25(el,e2) 
common  /vfunct25/  b(25) 

v25=b(l)+e  1  *(b(2)+e  1  *(b(3)+e  1  *(b(4)+el*b(5))))+ 

&e2*(b(6)+el*(b(7)+e  1  *(b(8)+el*(b(9)+e  1  *b(  1 0))))+ 

&e2*(b(ll)+el*(b(12)+el*(b(13)+el*(b(14)+el*b(15))))+ 

&e2*(b(16)+el*(b(17)+el*(b(18)+el*(b(19)+el*b(20))))+ 

&e2*(b(21)+el*(b(22)+e  1  *(b(23)+e  1  *(b(24)+el  *b(25)))))))) 

return 

end 

£*  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  %  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  **  *  4t  *  *  *  *  *  *  *  4c  *  *  *  *  *  **  *  *  *  *  4c  *  *  * 

c  SUBROUTINES  -  PLOTTING  (see  MASSCOMP  Data  Presentation 
c  Manual) 

c  Link  plot  routines  with  sp-45  library  (-Imp) 

g* * * * * * * * * * * * * * * * * * * * * * ♦ * ** * ** * * * * * ** * * * * * * * * * * * * * * * * * * * * * ** * * 

subroutine  plot(int,xa,y  a,xb,l  1 42,13  id) 
integer*4  NULL 
integer*4  SIZEOFGCA 
integer*  4  NULLINTRAY 
real*4  UNDEF 
parameter  (NULL=0) 
parameter  (SIZEOFGCA=1000) 
parameter  (NULLINTRAY—2) 
parameter  (UNDEF=-.  1  le39) 


ccc 
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integer*  4  gls(SIZEOFGCA) 
character*  30 11,12,13 
character*  9  In 

integer*4  mkpts(2),siz(2),li(2) 

real  xa(l),ya(l),xb(l) 

common  /plot2/  gls,ymaxx,umaxx 

ln='plottbl.g’ 

call  mpinit(gls) 

call  mpgrids(gls,0,0) 

call  mpaxtype(gls,3,-l) 

call  mptics(gls,3,-l,-l) 

call  mpdevice(gls,"mcdigh",2,0) 

call  mpaxvals(gls,  1 ,0.,umaxx,UNDEF,  1 0) 

call  mpaxvals(gls,2,0.,ymaxx,UNDEF,10) 

call  mptitle(gls,4,-l,-l,’  ’) 

call  mptitle(gls,4,-l,2,ll) 

call  mptitle(gls,2,-l,-l,12) 

call  mptitle(gls, 1,-1, -1,13) 

mkpts(l)=l 

mkpts(2)=l 

li(l)=NULL 

li(2)=NULL 

call  mplines(gls,2,li,mkpts) 

siz(l)=240 

siz(2)=240 

call  mplotchrs(gls,’x*’,2,NULLINTRAY,siz) 

call  mplotsrcy  (gls,l,int,0,ya(l),T'\l,l,NUIX,NULL) 

call  mplotsrcx  (gls,l,int,0^ca(l),''F',l,129ULL,NULL) 

call  mplotsrcy  (gls,2,int,0,ya(  1  ),"F',  1 , 1  .NULL, NULL) 

call  mplotsrcx  (gls,2,int,0,xb(l),"F',l,l,NUlX,NULL) 

call  mplot(gls, 0,1,0) 

if  (id.eq.2)  then 

call  mpfile(gls,ln,l,0) 

call  mpend(gls) 

endif 

return 

end 

£*******#*********+***************#******************************* 

PROGRAM  Non-Interactive  Coles  Fitting  Algorithm  for 

AUTOX COLESFIT:  Unmanipulated  Turbulent  Boundary  Layer  Mean 

Velocity  Profile  Data 

Input  Output 

profile  data  file  name  boundary  layer  parameters  data  file 

temperature  Coles  fit  comparison  plot  file 

stream  location 

c***************************************************************** 
program  autoxcolesfit 
parameter  (npo=200,mpo=30) 
parameter  (np=5,mp=6) 
dimension  pam(mp,np),yam(mp),pr(np) 
dimension  ycor(npo),umean(npo),ufluc(npo) 
real  xplot(npo),yplot(npo),uu(npo),yy(npo) 
character*20  name.oname 
character*  20 11,12,13 
character*?  In 

common  /fu/  ycor,umean,visc,istart  jend 
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£.*  **********************  ************  ********************* 
c  program  shear 

c  reads  data-file  with  mean-velocity  and  corresponding 
c  distance  from  wall,  finds  coefficients  for  arbitrary  order 
c  polynomial  for  best  curve  fit  (usually  best  values  are: 
c  order=20  or  30  for  about  6C  data  points.) 
c  Make  sure  that  the  data  points  are  closely  spaced  . 
c 

c  for  more  info  see  Applied  Numerical  Meth., 
c  B.  Carnahan...  page  574 
c 

c  gam,del,b  coefficients  for  polynomial 
c  icount  order  of  fit 
c 

ckay  herbert,  1988 

£******************************************************** 
icount=mpo 
c==>  read  data 

write(6,*)  ’data-file-name?  ’ 
read(5,*)  name 

open  (51,file=name,status=’old’) 
i=0 

umx=0. 

ymx=0. 

10  read(51,*,end=100)  yco,um,zzzl,ur,zzz2,zzz3 
if  (um.lt.-l..or.um.gt.55.)  goto  10 
i=i+l 

ycor(i)=yco*.0254 

umean(i)=um 

ufluc(i)=ur 

if  (yco.gt.ymx)  ymx=yco 
if  (um.gtumx)  umx=um 
goto  10 
100  continue 

125  write(6,*)  ’temperature(F)  ?  ’ 
read(5,*)  tempr 

write(6,*) 

write(6,*)  ’xposition  (in)  ?  ’ 
read(5,*)  xposi 

126  close(51) 
xposi  =xposi  * . 0254 
imax=i 

m  count =i 
imax2=imax/2 

if  (ycor(imax).lt.ycoi(l))  then 
do  20  i=14max2 
ii=imax-i+l 
dum=ycor(i) 

ycor(i)=ycor(ii) 
ycor(ii)=dum 
dum=umean(i) 
umean(i)=umean(ii) 
umean(ii)=dum 
dum=ufluc(i) 
ufluc(i)=ufluc(ii) 
ufluc(ii)=dum 
20  continue 
endif 

call  param(tempr,vnu,visc,rho,cnot) 

istart=4 

iend=25 


143 


eps=0. 

deps=l.e-4 

utau=umx/20. 

dutau=utau/10. 

pr(l)=eps 

pi(2)=utau 

yam(l)=funk(pr) 

pam(l,l)=pr(l) 

pam(l,2)=pr(2) 

pr(l)=eps+deps 

pr(2)=utau 

yam(2)=funk(pr) 

pam(2,l)=pr(l) 

pam(2,2)=pr(2) 

pr(l)=eps 

prf2)=utau4dutau 

yam(3)=funk(pr) 

pam(3,l)=pr(l) 

pam(3,2)=pr(2) 

ndim=2 

ftol=l.e-6 

call  amoeba(pam,yam,mp,np,ndim,ftol,funk,iter) 

eps=pam(l,l) 

utau=pam(l,2) 
uinf=umx 
do  101  i=l,imax 
dd=(ycor(i)+eps)*utau/visc 
if  (dd.le.15)  istart=i 
101  if  (dd.le.330)  iend=i 
deps=l.e-5 
dutau=utau/30. 
pr(l)=eps 
pr(2)=utau 
yam(l)=funk(pr) 
pam(l,l)=pr(l) 
pam(l,2)=pr(2) 
pr(l)=eps+deps 
pr(2)=utau 
yam(2)=funk(pr) 
pam(2,l)=pr(l) 
pam(2,2)=pr(2) 
pr(l)=eps 
pr(2)=utau+dutau 
yam(3)=funk(pr) 
pam(3,l)=pr(l) 
pam(3,2)=pr(2) 
ndim=2 
ftol=l.e-8 

call  amoeba(pam,yam,mp,np,ndim,ftol,funk,iter) 
write(6,"')  iter,pam(l  ,1  ),pam(  1 2),umx 

c===>  plotting 
eps=pam(l,l) 
utau=pam(l,2) 
vmax=25.*utau 
vmin=2.*utau 
kz=0 

do  300  i=l,npo 

uI=vmin+(vmax-vmin)*(float(i- 1  )/float(npo- 1  ))**2 
call  spal(yl,ul,.41,5.0,utau,visc) 
if(yl.le.O.)  kz=kz+l 
yplot(i)=yl 
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xplot(i)=ul 

300  continue 
kz=kz+l 

mcount=mcount-istart 
do  301  i=l,mcount 
uu(i)=umean(i+istait)/utau 

301  yy(i)=(ycon(i+istart)+eps)/visc*utau 
do  302  i=kzjipo 

xplot(i+ 1  -kz)=xplot(i)/utau 

302  yplot(i+l-kz)=yplot(i)Msc*utau 
npp=npo+l-kz 

ll=name 

12=’u/u-tau’ 

13=’u-tau*y/visc’ 
ln=’plot2.g’ 
write(6,*)  ’done’ 

c  call  plot2(npp,yplot,xplot,mcoiint,yy, 

c  *  uujl,12,l34n,l) 

do  401  i=l,mcount 
uu(i)=ufluc(i+istart)/utau*20. 

401  yy(i)=(ycor(i+istart)+eps)/visc','utau 

ll=name 

12=’u-p*20/u-tau' 

13=’u-tau*y/visc’ 

ln=’plotl.g’ 

write(6,*)  ’done’ 

c  call  plot2(npp,yplot,xplotjncount,yy, 

c  *  uud  1, 12,13  4n,l) 

c  write(6,*)  ’output  filename  ?’ 
c  read(5,*)  oname 
oname  =’xspalinf.dat’ 
open  (51,file=oname,status=’fresh’) 
do  501  i=2,mcount+l 
uu(i  )=umean(i+istart- 1 ) 

501  yy(i)=(ycor(i+istart- 1  )+eps) 
uu(l)=0. 

yy(l)=0. 

delstar=0. 

phmom=0. 

do  502  i=l,mcount 

ddy =(yy(i+ 1  )-yy(i))/2. 

phmom=phmom+ddy*uu(i)*(uinf-uu(i)) 

phmom=phmom+ddy*uu(i)*(uinf-uu(i+l)) 

delstar=delstar+ddy*(uinf-uu(i)) 

502  delstar=delstar+ddy*(uinf-uu(i+l)) 
delstar=delstar/uinf 
phmom=phmom/(uinf*  *2) 
h=delstar/phmom 
red=uinf*delstar/visc 

cf=2.  *utau*  *2/uinf  *2 
yplus=eps*utau/visc 
twall=iho*utau*  *2 
do  503  i=l,mcount 
if(uu(i).ge.(.995*uinf))  then 
delta=yy(i) 
i=m  count 
else 

goto  503 
endif 

503  continue 
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write(51  .*)  name,'  file-name’ 

write(51,*)  tempr.’  temperature  (K)’ 

write(51,*)  xposi.’  x-position  (m)’ 

write(5 1  ,*)  vise,  ’  kin.  viscosity  (mA2/sec)’ 

write(5 1  ,* )  rho,  ’  air  density  (kg/mA3)’ 

write(51,*)eps,  ’  y-offset  (m)7  =’,yplus,’  viscous  units.’ 

write(51,*)  utau,  ’  shear  velocity  (m/sec)’ 

write(51,*)  twall,’  tauwall  (Pa)’ 

write(51  *)  delstar,’  displ.  thickness  (m)’ 

write(51,*)uinf,’  u-infinity  (m/sec)’ 

write(51*)phmom,’  momentum  thickness  (m)’ 

write(51,*)h,’  Shape  Factor’ 

write(51,*)cf,’  cF 

write(51,*)  red,’  Reynolds  number,  delstar’ 

write(51,*)  delta,’  delta  (m)’ 

close(51) 

open  (50,file=’spall.dat’,status=’old’,access=’sequential’) 

write(50,*)  uinf,utau,red,cf, delstar, twall 

close(50) 

call  system("cat  spall.dat »  cfdata.dat") 

stop 

end 

subroutine  plot2(int,xa,ya,mt'j ,  xb,yb  ,1 1 ,12  J3,lnad) 

include’/usrAnclude/libmp.F 

integer*4  gls(SIZEOFGCA) 

character*20 11,12,B 

charaaer*7  In 

integer*4  mkpts(3),siz(3) 

real  xa(  1  ),ya(  1  )^cb(  1  ),yb(  1  ),y  c(200) 

j,***********#***************************** 

c 

c  plot2,  plot  2  x-y  graphs 
c  link  with  sp-45  library  (Imp) 
c 

c  kay  herbert  1987 
c 

****************************** ****** 

call  mpinit(gls) 
call  mpgrids(gls,0,0) 
call  mpaxtype(gls,3,-l) 
call  mptics(gls,3,-l,-l) 

call  mplotsrcy  (gls,l.int,0,ya(l),"F’,l,l,NULL,NULL) 
call  mplotsrcx  (gls,Unt,0,xa(l),"F',l,l,NULL,NULL) 
call  mplotsrcy  (gls,2,intb,0,yb(l),"F",l,l,NULL,NULL) 
call  mplotsrcx  (gls,2,intb,0,xb(  1  ),"F",  1 , 1  ,NULL2SHJLL) 
mkpts(l)=0 
mkpts(2)=l 

call  mplines(gls^2,NULLINTRAYjnkpts) 

siz(l)=240 

siz(2)=240 

call  mplotchrs(gls,’*’,l,NULLINTRAY,siz) 

c===>  take  comment  out  if  screen-plot  desired 
c  call  mpdevice(gls,"mcdigb",2,0) 


if  (id.eq.l)then 
call  mplogax(gls,l,2) 
call  mpaxvals(gls,l,l.,5000.,UNDEF,6) 
call  mpaxvals(gls,2,0.,50.,UNDEF,10) 
else 
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call  mpaxvals(gls,1.0.,l.,UNDEF.10) 
call  mpaxvals(gls,2,0.,10.,UNDEF,10) 
endif 

call  mptitle(gls,4,-l,-l,’  ’) 
call  mptitle(gls,4, -1,2,11) 
call  mpiille'gls  ,2,- 1 1 ,12) 
call  mptitle(gls,  1,-1,. -1,13) 

call  mpfile(gls,ln,l,0) 
call  mplot(gls,0,l,0) 
call  mpend(gls) 
return 
end 

£*********************************************************** ****** 
subroutine  sp al(y , uJc.B.ut, vise, del tal) 
real  k,uu(100),yy(100) 
up=u/ut 
uu(l)=0.0 
yy(l)=0 
uu(2)=.99 
yy(2)=l 
uu(3)=1.96 
yy(3)=2 
uu(4)=2.90 
yy(4)=3 
uu(5)=3.80 
yy(5)=4 
uu(6)=4.65 
yy(6)=5 
uu(7)=5.45 
yy(7)=6 
uu(8)=6.19 
yy(8)-7 
uu(9)=6.87 
yy(9)=8 
uu(10)=7.49 
yy(10)=9 
uu(ll)=8.05 
yy(H)=iO 
uu(12)=9.00 
yy(12)=12 
uu(13)=9.76 
yy(13)=14 
uu(14)=10.4 
yy(14)=16 
uu(15)=10.97 
yy(15)=18 
uu(  16)=  11.49 
yy(16)=20 
uu(17)=12.34 
yy(17)=24 
uu(18)=12.99 
yy(18)=28 
uu(19)=13.48 
yy(19)=32 
uu(20)=13.r  j 
yy(20)=36 
uu(21)=14.22 
yy(21)=40 
uu(22)=14.51 
yy(22)=44 
uu(23)=14.87 
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vy.:?‘-so 

iiu(  24  >=15.33 

yy(24»=60 

uu(25)=16.04 

yy(25)=80 

uu(26)=16.60 

yy(26)=100 

uu(27)=17.61 

yy(27)=150 

uu(28)=18.33 

yy(28)=200 

uu(29)=  19.34 

yy(29)=300 

uu(30)=20.06 

yy(30)=400 

uu(31)=20.62 

yy(31)=500 

uu(32)=21.08 

yy(32)=600 

uu(33)=21.79 

yy(33)=800 

uu(34)=22.35 

yy(34)=1000 

uu(35)=23.36 

yy(35)=1500 

uu(36)=24.08 

yy(36)=2000 

uu(37)=25.09 

yy(37)=3000 

uu(38)=25.81 

yy(3g)=4000 

uu(39)=26.37 

yy(39)=5000 

uu(40)=26.83 

yy(40)=6000 

uu(41)=27.54 

yy(41)=8000 

uu(42)=28.10 

yy(42)= 10000 

i=2 

10  if  (uu(i).lt.up)  then 
i=i+l 
goto  10 
endif 

xint=(uu(i)-up)/(uu(i)-uu(i-l)) 

yp=xint*(yy(i- 1  )-yy(i))+yy(i) 

y=yp*visc/ut 

return 

end 

function  funk(pr) 
real  k 

dimension  pr(20) 

dimension  y(200),u(200) 

common  /fu/  y, u, vise, ii,iii, del tal 

k=.41 

B=5.0 

funk=0. 
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do  1  i=ii,iii 

uu=u(i) 

yl=y(i) 

call  spal(y2,uu,k  JB  ,pr(2),visc,delta  1 ) 
tunk=funk+((yl+pi<l>y2)/l.)**2 

1  continue 
return 
end 

subroutine  amoeba(p,y,mp,np,ndim,ftol,funk,iter) 

parameter  (nmax=20.alpha=1.0,beta=0.5,gamma=2.0,itmax=1500) 

dimension  p(mp,np),y(mp),pr(nmax),prr(nmax),pbarfnmax) 

mpts=ndim+l 

iter=0 

I  ilo=l 

if(y(l).gt.y(2))then 

ihi=l 

inhi=2 

else 

ihi=2 

inhi=l 

endif 

do  11  i=l,mpts 
if(y(i).lty(ilo))  ilo=i 
if(y(i).gt.y(ihi))  then 
inhi=ihi 
ihi=i 

else  if(y(i).gt.y(inhi))then 
if  (i.ne.ihi)  inhi=i 
endif 

I I  continue 

rtol=2.*abs(y(ihi)-y(ilo))/(abs(y(ihi»+abs(y(ilo))) 

if(rtol.lLftol)retum 

if(iter.eq.itmax)  then 

write(6,*)  ’exceeding  maximum  iterations’ 
return 
endif 
iter=iter+l 
do  12  j=l,ndim 
pbar(j)=0. 

12  continue 

do  14  i=l,mpts 
iffi.ne.ihi)  then 
do  13  j=l,ndim 
pbar(j)=pbar(j)+p(ij) 

13  continue 
endif 

14  continue 

do  15  j=  1  ,ndim 
pbar(j)=pbar(j)/ndim 
pr(j)=(  1  .+alpha)*pbaitj)-alpha*p(ihi,j) 

15  continue 
ypr=funk(pr) 
if(ypr.le.y(ilo))then 

do  16  j=l,ndim 

prr(j)=gamma*pr(j)+(l.-gamma)*pbar(j) 

16  continue 
yprr=funk(prr) 
if(yprr.lt.y(ilo))then 

do  17  j=l,ndim 
p(ihi,j)=prr(j) 

17  continue 
y(ihi)=yprr 
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else 

do  18  j=l,ndim 
p(ihi,j)=pr(j) 

18  continue 
y(ihi)=ypr 

endif 

else  if(ypr.ge.y(inhi))  then 
if(ypr.lt.y(ihi))then 
do  19  j=l,ndim 
p(ihi,j)=pr(j) 

19  continue 
y(ihi)=ypr 

endif 

do  21  j=l,ndim 

prr(j)=beta*p(ihi,j)+(  1  ,-beta)*pbar(j) 

21  continue 
yprr=funk(prr) 
if(yprr.lt.y(ihi)>  then 

do  22  j=l,ndim 
p(ihij)=prr(j) 

22  continue 
y(ihi)=yprr 

else 

do24i=l,mpts 
if(i.ne.ilo)then 
do  23  j=l,ndim 
pr(j)=0.5*(p(i,j)+p(ilo,j)) 

P(ij)=pr(j) 

23  continue 
y(i)=funk(pr) 

endif 

24  continue 
endif 

else 

do  25  j=l,ndim 
p(ihi,j)=pr(j) 

25  continue 
y(ihi)=ypr 

endif 
goto  1 
end 

subroutine  param(temp,visc,vkin,iho,cnot) 
c********************************************** 
c  subroutine  param.f  by  Kurt  W.  Roth,  8/88. 
c  modified  by  kay  herbert 
c 

c  routine  to  compute  viscosity  of  air  at  any  temperature 
c  Farenheit,  Celsisus,  or  Kelvin;  it  also  computes 
c  pressure  drop  for  pitot  at  any  speed  (m/s). 

c  for  viscosity  formula,  see  VISCOSITY  OF  GASES  by  WATSON; 
c  see  CRC  PHYSICS  and  CHEMISTRY  for  density  formula, 
c 

c  Error  of  formula  270-600  K  is  .25%. 
c 

c  input  temperature  in  fahrenheit(F) 
c  output  viscosity(kg/m-s) 
c  kinematic  viscosity(mA2/s) 
c  density(kg/mA3) 
c  speed  of  sound(m/s) 
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c  temperature  (K) 
c 

£*********************************************** 
integer  icaler,krunch 

real  TEMP,  CONVER,  VISC,RE,RHO,VKIN,VEL,DP,CON 
icaler=l 


if  (icaler  .eq.  2 )  goto  2 
if  (icaler  .eq.  3)  goto  1 
TEMP  =(  TEMP  -  32)/1.8 
if  (icaler  .eq.  1 )  goto  2 

1  TEMP  =  TEMP  -273.15 

2  RHO  =  1 .293/(  1  +  ,00367*TEMP) 
TEMP  =  TEMP +  273.15 


s  =100/(TEMP) 
aO  =  .55279544 
al=  2.818916 
a2  =  -13.508340 
a3  =  39.353086 
a4  = -41.419387 
RE  *  sqrt(TEMP) 

VISC=(RE)/(aO  +al*s  +a2*s**2  +a3*s**3  +a4*s**4)/10**6 

VKIN  =  VISC/RHO 

cnot=l  .4*287*temp 

cnot=sqrt(cnot) 

return 

end 

c*************************** ************************************** 

PROGRAM  Interactive  Coles  Fitting  Algorithm 

MANLIALXCOLE: 


Input  Output 

profile  data  file  name  Coles  fit  comparison  plot  file 

Coles  data  file  Values  for  ux  and  y+  offset 

temperature 
ux  —  estimate 
y*  offset  —  estimate 

g***  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ★  *  **  %  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  *  *  *  *  *  *  * 

c  This  program  uses  historical  data  from  Coles  for  a  turbulent 
c  boundary  layer  mean  velocity  profile  and  compares  the 
c  experimental  mean  velocity  profile  to  that  of  Coles.  After 
c  enterin  appropriate  values  for  utau  and  y-offset,  a  visual 
c  fit  to  Coles  can  be  performed  (Note:  utau  may  be  obtained 
c  from  surface  fence,  preston  tube,  or  estimated.).  The  fit 
c  should  be  performed  such  that  data  for  y-plus  of  20-200 
c  either  lies  directly  on  Coles  data  (in  the  case  where  no 
c  determination  of  utau  was  made)  or  parallel  to  the  same  (if 
c  utau  is  otherwise  determined), 
c  -  frank  camelio,  june  1989 

£  ft  %  *  £  *  ft  *  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ft  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ft  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  £  *  *  Jfr  *  £  £  *  ♦  *  *  *  £  £  *  * 

program  manualxcole 

real  xu(300),xv(300),y(300),ya(50),ua(50) 

real  xuu(300),xuv(300),xw(300),xup(300),yp(300) 

real  temp,roh,mu,nu 

character*  30 11, 12,13,  file  1 

character*  15  lnl 
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12='  UMEAN/UTAU ' 

13='YPLUS  (y*utau/nu)’ 
lnl=xtunu.g’ 

write(6,*)  'Enter  file  name  containing  profile  data.’ 

read(5,*)  filel 

write(6,*) 

ll=filel 

i=l 

open(5 1  ,file=file  1  .statu  s=  ’  old  ’ ) 

50  read(51,*,end=100)  y(i),xu(i),xv(i),xuu(i),xuv(i),xvv(i) 
i=i+l 
goto  50 
100  continue 
close(51) 
int=i-l 
k=l 

open(52,file=’coles.dat\status=’old’) 

125  read(52,*,end=130)  ua(k),ya(k) 
k=k+l 
goto  125 
130  continue 
close(52) 
inta=k-l 

c*******  Reassign  values  to  select  data  for  plotting  *************** 
write(6,*)  ’Enter  temperature  (deg  F)  at  time  data  collected.’ 
read(5,*)  temp 
temp=(temp-32.)/l  .8 
roh=1.293/(l+.00367*temp) 
mu=1.853e-05-(.005e-05*(27-temp)) 
nu=mu/roh 

160  write(6,*)  ’Enter  value  for  utau  (m/sec).’ 
read(5,*)  utau 
write(6,*) 

write^,*)  ’Enter  offset  in  y-plus  (viscous  units).’ 

read(5,*)  yoffset 

write(6,*) 

yoff=yoffset*nu/utau 

jp=l 

do  200  j=l,int 
xupfj)=xu(j)/utau 
ypp=.0254*y(j)*utau/nu+yoffset 
if(ypp.lt.l.)  goto  200 

yp(jp)=ypp 

xup(jp)=xup(j) 

jP=JP+l 
200  continue 
jp=jp-l 

call  plot2(inta,ua,ya,jp,xup,yp,l  l,12,13Jnl) 
write(6,*)  ’utau  =  \utau,’  m/sec’ 

write(6,*)  ’yoffset  (m)=  ’.yoff,’  =  ’.yoffset,’  viscous  units.’ 

250  write(6,*) 

write(6,*)  ’Enter  INTEGER  1  to  try  a  different’ 
write(6,*)  ’y-offset  and/or  utau.  Else  any  integer.’ 
read(5,*)kl 
if(kl.eq.l)goto  160 
1000  stop 
end 


£**  *  *  4c  **  *  *  **  *  *  *  *  *  *  *  *  *  *  *  *  **  **  *  *  *  *  *  *  *  *  **  *  *  *  ***  *  *  *  4c  **  *  **  *  **  *  ***  **  4c  *£ 


subroutine  plot2(int,ua,ya,intp,xup,yp,I  1 , 12,13 ,lnl ) 
include’/usr/include/libmp.f 
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integerM  gls(SIZEOFGCA) 
character*30 11,12,13 
character*  10  lnl 

integerM  mkpts(2).siz(2),LINRAY(2) 

real  ua(  l),ya(  1  ),xup(l  ),yp(  1 ) 

call  mpinit(gls) 

call  mpgrids(gls,0,0) 

call  mpaxtype(gls,3,-l) 

call  mptics(gls,3,-l,-l) 

call  mplotsrcy  (gls ,  1  ant,0,ua(  1  ),"F' ,  1 , 1 ,NULL,NULL) 

call  mplotsrcx  (gls ,  1 ,int,0,y  a(  1  ),”F’ ,  1 ,1 ,NULL,NULL) 

call  mplotsrcy  (gls,2,intp.0,xup(l),"F',l,l  .NULL, NULL) 

call  mplotsrcx  (gls^,intp,0,yp(l),"F", 1,1, NULL J^TJLL) 

mkpts(l)=0 

mkpts(2)=l 

LINRAY(1)=1 

LINRAY(2)=1 

call  mplines(gls,2,LINRA  Y  jnkpts) 

siz(l)=200 

siz(2)=200 

call  mplotchrs(gls,’*’,l,NULLINTRAY,siz) 

call  mpdevice(gls,"mcdigh",2,0) 
call  mplogax(gls,l,2) 

caU  mpaxvals(gls,  1  ,UNDEF,5000.  ,UNDEF,6) 

call  mpaxvals(gls,2,0.,40.,UNDEF,6) 

call  mptitle(gls,4,-l,-l,’  ’) 

call  mptitle(gls,4, -1,2,11) 

call  mptitle(gls,2,-l,-l,12) 

call  mptitle(gls,l,-l,-l,13) 

call  mpfile(gls,lnl,l,0) 

call  mplot(gls, 0,1,0) 

call  mpend(gls) 

return 

end 

c**************************************************************** 

DATA  FILE:  Data  of  Coles  [3]  to  Which  Experimental  Mean 

Velocity  Data  is  Fit  (Preceding  Program) 

*  *  4c  *  *  **  *  *  *  *  *  *  **  *  *  *  *  *  *  **  *  ★  *  *  *  *  %  *  *  *  *  *  *  *  4c  *  *  *  %  *  *  *  *  **  *  **  *  *  *  *  *  *  *  *  ★  *  *  *  * 

.991. 

I. 96  2. 

2.9  3. 

3.8  4. 

4.65  5. 

5.45  6. 

6.19  7. 

6.87  8. 

7.49  9. 

8.05  10. 

9. 12. 

9.76  14. 

10.4  16. 

10.97  18. 

II. 49  20. 

12.34  24. 

12.99  28. 

13.48  32. 

13.88  36. 

14.22  40. 

14.51  44. 
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14.87  50. 

15.33  60. 
16.04  80. 
16.6  100. 

17.61  150. 

18.33  200. 

19.34  300. 
20.06400. 

20.62  500. 
21.08  600. 
21.79  800. 

22.35  1000. 

23.36  1500. 
24.08  2000. 
25.09  3000. 
25.81  4000. 

26.37  5000. 
26.83  6000. 
27.54  8000. 
28.1  10000. 
100.  10001. 


.***************************************************************** 


PROGRAM 

XTBLPARAMETERS: 


Coles  Interactive  Fitting  Algorithm  and  Calculation  of 
Various  Boundary  Layer  Parameters 


Input 


Output 


profile  data  file  name  boundary  layer  parameters  data  file 

temperature 
uz  (Coles  fit) 
y*  offset  (Coles  fit) 


c************************************************************** 
c  This  program  reads  boundary  layer  data  as  measured 
c  with  an  X-probe  (reads  output  of  data  acquisition  program 
c  such  as  xtbldatacq.f).  The  values  of  utau  and  y-offset 
c  from  either  a  manual  or  auto  Coles  fit  program  are 
c  required  to  calculate  certain  parameters.  Note  that  the 
c  boundary  layer  thickness,  delta,  is  obtained  approximately 
c  using  the  .995  definition.  This  was  chosen  so  that  data 
c  can  be  compared  with  historical  results  (Klebanoff)  when 
c  miscellaneous  plots  are  obtained.  Another  option  would  be 
c  to  use  the  law  of  the  wake  and  incorporate  determination  of 
c  delta  in  Coles  fit  routine. 

£*  £  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦♦  ★  ♦  ♦  *  £  ft  %  ft  ♦  ♦  ♦  £  ♦  ♦  ft  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦  ♦♦  ♦  ♦  *  £  ft 

program  xtblparameters 

real  y(300),u(300),v(300),uu(300),uv(300),w(300) 
real  del(300),theta(300),nu 
character*  20  filel 

open(7  ,file=  ’xtblanal  ,dat  ’  ,status=  ’fresh  ’) 
write(6,*)  ’Enter  tbl  data  file  name.’ 
read(5,*)  filel 
write(6,*) 

write(6  *)  ’Enter  value  for  y-offset  based  on  Coles  fit’ 

read(5,*)  yoff 

write(6,*) 

write(6,*)  ’Enter  value  for  utau  based  on  Coles  fit’ 

read(5,*)  utau 

write(6,*) 
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write(6,*)  ’Enter  temp  (deg  F)  at  time  of  b.l.  measurements.’ 

read(5,*)temp 

write(6,*) 

temp=(temp-32.)/l  .8 

iho=1.293/(l.+.00367*temp) 

nu=(1.853e-05-.005e-05*(27-temp))/ifao 

yO=nu*yoff/utau 

y(l)=0. 

u(l)=0. 

v(l)=0. 

uu(l)=0. 

uv(l)=0. 

w(l)=0. 

i=2 

open(5 1  ,file=file  1  ,status=  ’old  ’ ) 

100  read(51,*,end=200)  y(i),u(i),v(i),uu(i),uv(i),w(i) 
y(i)=.0254*y(i)+y0 
uinl=max(uinf,u(i)) 
i=i+l 
goto  100 
200  close(51) 
n=i-l 

do  250  i=l,n 

if(u(i)/uinf.ge.0.995)  then 
delta=y(i) 
i=n+l 
goto  250 
endif 

250  continue 
do300i=l,n 
del(i)=l-u(i)/uinf 
il(del(i).lt0)del(i)=0. 
theta(i)=del(i)*u(i)/uinf 
300  if(theta(i).lt.0)  theta(i)=0. 

delstar=(  1  -del(2))*yO*.5-Klel(2)*y(2) 
thestai=theta(2)*y(2)*  .5 
do  400  i=2,n-l 

delstan=delstar+0.5*(del(i+l)+del(i))*(y(i+l)-y(i)) 

400  thestar=thestar+0.5*(theta(i+l)+theta(i))*(y(i+l)-y(i)) 
h=delstar/thestar 
Redel=uinf*'delstar/nu 
01=2*0130**2/0101**2 
tauwall=ifio  *utau  *  *2 

write(6,*)  ’displacement  thickness  (m)  =  ’,  delstar 

write(6,*)  ’momentum  thickness  (m)  =  ’,  thestar 

write(6,*)  ’shape  lactor  =  ’,  h 

write(6,*)  ’yoflset  (m)  =’,y0,’  =  \yofl,’  viscous  units’ 

write(6,*)  ’kinematic  viscosity  (mA2/sec)  =’4iu 

write(6,*)  ’air  density  (kg/mA3)  =  ’,rtio 

write(6,*)  ’uinl  (m/sec)  =’,uinl 

write(6,*)  ’temp  C  =  ’.temp 

write(6,*)  ’Reynolds  #  (delstar)  =  ’.Redel 

write(6,*)  ’delta  (m)  -’.delta 

write(6,*)  ’cl  =  ’,  cl 

write(6,*)  ’tauwall  (Pa)  =  ’.tauwall 

write(6,*)  ’friction  velocity  (m/sec)  =  ’,  utau 

write(7,*)  ’displacement  thickness  (m)  =  ’,  delstar 

write(7,*)  ’momentum  thickness  (m)  =  ’,  thestar 

write(7,*)  ’shape  lactor  =  ’,  h 

write(7,*)  ’yollset  (m)  =’,y0,’  =  ’.yoll,’  viscous  units’ 

write(7,*)  ’kinematic  viscosity  (mA2/sec)  =’nu 

write(7,*)  ’air  density  (kg/mA3)  =  \iho 
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write(7,*)  'uinf  (m/sec)  ='.uinf 

write(7,*)  ’temp  C  =  ’.temp 

write(7,*)  ’Reynolds  #  (delstar)  =  ’.Redel 

write(7,*)  ’delta  (m)  ~ ’.delta 

write(7,*)  ’cf  =  ’,  cf 

write(7,*)  ’tauwall  (Pa)  =  ’.tauwall 

write(7.*)  ’friction  velocity  (m/sec)  =  ’,  utau 

close(7) 

stop 

end 

c***********************m«*********************************** 

PROGRAM  XTBLPLOT:  Plotting  of  Boundaiy  Layer  Parameters 


Output 


Input 

profile  data  file  name 
Klebanoff  data  file 
ux  (Coles) 

y  offset  in  meters  (Coles) 


five  plot  files 
mean  velocity  profiles 
Klebanoff  comparison  profiles 
stream  rms  fluctuating 
vertical  rms  fluctuating 
Reynolds  stress  (normalized) 
shear  correlation 


c********************************************************** 
c  This  program  generates  five  plots  of  mean  velocities, 
c  fluctuating  velocities,  Reynolds  stress,  and  shear 
c  correlation  as  a  function  of  distance  or  non-dimensional 
c  distance  from  the  wall.  Fluctuating  velocities  are 
c  plotted  with  historical  data  of  Klebanoff.  Input  data 
c  is  obtained  from  output  of  data  acquisition  program 
c  purged  of  non-numeric  terms,  from  Klebanoff  data  file, 
c  and  interactively  from  results  of  Coles  fitting  routine 
c  (utau,  y-offset). 

c  -  frank  camelio,  june  1989 

c********************************************************** 
program  xtblplot 

real  xu(300),xv(300),y(300),yd(300) 
real  xuu(300),xuv(300),xw(300),xut(300) 
real  vk(50),ykv(50).uk(50),yku(50) 

character*30 11,112,113422,123,132,141,142,151452 

character*?  Inl,ln2,ln34n4,ln5 

character*  15  fname 

112=’  HEIGHT  ABOVE  WALL  (IN.)’ 

113=’(*)UAVG  (x)VAVG  (M/SEC)’ 

122=’(*)u/U  (x)u/U  (Klebanoff)’ 

123=’  y/delta’ 

132=’(*)v/U  (x)v/U  (Klebanoff)’ 

141=’NORMALIZED  REYNOLDS  STRESS’ 

142=’  -uv/utauA2’ 

151= ’SHEAR  CORRELATION’ 

152=’  -uv/(uims*vrms)’ 
lnl=’xtunl.g’ 

In2=’xtun2.g’ 

In3=’xtun3.g’ 

In4=’xtun4.g’ 

In5=’xtun5.g’ 
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open(52.fiIe=’kJebnoff.dat',status='old’) 
read(52,*)  nku,nkv 
do  25  i=l,nku 

25  read(52,*)  uk(i),yku(i) 

do  50  i=l,nkv 

50  read(52,*)  vk(i),ykv(i) 

close(52) 

write(6,*)  ’Enter  name  of  tbl  complete  velocity  data  file.  ’ 

read(5,*)  fname 

ll=faame 

write(6,*) 

write^,*)  ’Enter  utau  in  m/sec  based  on  Coles  fit.’ 
read(5,*)  utau 

write(6,*)  _ 

write(6,*)  ’Enter  yoffset  in  METERS  determined  from  analysis  or’ 

write(6,*)  ’an  impropriate  fitting  routine.’ 

read(5,*)  yoff 

write(6,*) 

i=l 

open(51,file=friame,status=’old’) 

70  read(51,*,end=100)  y(i),xu(i),xv(i),xuu(i),xuv(i),xvv(i) 
c  data  in  y  originally  in  inches  converted  to  meters. 
y(i)=.0254*y(i)+yoff 
i=i+l 
goto  70 
100  continue 
close(51) 
int=i-l 
uinf=0. 
do  150j=l,int 
uinf=max(uinf,xu(j)) 

150  continue 

c  Note  the  .995  method  of  calcualdng  delta 
do  155  i=l,int 
if(xu(i)/uinf.ge.0.995)  then 
delta=y(i) 
i=int+l 
endif 

155  continue 
do  160  j=l,int 
xut(j)=-xuv(j)/(utau**2) 

160  continue 
do  170j=l,int 

xuv(j)=-xuv(j)/(xuu(j)*xw(j)) 

170  continue 
do  200  j=l,int 
yd(j)=y<j)/delta 
xuu(j)=xuu(j)/uinf 
200  continue 
do  400j=l,int 
xw(j)=xw(j)/uinf 
400  continue 

£**** *************************************************** ****** 
call  plot  l(int,y ,xu,int,y,xv,l  1 ,1 1 2,1 1 3, In  1 , 1 ) 
call  plot  1  (int ,xuu,yd,nku ,uk,yku,l  1 ,122,123,ln2,2) 
call  plot  1  (int,x  w,yd,nkv,vk,ykv,ll ,132,123  ,ln3,3) 
call  plot2(int,xut,yd,14 1 ,142,123, ln4,4) 
call  plot2(int,xuv,yd,151,152,1234n5,5) 


! 
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stop 

end 


£***************************************************************(; 


subroutine  plotl(intl,yl,xl,int2,y2ot2411, 112413 JnJm) 

include  Vusr/include/libmp.f 

integer*4  gls(SIZEOFGCA) 

integer*4  linray(2),mikray(2),siz(2) 

real  x  1  (int  1  ),y  1  (int  1  ),x2(int2),y2(int2) 

character*  30 111412,113 

character*?  In 

call  mpinit(gls) 

call  mpgrids(gls,0,0) 

call  mpaxtype(gls,3,-l) 

call  mplotsrcy  (gls,14ntl,0,y  1(1  ),"F", 1,1, NULL J^IJLL) 

call  mplotsrcx  (gls,14ntl,0,xl(l),,,F,\l,l,NULL4sJULL) 

call  mplotsrcy  (glsr24nt2,0,y2(l),"F',l,l,NULL4^ULL) 

call  mplotsrcx  (gls,2,int2,0,x2(  1  ),"F",1 ,1  ,NULL,NULL) 

call  mpdevice(gls,"mcdigh",2,0) 

linray(l)=l 

linray(2)=l 

mrkray(l)=l 

mikray(2)=l 

call  mplines  (g!s,2,LINRA  Y  ,MRKRA  Y) 

siz(l)=150 

siz(2)=180 

call  mplotchrs(gls,’*x’^,NULLINTRAY,siz) 

call  mpaxvals(gls,  1  ,UNDEF,UNDEF,UNDEF,6) 

callmptitle(gls,4,-l,-l,’  ’) 

call  mptitle(gls,4, -1,2,111) 

call  mptitle(gls,2,- 1,-1 ,112) 

call  mptitle(gls,l,-l, -1,113) 

call  mpfile(gls,ln,l,0) 

write(6 ,*)  ’To  view  plot  ’,kn,’  type  ’Jen,’  and  (return)’ 
read(5,*)  kl 
if(kl.eqJcn)  then 

call  mplot(gls,0,l,0) 
else 

goto  778 
endif 

778  continue 
call  mpend(gls) 
return 
end 


subroutine  plot2(int,xa,y,ll  1  ,U2,113,ln^kJ) 

include’/usr/include/libmp.F 

integer*4  gls(SI2£OFGCA) 

integer*4  linray(  1  )jnrkray  ( 1  ),siz(  1 ) 

real  xa(int),y(int) 

character*30 111412,113 

character*?  In 

call  mpinit(gls) 

call  mpgrids(gls,0,0; 

call  mpaxtype(gls3.-l) 

call  mplotsrcy  (gls,14nt,0^a(l),"F',l,14^ULL,NULL) 
call  mplotsrcx  (gls,  1 4nt,0,y(l  ),"F',1 , 1  ,NULL,NULL) 
call  mpdevice(gls,”mcdigh",2,0) 

Iinray(l)=l 

mikray(l)=l 

call  mplines  (g!s,l,LINRAY,MRKRAY) 
siz(l)=200 

call  mplotchrs(gls,’*’,l,NULLINTRAY,siz) 
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if(kl.eq.4)  call  mpaxvals(gls.l,UNDEF.LTNDEF,UNDEF,6) 

if(kl.eq.5)  call  mpaxvals(gls,2,UNDEF,1.0,UNDEF,6) 

czdl  mplitie(gls,4,-l,-l,’  ’) 

call  mptitle(gls  ,4,- 1,2,111) 

call  mptitle(gls  ,2, -1,-1 ,112) 

call  mptitle(gls, 1,-1, -1,113) 

call  mpfile(gls,ln,l,0) 

write(6,*)  ’To  view  plot  \kl,’  type  \kl,’  and  (return)’ 
read(5,*)k2 
if  (k2.eq.kl)  then 
call  mplot(gls, 0,1,0) 
else 

goto  779 
endif 

779  continue 
call  mpend(gls) 
return 
end 

c***************************************************************c 
DATA  FILE:  Klebanoff  (TN  3178)  Historical  Data  of  Fluctuating 

Stream  and  Vertical  Velocity  Components 

C  *  *  *  *  *  *  *  *  4c  *  4c  *  *  *  *  *  *  *  *  *  4c  4c  *  %*  *  *  *  *  4c  *  *  *  *  *  *  *  **  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  4c  *  *  * 

33  20 

0.  0, 

.061  .0015 
.075  .00188 
.084  .00225 
.087  .00275 
.092  .00313 
.103  .00413 
.107  .005 
.110  .00675 
.112  .00838 
.103  .0146 
.0865  .0213 
.0795  .040 
.0748  .083 
.071  .133 
.0715  .143 
.0643  .230 
.0625  .258 
.061  .295 
.0605  .330 
.055  .405 
.0535  .430 
.050  .493 
.0508  .495 
.0425  .585 
.040  .6175 
.0383  .655 
.0369  .655 
.0289  .785 
.013  .89 
.0065  1.0 
.0035  1.15 
.0015  1.33 
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0.  0. 

.0325  .017 
.033  .020 
.034  .023 
.036  .040 
.039  .103 
.0394  .23 
.0384  .30 
.0345  .39 
.033  .49 
.031  .506 
.032  .527 
.0226  .673 
.024  .693 
.0226  .723 
.0186  .793 
.0156  .84 
.0078  1.003 
.0034  1.17 
.0020  1.33 

£********************************************  V******************** 
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